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EFFECTS OF SIRAS -INDUCED OPTICAL .jCIIVin AND 
DEAMINATION OF OPTICALLY SOOI/ILF'E .iYDELS 
IN THRES-DIMENS IONAL PHOPOSLaSTICIIY 


A Thesis Submitted 

In Partial Fulfilment of the Requirements 
for the Degree of 

Ph.D. 

by 

A.V.S.S.S.R.SARMA 

DEPARTMENT OF MECHANICAL ENGINEER DIG 

INDIAN INSTITUTE OF TECHNOLOGY, IAN PUR 

JULY A972 

Among the methods available to determine the complete 
state of stress at any point in a three-dimensional body, photo 
elasticity ranks very high. In photoelasticity, there are 
essentially two methods, (a) using transmitted, light, and 
(b) using scattered light. The transmitted light technique 
presently needs a frozen stress model which has to be cut into 
slices and sub-slices to obtain the required data. Hence, the 
method is destructive in nature. The scattered light technique 
however is nondestructive and has certain other advantages in 
the sense that no frozen stress model is necessary. The models 
can be subjected to live loads at room temperatures and 
analysed. 


In scattered light technique there are essentially 
two methods, one developed by Frocht and Srinath, and the other 
developed by Robert and Guillemet. The method developed by 
Frocht and Srinath is limited to cases where the effects of 
stress -induced optical activity is not high. This is equi- 
valent to saying that the rate of rotation of the secondary 
principal-stresses along the light path is small compared to 
the rate of retardation. The method developed by Robert and 
Guillemet does not involve any assumption regarding optical 
activity, but requires certain characteristic parameters which 
are indirectly related to the secondary principal-stresses and 
their directions. 

vjy' 

The present investigation is concerned with the 
following topics: 

1-) Development of relevant optical equations using Stokes 

parameters for the analysis of problems involving optical 
activity, 

2) Experimental techniques to determine the characteristic 
parameters, that is the primary and secondary characteristic 
directions and the characteristic retardation; 

3) Reduction of a three-dimensional photoelastic model to an 
optically equivalent system; 

*+ ) New methods for scattered light analysis without involving 

» 

assumptions regarding stress-induced optical activity; 



5) Development of new graphical methods to determine the 
integrated effect in problems involving optical activity; 

6) Development of a large-field compensator using a few 
plates of arbitrary but equal retardation -an extension of 
Pancharatnam ' s andTuzi’s work; 

7) Determination of integral fringe orders using different 
wavelengths of light; 

8) Experimental verification of the above methods using a 
model under combined torsion and tension; 

A short write-up on the above topics is giveb below: 

The usual photoelastic analysis consists of trigono- 
metrical resolutions of the light vector at various stages, 
as it passes through the model* However, the study of light 
ellipse and its properties as it passes through the various 
optical elements provides a better and elegant method of 
analysis. This concept has been used throughout this inves- 
tigation. To start with, various previous systems of equation 
of three -dimen si onal photbelasticity are studied. Then these 
optical equations are kinematically derived in terms of Stokes 
parameters. Since the transformations of these parameters 
when the light ellipse passes through a retarder are linear 
in these parameters, the resulting optical equations derived 
in terms of Stokes parameters turn out to he a set of ordinary 
linear differential equatiohs which are very convenient for 



analysis and application. Using these equations, the analytical 
solution for the case where the ratio of the rate of rotation 
to the rate of retardation is a constant, is obtained (The rate 
of rotation itself can be any arbitrary function). This case is 
more general compared to the case investigated by Drucker and 
Mindlin where each of these rates is a constant. 

From these equations, the equivalence of any optical 
path in a three-dimensional photoelastic model to a system 
consisting of a retarder and a pure rotator follows directly. 

Two new direct methods of finding these characteristic para- 
meters experimentally in a transmission polariscope are presented. 
Extension of those methods to the scattered light technique is 
described. These methods of scattered light analysis do not 
involve any assumptions regarding optical activity. Also 
these do not require the rotation of the analyser at a parti- 
cular speed as is required in the method of Robert and 
Guillemet. A direct and whole field method which is derived 
from the optical equations in terms of Stokes parameters is 
presented. This method again is general and it does not in- 
volve any assumptions. 

Further some investigations are carried out on a few 
related topics. The graphical methods are very convenient to 
find the integrated optical effect when the secondary principal- 
stress difference and their directions are given. The existing 
graphical methods are essentially based on stereographic or 
parallel projection of the Poincare sphere or the quaternions. 


A brief review of these methods is presented. A more direct 
and simpler approach to develop the graphical methods which 
is based on the properties of the light ellipse, is then 
presented. Two simple graphical methods which have certain 
advantages over the previous methods are described. 

The compensator is an indispensable optical element 
in the photoelastic work. Construction of a large field com- 
pensator using a few plates of arbitary hut same retardation 
is described. This is essentially an extension of the works 
of Pancharatnam and Tuzi. An extension of Pant’s method to 
find the integral fringe order using the fractional fringe 
order for different wavelengths of light, is presented. 

Some experimental results obtained from a rectangular 
photoelastic model under combined tension and torsion in a 
transmission polariscope are presented as supporting evidence 
for the existence of an optically equivalent system, the two 
methods of determining these parameters experimentally, limit a 
tions of the previous scattered light methods and the determi- 
nation of the integral frihge orders. 



CHAPTER 1 


INTRODUCT I ON 

Phot oelastieity and Electromagnetic Theory of Light 

Photoelasticity deals with the study of propaga- 
tion of light through elastic bodies under load in order 
to determine the stress distribution inside these bodies. 
Photoplasticity, photoviscoelasticity, dynamic photoelas- 
ticity etc. are immediate extensions of the basic princi- 
ples of photoelasticity. The basic phenomenon of double 
refraction on which photoelasticity stands was discovered 
in 1816 by David Brewster. Neumann formulated the stress- 
- optic law in terms of strain in 1 84-1 . later after ten 
years, the same law was formulated in terms of stress by 
Maxwell. The first engineering applications were tried 
by French engineer Mesnager and the British scientists 
Coker and Filon in the beginning of the present century. 

The electromagnetic theory of light explains all 
the photoelastic phenomena. According to this theory' 1 , 

it is assumed that the transmission of all forms of radiant 

/ 

energy (including light) through space takes place in the 
form of electromagnetic disturbances. The radiant energy 
of light is equally divided between electric energy and 
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magnetic energy. At any instant, these electric and. mag- 
netic fields at a given point in a ray of light may he repre- 
sented by two equal and mutually perpendicular vectors called 
electric and magnetic vectors in a plane perpendicular to 
the ray of light . Because of the interdependence of these 
two vectors, we can for all practical purposes in the photo- 
elastic analysis, consider only one vector say, the electric 
vector. The electromagnetic spectrum of radiant energy 
covers in addition to the visible light rays, the heat rays, 
X-rays, radio waves and the shorter gamma rays. In ordinary 
or unpolarized light, the light vector (electric vector) has 
no directional preference and its changes are chaotic. When 
the electric' vector changes its azimuth (orientation with 
respect to some fixed axes) and magnitude according to a 
certain law, we have polarized light. When the tip of the 
electric vector traces out an ellipse or a circle, it is 
elliptically polarized light or circularly polarized light 
respectively. When the azimuth of light vector does not 
change, that is, the electric vector stays in a particular 
plane, it is plane polarized light. This phenomenon of 
polarization is observed throughout the electromagnetic 
spectrum. 

Stress-Optic Law in Two and Three -Dimensions 

In a transparent isotropic plate, in which the 
stresses are two-dimensional, and within the elastic limit, 
the phase difference or relative retardation p in wavelengths 
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between the rectangular wave components of light travelling 
through it and produced by temporary birefringence is given 
by 

p = C A h (p - q) (1 .1 ) 

in which is a constant known as the stress- 

optic coefficient for the particular material for wavelength 
X , h is the thickness of the plate and p and q the prin- 
cipal-stresses. This temporary retardation or phase diffe- 
rence (this vanishes when the loads are removed) is intro- 
duced along the principal-stress axes. 

Stress-Optic Law in Three-Dimensions 

In a three-dimensional photoelastic model, a retar- 
dation (in wavelengths or fringes) equal to 

dn = dh (p T - q' ) (1 .2) 

will be added along the secondary principal-stress 
axes to the phase of the components of light vector along 
these directions. is the stress-optic coefficient as 

in two-dimensional case, dh is the actual length of light 
path (dh is so small that the rotation of the secondary 
principal-stress axes over this distance is negligible), p' , 
q' are the secondary principal-stresses for the direction 
of the given ray. The secondary principal-stresses for a 
given direction of light propagation are defined as the 
principal-stresses resulting from the stress components 
which lie in a plane normal to the direction of light ray. 

In both the above cases, the stress-optic coefficient may 




also depend upon the temperature of the model material. 

Thus, the photoelastic analysis essentially involves 
the measurement of retardation introduced in the model and 
finding the directions along which this retardation is intro- 
duced. A photoelastic model under an isotropic state of 
stress will not introduce any retardation and a general 
state of stress can he broken into two parts, the isotropic 
part and the deviatoric part. So, the photoelastic analysis 

will yield five stress components at a point out of the total 
o 

six components . The other unknown can be found by using 

either the shear difference method ^ , oblique incidence^ , 

2 S 6 

electric analogy or by an interferometric method 
The photoelastic technique in the two dimensional case is 
almost perfected. In this case, the model making and the 
simulation of the loading of the prototype take much of the 
time while the experiment is direct and simple. The bire- 
fringent coating method and photoelastic study with a ref- 
lection polarise ope eliminate model making without loosing 
the advantages of the transmission method^ . 

Three-Dimensional Phot oelastic Techniques and Rotational 
Effects; 

In the general three-dimensional cases of stress 
distribution, the secondary principal -stresses vary conti- 
nuously along the wavenormal and this makes the analysis 
difficult. The three-dimensional photoelastic techniques 
in vogue are essentially two, (i) the stress freezing 
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technique where the stresses are frozen into the model 
and the model is carefully cut into very thin slices without 
disturbing the frozen stress pattern, so that in each slice 
the rotational effect is negligible (to freeze the stresses 
into it, the model is made to undergo a special thermal 
cycle) and (li) the scattered light methods where 

measurements on light scattered in a plane perpendicular 

to the light path are made. 

17 

In 1939 Weller made use of the scattered light 

phenomenon in his analysis of three-dimensional state of 

stress in a photoelastic model. Using the ©ther theory of 

18 

light Drucker and Mindlin investigated the effect of 
rotation of secondary principal-stress axes along the wave 
normal, on the retardation in a photoelastic model. They 
solved the problem where the principal-stress magnitude 
remains constant and their directions change linearly along 
the wavenormal. They concluded that (i) the rotation incre- 
ased the relative retardation by a factor S = + *+R 2 )^ 

where R is the ratio of the rate of rotation to the rate of 
retardation and that (ii) at high stress the light vector 
rotated with the principal-stress directions, plane polarized 
light entering along a direction of principal-stress remained 
practically plane polarized in the direction of stress 
through the entire model. 

Mindlin and Goodman 19 derived optical equations . 
of three-dimensional photoelasticity from Maxwell’s electro- 
magnetic equations by neglecting the inplane (plane perpendi- 
cular to the wavenormal) variations of the indeces of 
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refraction and the displacements and some other quantities 
of very low magnitude. They derived Maxwell -Neumann equa- 
tions^ and showed the order of magnitudes of necessary 
approximations involved. They also showed that the Para- 
meter R was the essential one controlling optical wave pro- 
pagation in three -dimensional photoelasticity. 

PI 

Ginsburg derived yet another system of differen- 
tial equations governing the propagation of electromagnetic 

pp 

waves in an anisotropic, inhomogeneous medium. Aben 
simplified these equations by neglecting certain quantities 
and assuming a particular form for the solution. His equ- 
ations are finally reduced to two first order ordinary diffe- 
rential equations with variable coefficients. He proved 
the equivalence of any length of optical path in the model 
to a system containing three parameters called the charac- 
teristic parameters in analogy with the parameters for a 

23 

system of birefringent plates. He suggested that the 

experimental data of these parameters for various wavelengths 

of light would yield additional information regarding the 

* 

stress distribution inside the photoelastic model. 

Frocht and Srinath 10,11,12 £ eve q 0 p e £ a scattered 
light method in which they assumed that; when the amplitudes 
of light vibration along the secondary principal-stress 
axes at the entrance were equal, with the addition of re- 
tardation along these axes at the entrance (i) the amplitudes 
along the secondary principal-stress axes at any point on 
the light path did not vary and (2) whatever retardation 
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was added at the entrance, the same amount appeared as an 
additional retardation along the secondary principal-stress 
axes at any point on the optical path. This is equivalent 
to assuming that the retardation vs. distance (along the 
wavenormal) curve shifts parallel to itself when a certain 
amount of retardation is added at the entrance. These as sum 
ptions however are valid for the cases where the rate of ro- 
tation is not high, and these are "based on the findings of 
Drucker and Mindlin 1 ^ . 

Aderholdt, Mckinney, Ranson and Swinson 1 *^ in their 
scattered light method treated the light components to be 
following the secondary principal -stress axes in cases where 
the rotation to retardation ratio is small. This was also 
based on the findings of Drucker and Mindlin 1 ^ . 

Robert and Guillemet 14 > 15 » 16 followed an altogether 
different approach to overcome the rotational effect. They 

used the elliptic it y measurement technique that was common 

23 

in crystallographic investigations . In this technique, 
the analyser is rotated at a particular speed and the ratio 
of the amplitudes of alternating and direct components of 
light from the analyser is measured with suitable electronic 
instrumentation. The analysis in this method is not direct 
but through the characteristic parameters. This method 
however, does not involve any assumptions as in the above 
methods. 
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Present Investigation 

The objective of the present investigation is 

(i) to develop the necessary equations describing 
optical phenomena in a photoelastic medium in the 

-presence of rotation of the secondary principal-stress 

oq 

axes, in terms of Stokes parameters , 

(ii) to solve these equations analytically or numerically 
for various cases and study the validity of assump- 
tions of previous investigations, 

(iii) to use these equations for proving the existence 
of an optically equivalent system consisting of a 
retarder and a rotator, 

(iv) to devise methods (both with transmitted and 
scattered light) for the determination of the 
parameters of the optically equivalent system 
(characteristic parameters) experimentally, 

(v) to devise new graphical methods to describe 
passage of light through a photoelastic model in 
which the rotation of secondary principal-stress 
axes is present and to apply these graphical methods 
for the determination of the integrated optical 
effect and the characteristic parameters, 

(vi) to develop new techniques of compensation 

(vii) to develop methods for the determination of integral 
fringe orders (upto any value, when there is no 
stress free corner to count from), 
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(viii) to devise and perform experiments to check the 
validity of (a) the governing optical equations 
Cb) the existence of an optically equivalent system 
and the techniques for the determination of the 
parameters of this equivalent system (c) the assump- 
tions of the previous scattered light methods and 
(d) to find the variation of characteristic para- 
meters with the wavelength of light. 

A brief outline of the contents of each chapter is 
presented below. 

In Chapter 2, a discussion on the available systems 
of optical equations in three-dimensional photoelasticity 
is presented. A new set of ordinary, linear simultaneous 
equations are derived kinematic ally in terms of Stokes para- 
meters . These .are very convenient for analysis and application. 
Using these equations, the analytical solution for the case 

where the ratio R of the rate of rotation (0 f ) to the rate 
of retardation p* is a constant (0' can be any arbitrary 
function) is derived. Also presented is a discussion on the 
concept and properties of light ellipse. 

In Chapter 3, the existence of an optically equiva- 
lent model for any length of light path in the photoelastic 
model is proved. The importance of the concept of an opti- 
cally equivalent system is appraised. Some special cases 
in which there exists some relation between the parameters 
of the optically equivalent system (called characteristic 
parameters) are discussed. Two direct and simple methods of 
determining these parameters experimentally are described* 
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In Chapter 4 , a discussion on the previous methods 
of scattered light analysis is presented. The validity of 
the assumptions involved in these methods is investigated 
using the analytical solution derived in Chapter 2, the 
numerical solutions for some cases of loading ( of a rectan- 
gular har under combined tension and torsion) and the concept 
of 'optically equivalent system'. Modifications of the 
methods of determining the characteristic parameters with 
transmitted light (described in Chapter 3 ) for use with 
scattered light are suggested. 

Chapter 5 on Miscellaneous Topics ’deals with the 
following: . 

(A) 'Discussion on various graphical methods in three - 

-dimensional photoelasticity' - In this section a 
brief discussion of various existing graphical 

OA 

methods developed by Jerrard ^ , Robert and Guille- 
met'-^ , Kuske ^5 y Schwieger26 5 Cemosek27 1 is pre- 
sented. A new, direct approach based on the concept 
of light ellipse with two graphical methods is des- 
cribed. Use of these graphical methods in finding 
the characteristic parameters is demonstrated. 

(B) 'Discussion of the existing methods of compensation 
and development of a large field compensator' - In 
this section use of (a) a system of four retarders 
each having the same retardation, with the pola- 
rizer or analyser and (b) a system of eight retar- 
ders each having the same retardation, as a 
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compensator is described. This is essentially an 

QO 

extension of the works of Pancharatman and 
Tuzi 

(C) ’Determination of integral fringe orders in . 

phot oe la stic it y’ - In this section the construction 
of a nomogram and a ready reckoner table which can 
be used to determine the integral fringe order 
(virtually upto any value) from the fractional 
fringe orders for two wavelengths of light, is des- 
cribed. This method is essentially an extension of 

30 

the method developed by Pant . Some experimental 
values (for a circular disc under diametral compre- 
ssion) are presented to demonstrate the value of 
these methods. 

In Chapter 6 are presented results from experiments 
conducted on a rectangular bar under combined tension and 
torsion. These form the supporting evidence for the validity 
of the governing optical equations and the existence of an 
optically equivalent system. The techniques described in 
Chapter 3 are successfully employed to determine the chara- 
cteristic parameters experimentally for various values of 
tension and torsion. These experimental values are compared 
with the values obtained by numerically integrating the 
governing optical equations derived in Chapter 2. The ' 
agreement between these two sets of values is found to be 
good. The characteristic parameters are also found experi- 
mentally for three wavelengths of light for the same loading. 
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The change in the characteristic directions is 3*+° when the 
v, wavelength is- changed from 5004- °A to 5950°A (the correspond- 
ing change' in the material fringe constant is from 128 - 158 
Psi. in. /fringe ) . For the particular values of tensile and 
torsional loads considered, the assumptions of the previous 
scattered light methods are not found to be valid. 




CHAPTER 2 


THE uIGHT ELLIPSE AND OPTICAL EQUATIONS 
IN TERMS OF STOKES PARAMETERS 

Introduction s 

The photoelastic analysis consists of observing 
the changes in the properties of the polarized light in the 
photoelastic model experimentally and relating these to the 
(secondary) principal-stresses and their directions using the 
stress-optic law. This is usually done by resolving trigo- 
nometrically the light vector along different directions and 
examining the amplitudes and phase angles of each component. 

A simpler and more convenient way is to study the properties 
of the light ellipse and its transformations as the light 
passes through a model or a series of optical elements and 
then to use this information in the photoelastic analysis* ^ 

Light Ellipse 

The most general type of polarized light is the 
elliptically polarized light which can be represented by 
an ellipse called the light ellipse. This can be represented 
by two simple harmonic vibrations along two orthogonal refe- 
rence axes OX and OY, Fig. 1 as: 
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E X = ax Coswt 

and (2.1) 

Ey = a Y c os(“ t +e) 

where C is the relative phase difference between the two 
vibrations, and ax, aY are the' amplitudes along these two 
axes. 


In the case of a stressed phot oel istic medium, it 
is conventional to make the reference axes coincide with the 
secondary principal-stress axes. 


Eliminating wt between the above two equations, the 
equation of the light ellipse is obtained as 


pc 


+ 


2 


hi _ 2 % % cose = si„2e 

ay l H V 


( 2 . 2 ) 


The major axis of the light ellipse makes, an angle 
with the X-axis such that 

2 Uy- ay 

tan 2* = - 2 J/Ose 

^ - a Y 


= tan2 e Cos G 


( 2 . 3 ) 


ciy* 

where tane= = amplitude ratio 


( 2 .>+) 


o( is called the azimuth of the light ellipse with respect 
to the reference axis -OX. 

If a and b are the lengths of the semi -major and 
semi -minor axes of the light ellipse, their values are 
given by 
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pop p p 

a = af Cos °( + ay Sin ^ + 2ax ay Cos^ Sino( Cosg 

pop 2 ? (2*!?) 

= a^ Sin <* + ay Cos o( - 2ax ay Cos^c Since CosG 
2 2 2 2 

so that a + b = ay + ay (2*6) 

The ratio ; 

tan a) = ~ = e (2*7) 

cl 


e is called the ellipticity. 

The light ellipse can also be referred to its own 
major and minor axes 5 and n - giving 
= a Cos(“t +C 0 ) 

5 ( 2 . 8 ) 


E^ = b Cos (“t +e o ±rn;/2) 


= ± b 8in(o3t + e o ) 

a.j SinG Sinc< 

Where tan£ 0 = ax Co£W ~ aj S i n « C os e 


(2-9) 


The plus or minus sign indicates left handedness 

or right handedness of the light ellipse. 

a 

The following points are observed and used in the 
subsequent analysis. 

1) When the light ellipse is referred to its own axes, the 
relative phase difference is +v/2* 

2) If E^ is ahead of F^ by ir/2, the light ellipse is right 
handed. If Eg is ahead of E^ by tr/2, the ellipse is left 
handed, Fig .2. 

3) If a phase difference of + tr/2 is added along the major 
or minor axis of the light ellipse, it becomes a line ellipse 
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(linearly polarized), oriented at an angles such that 
tan « = + b/a. 

4) If the amplitudes ay and ay along the reference axes 
are equal, then tan2c< =oo , which means, that the axis of the 
ellipse will be at 45° to the X-axis. 

5) If ay = ay and G = ir/2, the light ellipse becomes a 
circle and a = b (circularly polarized - from Eq. (2.5))* 

The state of polarization as represented by Eqs . 
(2.5) a nd (2.6) can be expressed in terms of the azimuth 
and the eilipticity b/a. These two quantities are sufficient 
to express completely the state of polarization if one dis- 
regards the absolute intensity (as is the case in most of 
photoelastic work). 

Stokes Vector 

The components of light vibration along the refe- 
rence axes OX and OY, Fig. 1, are given by Ey and Ey in 
Eq. (2.1). The components of this vibration along a set of 
new axes CX* and OY 1 , Fig. 3, inclined at an angle ofcc* to 
the 0X,0Y axes are given by 

Ey-' = a Y Cos^t Cos /+ av Sinef 1 (Coswt Cosg - Sin“t Sing) 

A x ^ (2.10) 

Ey! = - ay Coswt Sine* + ay Cost* 1 (Coswt CosG - Sinwfc Sing) 

The amplitudes along the new axes are given by 

ayt = (a x Goa** + ay Sin* 1 Cosg) 2 + (ay Sin*' Sing) 2 

a|-t = (- ay Since' + ay CoSce' Cosg)^ + (ay Cose* 1 Sing) 2 


( 2 . 11 ) 
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. . a^i + a^r = 3x + a Y since the total intensity will not 
change with the rotation of the reference axes. 

From Eqs. (2.5") and using Eq. (2.3), 


a 2 - b 2 = (a^ - &y) Cos2=< + 2a^a Y Sin2o( Cos£ 

= (a| - a.y) Cos2o( + (ajr - a§0 tan2c< Sin2o( 

2 2 
a y - af 

= Cos2°C 

But, from Eqs. (2.4) and (2.7) 

a 2 - b 2 _ Cos2w 
a^a| C#s20 

. . Cos2a)Cos2o( = Cos2e (2.12) 

Multiplying this by Eq. (2.3), we obtain: 

Cos2a)Sin2°( = Sin20 Cos£ (2.13) 

Squaring and adding the above two Eqs., 

Cos 2 2 w = Cos 2 2 0 + Sin 2 2e Cos 2 £ 

or Sin 2 2 m = Sin 2 20 Sirf^G 

Taking the positive square root on both sides, 

Sin2 w = Sin20 Sin£ (2*14) 

The four quantities (a^ ■+■ 3>y)j (^K ** Cos £ 

32 33 34 2 3. 

and 2a^ Sin£ are referred to as Stokes parameters'.^ ’• ' ’’ 
which were introduced by the British physicist G.G. Stokes 
in 1852. The first quantity (a 2 + a|) corresponds to the 
total intensity of light vibration. The variations of this 
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total intensity due to absorption in the medium are not of 
any consequence in most of the photoelastic investigations. 
For convenience, we set this quantity equal to unity. We 
shall represent the Stokes vector by V = (Sq S-] ,$2,83). 


Throughout this discussion only the last three 
parameters will be used ignoring the first one as it will 
not affect the present analysis. The Stokes vector will 

;~sf 

either be written in the form of a column vector 


S- 


or m 


the form of a row vector (S i S Q S_) depending on the convex 

1 > 3 

ience. 


k Now, we have the following equivalent forms for 

the Stokes vector (with a|r + a|. = 1 ), 


"s 1 ” 


2 P 

A - a T 


Cos2 &■ 


Cos2w 

C\J 

CO 

o 

o 

S 2 

- 

2a X a Y 


Sin2© CosG 

5 

3 

C\J 

to 

o 

O 

Sin2c< 

_ S 3 

f 

2ay ay Sin C 


Sin2e Sine 

■*— — » 


Sin2w 



It may now be noted that these parameters are with 
reference to a set of perpendicular axes which make an angle 
of with the major and minor axes of the light ellipse'". If 
these reference axes are rotated, the quantities 0, 8» 

a x and ay change and so the Stokes parameters will also 
change cor re spending ly. 

Also sf + S§ + S§ = 1 (2.16) 

Since S<j = - a^, the first Stokes parameter S^ 
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depends only on the amplitudes of components along the 
reference axes. Also, S^ = Sin2uj and the third Stokes 
parameter depends only on the elliptic it y of the light 
ellipse . 

Transformations of Stokes Vector 

When a certain amount of retardation is a dded 
along one of the reference axes, say, OY-axis, Fig. 3, the 
phase differences , becomes (£ + £,), and the amplitudes 
ay and ay remain unaltered (i.e. © remains unaltered), i.e. 
the parameter S2 = Sin2 ©Cos € becomes Sp = Sin2 GCos( £ + £3) 
and S^ = Sin2 ©Sin £ becomes S^ = Sin2 ©Sin( £+ £3) 



S 1 


s i 


S 1 

i.e. the vector 

s 2 

becomes 

^2 

= 

s 2 Cos Cl - S3 Sin e : 


s 3 


S 3 


S^ Cos £1 + S 2 Sin £;. 

| 

-j 




L (O . 


If the reference axes are rotated in the anti- 
-clockwise direction through an angle Fig. 3, the azimuth 

^ becomes ( °f - =<i) and the ellipticity remains unaltered, 
so that S3 = Cos2 wCos2 a becomes S3 = Cos2“ Cos2(c( - oc^) 
and Sg = Cos2co Sin2=< becomes S2 = Cos2 “ Sin2( o< - <*3). 
i.e. the vector 


V 


S 1 Cos2°fi + S 2 Sin2=<2 


X 

S 2 

becomes 

S 2 Cos2 c ( 1 - S3 Sin2=<2 


^2 

S 3 

■ i 

S 3 

j 

% 


( 2 . 18 ) 
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Now, consider the addition of a retardation of £3. 
along the set of axes OX', OX', Pig. 3, which make an angle 
of =<qwith the original reference axes. For this process, 
the transformation is obtained by first decreasing the 
azimuth °( by an amount °<q without altering the elliptic ity, 
and then increasing the phase difference £.1 by an amount % 
without altering tile amplitudes. This is equivalent to 
the successive application of the above two transformations. 
In this case, the vector C Sq ? S 2., S3 ) becomes, 

~Sl 0os2 *1+ S2 Sin2 <*q "j [s^ 

(S 2 cos2 X,- Sq Sin2 c < 1 )Cos £q- S 3 Sin £q = S 2 (2.19) 

S 3 Cos G 1+ Sin £^S 2 Cos 2°^ q - S 1 SinS^q) S 3 

Here, the vector (S-i ? S 2j S3 ) is with respect to the 
new set of axes OX', OY 1 . If this vector is again referred 
to OX, OY axes, we will have 




(S-) Cos 2 °Cq + S 2 Sin2o(q)Cos2c3C 1 - Sin 2 o(q^S 2 Cos 2 oc -j 



- S 1 Sin 2 c<q)Cos Gq - S 3 Sin G]_ 

c\l 

m 

= 

[(S 2 Cos2c(^ - S^ S in2 oc-^C 0 s - S^ Sin g-JjoosScc^ 



+ Sin 2 c< 1 (S'| CosGoc-^ + S2 Sin 2 o(-j) 



S3 CosGq + Sin£ 1 (S2 Cos2a - S^ Sin 2 oCq) 


( 2 . 20 ) 


Intensity Calibration Method To Determine Ellipticity 


Consider a system consisting of a xA plate and 
an analyser, the angle between their axes being k- 5 °* Let 
the Stokes vector for the light ellipse before this system 
with respect to the axes of the X-A plate be A-^ ,= (Si ? S 2j S3). 
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Let be the Stokes vector for the light ellipse after 
the plate referred to the (V4) plate axes and , the 
Stokes vector of the light ellipse after the (*/+) plate 
referred to the axes of the analyser. Then the transfor- 
mation through this system is given by 


M 


B- 


S-| \ 

' 

n 

°1 


” ~ S 3~ 


X 

s 2 


"S3 


-Si 


S 2 

®3. 


*2 


s 2 


X 


( 2 . 21 ) 


Using the transformation laws given by Eqs. (2.17) 
and (2.18). If 




— 2 2 
and S^ = a x - a^ 


sf 


C os2 u Cos2 

S 2 

= 

Cos2w Sin2 °c 

S 3. 


Sin2w 


then, S-] = - Sin2e = 2a£ - 1 
or ag - 


since 


a| + ay = 1 


Also, the intensity along the OX -axis is given by (QX,OY 
now refer to the analyser axes) 

I = K a| = | (1 - Sln2w) (2.22) 

where K is the constant of proportionality. 

So, we have the following theorem; 

Theorem 1 ; Consider a system consisting of a ( V+) Plate 

and an analyser, the angle between their axes being 45° • 
The intensity of light I from the analyser is dependent 
only cn the elliptic it y of the light incident cn the 




s “ Source of light 
P ” Polarizer 

tt/4),,(V4) 2 - First and second quarter-wave P | ales . 
M — Model 
P ~ Fast axis 
S - Slow axis 

P,q~ Principal- stress directions 
A — Analyser 



OPTICAL SET-UP FOR TWO 
STRESS ANALYSIS. 


~ DIMENSIONAL 
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( X A ) plate and is given by Eq. (2.22). This theorem can 
be used with advantage in two-dimensional stress analysis. 
Consider the optical set up shown in Fig. U-. The light 
after the (X/+)i plate is circularly polarized so that the 
amplitudes along the p,q-axes are equal. Along p,q-axes 
the phase difference is ir/2 before the model and (ir/2 + p) 
after the model, if the amount of retardation added in the 
model is p . From the property of light ellipse, the ellip- 
ticity of light after the model is tan (irA + p/2). Using 
Theorem 1 , the intensity I from the analyser is given by 
Eq. (2.22), where tan to is the ellipticity. 

8 o 

1=^(1- Sin (tt/ 2 + p ) = If (1 - Cos p) (2.23) 

This equation relates the retardation introduced 
in the model which is proportional to the principal-stress 
difference (p - q) to the intensity I from the analyser. 

To find the constant of proportionality K, Tardy's method 
can be used. For this, remove the (* A ) 2 Plal e and rotate 
the analyser for minimum intensity, thus obtaining the 
directions of the major and minor axes of the light ellipse 
after the model. Now set the axes of the ( V *+) 2 Plate 
parallel to these major and minor axes so that the light 
after the ( V +) 2 Plate is plane polarized. Rotate the 
analyser for minimum intensity, thus finding the direction 
of plane polarization. The angle between this direction 
and the major axis will be w where tan “ is ellipticity of the 
light ellipse after the model and 0 = ( p + ir/2)/2. Now 
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arrange the ( ^/k plate and analyser as shown in Fig. if, and 
note the intensity I. Since the quantities I and p are known* 
the constant K can be evaluated using Eq. (2.23). Once the 
value of K is kno\>n, the principal-stress difference can be 
computed using Eq. (2.23) from the value of + for any other 
point on the model. A calibration graph showing the relation 
between I and (p - q) can be drawn for convenience. 

Optical Equations of Three-Dimensional Photoelasticity 

Several systems of equations have been put forth 
to describe the optical phenomena in the presence of the 
rotation of secondary principal-stress axes. We shall give 
a brief account of some of these derivations. The first of 
these are the Maxwell -Neumann equations; 
da = - a Cos djZ5; da = - a Cos 60 (2,2k) 

and d = (~ - -—) Sin 60 + dp (2.2^) 

a x 

where 60 is the change in the orientation of the secondary 

principal-stresses and dp is the rdtardation introduced 

(due to stresses only) in an infinitesimal length dz of the 

20 

optical path in the model. Coker and Filon derived these 
equations kinematically. They assumed the components of the 
light vector in the trigonometrical form. Also, in this de- 
rivation, after a distance dz along the optical path in the 
model, the reference axes were first rotated to coincide with 
the directions of the new principal-stress axes and then the 
retardation caused due to stresses alone was added to obtain 
the components of the new light vector. 
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Drucker and Mindlin"f ormulat ed the problem in terms 
of the elastic ether theory of light. The ether theory is 
based upon the ideas of stress , strain and displacement 
and upon the equations of motion. The governing differen- 
tial equations are exactly the same as those obtained from 
the electromagnetic theory of light if the symbols repre- 
senting displacement, density and shear modulus are inter- 
preted as magnetic force, magnetic permiability and reci- 
procal of dielectric coefficient respectively. They obtai- 
ned the analytical solution for the case, where the secondary 
principal-stress difference remained constant and their 
orientation varied linearly along the light path. They 
concluded the following ; (1) Rotation increases the i rela- 

r p ^ 2 t T 

tive phase retardation by a factor S = 1 ( ^ ) j 

which at times may be quite large, (pis the orientation 
of the secondary principal-stresses for the case under 
consideration. & is the relative retardation caused by 
the magnitude of stress alone). If the ratio -y-^is not 


constant, the ratio of the increment in 2p to the increment 
in A between fringes can be used in regions of high stress. 
(2) At high stress, the light vector rotates with the 
principal-stress directions, plane polarized light entering 
along a direction of principal-stress will remain practi- 
cally plane polarized in the direction of the stress through 
the entire model. 
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Mindlin and Goodman derived a set of dynamical 
equations of three-dimensional photoelasticity using the ; 
Electromagnetic theory, with the purpose of revealing -the 
approximations involved in the kinematieal derivation. 

For this they made an analysis of the orders of magnitudes 
of the quantities involved in the equations. The orders 
of magnitude of these qLiantities are presented below. 


(i) 

b 2 

"i 2 

a s 

= o 

X 2 

os 

(ii) 

rO|l 3 

i_i 
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(iv) £ bz 
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0 




( 2 - 26 ) 


(v) 




(vl) z 6c 2 6b" 1 

To os os 

S, n and b are representative lenghts in the 
X,Y plane (the plane of the wave front), index of . refrac- 
tion and wave velocity respectively, X is the wavelength 
of light in vacuum and 0 £ ^stands for n order of magnitude 
of”. To estimate the order of magnitude of these ratios, 
consider a wavelength of 6x1 0 -1 *' mm, a stress gradient of 


xz 


OS 
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Q 2 

10° dynes/cm2(~ I^OO lb /in 7) per mm. and a stress-optical 


coefficient of 100 Brewsters. Under these extreme conditions 
' ^ and, for a model 5 nan. thick, z = 0 1 10 | 

With regards to terms involving j# , if the rotation of the 




0 


10 


principal axes were 2 ir radians in one millimeter, 


0 


$S. 


10 


-2 


Since the ratios involve squares or 


products of terms not exceeding 0 


selves donot exceed 0 


~ 10 " 6 ] . 


10 


-3 


the ratios th em- 


it is true that in the 


neighbourhood of an isotropic piont , becomes very large. 

However, may be expressed in terms of ~ and the principal 
Q S 0’S 

stress difference through the use of the stress equations 
of equilibrium and the stress-optic law, following which 
it may be shown, by a limiting process, that the ratios 


involving off are, again, small. Also, the order of magni- 

3S 

tude of (n-j - n 2 )/n and * d& (hi ^ are the r efractive 

indeces) is|10 _ ^|. It is not impossible fordj# to become 

L J dz 

very large. An example is the case of combined tension 

and torsion of a cylindrical bar. For a wavenormal inter- 
secting the axis of the bar^! approaches infinity as 

dz 

the tension approaches zero and the product of the two has 
a zero limit in the neighbourhood of the axis of the bar. 
They simplified the governing equations by neglecting these 
quantities. They obtained the Maxwell-N'eumann equations by 

assuming the solutions in the trigonometric form. By 

. ' ■ . 

changing the independent variables, they obtained the 
following simplified form of the optical equations.. 
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R u'= 2-J 

9 t 

R ■y' =: — 

9 a 


where 


R 




(2k /\) (n 2 - n«| ) 
cr -H (t - /KLj'"' dz) 


to(t - JK 2 dz) 



-1 -2 ? -1 -2 
v n-j , KZ, ~ y n 2 


(2.27) 

(2.28) 

(2.29) 

( 2 . 30 ) 

(2.31) 

( 2 . 32 ) 


In the electrostatic system of units y = v q 
where v q is the velocity of light in vacuum. u r and v* 
are the components of the light vector along the secondary 
principal-stress directions. 

Eqs. (2.27) and (2.28) show that R is the essential 


parameter controlling optical wave propagation in three- 
-dimensional photoelasticity. Thus, we see that the kine- 
matical equations for all practical purposes are accurate 
enough as can he observed from the assumptions involved in 
the dynamical derivations. 


Aber?2 derived yet another system of differential 
equations governing the optical phenomena in three- 
dimensional photoelasticity. He used, the system of equa- 
tions for the propagation of electromagnetic waves in an 
anisotropic inhomogeneous medium obtained by G-insburg 21, 
Using the fact that the photoelastic medium is cnly 
weakly anisotropic, ' he obtained the following system of 
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ordinary linear differential equations. 

dBi Ar* 

aT = "I C ('i - a 2> si +fj bJ (2.33) 

and 

i®£ = - J* bI 4c' ( cr a ) B ' ? (2.3*0 

dz dz 1 2 1 2 a 

where C* = CC Q (2.35) 

C is the velocity of light in vacuum and C 0 is the photo- 
elastic constant, aq are the secondary principal- 

"Stresses and bIj , are the components of light vector 
along the directions of the secondary principal-stresses. 

Equations Using Stokes Parameters: 

We shall now obtain a kinematical derivation of 
the governing optical equations in the presence of the 
rotation of the secondary principal-stress axes to faci- 
litate a better understanding of the optical phenomena. 

An observation of the transformation Eq. (2.19), 
shows that it is linear (i.e. no products, squares or 
trigonometric functions of the Stokes parameters are invol 
ved). This prompts us to use the Stokes parameters direct- 
ly as the variables in the kinematic derivation. 

Let the Stokes vector referred to the secondary 
principal-stress axes at points 0 and o'be (S^ , S 2 ? S 3 ) 
and (S 1 , S 2 , S 3 ) respectively. The point 0 is at a 
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distance z from the origin -which is at the point of entrance 
to the model, Fig. 5* Let the change in the orientation of 
the secondary principal-stress axes over the distance dz he 
d 0 and the retardation introduced due to stresses alone he 
dp . Then using the transformation represented by Eqs. 
( 2 . 19 ), and noting that for sufficiently small distance dz, 
the quantities and dp will be very small so that 

Sin 2 $ 2 dp, Sin dp ^ dp, Cos2d^ ^ 1,Cos dp ^ 1, dp dp * 0, 
we obt ain 


~Si‘ 


s 1 + dS-f 


S 1 + 2S 2 dp 

s 2 

~ 

S 2 ■+■ dS 2 

= 

S 2 - 2S 1 dp - S 3 dp 

S3 

H 


S 3 + dS 3 


S3 + S 2 dp 


( 2 . 36 ) 


Dividing each of the above equations by dz and taking the 
limit as dz •+• 0 , we obtain, 

! dS 1 H0 

S 1 = = 2S 2 tz = 2S 2^ (2-37) 

s' = = - 2S- M - s. = - 2Sj' - S. P' (2.38) 

2 dz 1 dz 3 dz 1 8 

and 

S 3 = = S 2 = S 2 p* (2.39) 

D dz ^ dz ^ 

If the Stokes vector (S-j , S 2 , S 3 ) is referred to 
the fixed axes CK, OY, we have from Eqs. (2.20) (replace 
=< ^ by p + d |5 and £ -] by p) 

5 1 + dS a “ 

5 2 + dS 2 

5 3 + dS 3 


S 1 

®2 
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(S 1 Cos2Q# + 60) + S 2 Sin2(d^ +0)) Cos2(J# + 60) 

- Sin2(^ + 60) $(S 2 Cos2<$ + 60) - S 1 Sin2(j2l + d0)) Cos dp 

- S 3 Sin dpj 

Cos2( 0+60) ^ (S 2 Cos2( 0+60) - S 1 S±n2( 0+ 60) Cos dP 
= - Sin dpj + Sin2 ( 0+60) (S 1 Cos2( 0+60) \ 

+ S 2 Sin2( 0+60) 

S 3 Cos dp + Sin dp (Sg Cos2( 0+60) - S Sin2( 0+ d ft 
~ (2.40) 

As above, for sufficiently small distance dz, 6 0 
and dff are very small so that Sin dp ^ d , Cos: 260^1, 

Cos dp^ 1, 60 dp ^ 0, Cos2( 0+600 Cos2 0 -y2d 0 Sin2 0 ; 
Sin2(dj2 +0) ^260Cos20 + Sin2 ^ so that 


Si 


S 1 + dS 1 


S-] + S 3 dp Sin2j# 

s 2 

= 



S 2 - S 3 dp Cos2j# 

V 


s 3 + ds 3 


S 3 + S^ dP Cos2$ - S^ dp Sin2j# 



— _ 




( 2 .4 1 ) 


Dividing each of the above equations by dz and talcing the 
limit as dz — *» 0, we obtain 


IS-, 
S 1 ~ dz 
, aso 


S 2 

and 

s' 

S 3 “ dz 


S 3 P* Sin 20 

- S^>’ Cos2jZi 

- S 2 p* Cos2p - S-j P* Sin2^S 


(2.42) 

( 2 . 43 ) 

( 2 . 44 ) 


These equations can also be obtained from Maxwell - 
-Neumann equations by converting the dependent variables 
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into Stokes parameters. 

Eqs. (2„*f2) - (2**40 will be particularly suitable 
for scattered light analysis. This will be discussed in 
a later chapter. The above sets of equations are. linear, 
ordinary simultaneous first order differential equations 
which are very convenient for analysis and application. 
Using the Eqs. (2.37) - (2*39), we shall obtain the analyti 
cal solution for the case where the rate of rotation fl can 
be any arbitrary function, but the ratio R =^- is a cons- 
tant, (The case studied by Drucker and Mindlin where both 
and p' are separately constants is a particular case of 
this) . 


Let K 2 = 

4r 2 


from Eqs. (2.37) and (2-39) 

s', 

= — 

3 2R. 


( 2 .* 0 ) 


Integrating both sides, 

S. 


S. 


2R 


+ A, 


where A-j is an arbitrary constant. 

Substituting this expression for S^ into the equation 

O O O 

a, + S2 + s 3 — 1 

(instead of using Eq. (2.38)) 
we obtain an expression for S in terms of S ^ . 
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Substituting this expression in Eq. (2.37) and integrating, 
we obtain 


A-i 1 A-i ? 

S 1 =“ i p + K (1 ""3 ) Sln(2Kp+A 2 ) (2.4-6) 

2RK C ~ K^ 


where A 2 is an arbitrary constant. 

Using Eq. (2.37) again, we obtain 

S„ =% = (1 Cos (2KJ3+ Ap) (2.4-7) 

2 2 ^ 2 

and (1 - h )i 

S 3 = 2R + ^ ~ ^2 + 2RK Sln *“ 2K (2.48) 

The arbitrary constants A^ and A 2 are to be evalua- 
ted from the initial conditions. 


Jones Vector 


It is sometimes more convenient to use Jones 


op 

vector 0 in photoelastic analysis. 

This vector has two 

e lement s 


i(e Y + 7t) 

P = a x e 

(2.49) 

and _ 

i(£ +«t) 

Q = a e x 

( 2 . 50 ) 


The quantities ay are the amplitudes along two 

reference axes OX, OY; and are the P^ a se angles of 

the instantaneous components along these axes. 

i(e^ + wt) iifc i% 

_ Since e = e e and the term, 

in 

e occurs in the expressions for both P and Q this can 



38 


be ignored in the photoelastic analysis where the varia- 
tion with time is of no consequence. So, we have 

P = a x e 1WX (2.51 ) 

i£ 

Q = a y e (2.52) 

It can be noted that the transformations corres- 
ponding to the rotation of the reference axes will be 
linear in P and Q. So, if we derive the optical equations 
of three-dimensional photoelasticity in terms of the ele- 
ments of the Jones vector, -re should obtain ordinary linear 
differential equations in P and Q. 

Now, referring to Fig. 5 ? let (P,Q) be the Jones 
vector -with respect to the axes ON, OV at the point 0.' If 
this vector is referred to the axes OU 1 , OV* its elements 
b ec ome 

P' = P Cos A 0+ q SinAjZf 

Q r = - P Sin A 0+ Q Cos A fi 

Now, if a retardation of Ap is added to the P 1 

component, the elements of the new vector are given by 

iAp iAp iAp 

P l * = P ! e = P e CosA fS + Q e Sin A 0 

and Q m = Q f 

For sufficiently small distance Az, Ap and A 

are very small so that Sin Ap ^ Ap , Sin A fi ^ ft, Cos Ap ^ 1 , 

iAp 

CosA^vl, Ap A4>^0, e v 1 + iAp . We, then have 

P M = P + AP = P(i + iAp ) + Q 

and q" = Q + a Q = - P A $ + Q 
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Dividing these equations by az and taking the limit 


tS3 

CO 

c3 

0, ■ we 

■ obtain 

SS = P=L 

dp p 

+ Q — ■ 

dz 

dz 

^ dz 

and 



IS _ Q»_ 
dz “ ^ ~ 

~ d0 
dz 

1 

li 


(2.53) 


(2.5V) 


Instead of adding a retardation of Ap to one com- 
ponent only, if we add retardations £ 2 - and -Ap/2 to the two 
, „ . iAp/2 2 -i Ap/2 

.ompcnents (i.e. P = P e , Q = Q e ), we obtain 

P'= O'P + 0' Q (2.5?) 

and Q* = - 0 ' P + | p* Q (2»56) 

These equations are same as those derived by 

Aben . 


In this chapter, we have introduced the concept 
•f light ellipse as an elegant tool in photoelastic analy- 
sis. We have given a brief review of some of the previous 
derivations of the optical equations of three-dimensional 
photoelasticity to show that, for all practical purposes, 
the kinematically derived equations are sufficiently accurate. 
We have also derived two new sets of equations for three- 
-dimensional photoelasticity, which are convenient for 
analysis and application. Using these equations., we obtained 
the analytical solution for the case where the rate of rota- 
tion 0* can be any arbitrary function, but the ratio 
R = 0 1 / p * is a constant. The contents of this chapter form 
the basis for subsequent chapters. 



CHAPTER 3 


CHARACTERISTIC PABAMETERS IN THREE-DIMENSIONAL PHOTO- 
ELASTICITY AND THEIR EXPERIMENTAL DETERMINATION 

Poincare 1 Sphere Representation ; 

In this chapter, we shall discuss the existence of 
an optically equivalent system, its theoretical and experi- 
mental determination and its importance in three-dimensional 
photoelastic analysis. In crystallography, it is establi- 
shed that a series of retarders is equivalent to a system 
consisting of a single retarder (an optical element which 
adds a certain amount of retardation to the components of 
the incident light along a particular set of axes) and a 
pure rotator (an optical element which rotates the incident 
light ellipse through a certain angle known as the rotatory 
power). This has been proved geometrically by using the 
Poincare* sphere representation of polarized forms. We shall 
present this proof below. 

An elliptically polarized light can be characterised 
by two quantities (i) the azimuth and (ii) the ellipti- 
city tan o. This can be represented by a point P on the 

t 

surface of a sphere of unit radius called the Poincare 
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sphere with co-ordinates : x = Cos2w Cos2°<, y = Cos2t») Sin2°c, 
and z = Sin2w(these three correspond to the three Stokes 
parameters , S 2 and S^ of the light ellipse), as shown 
in Fig. 6. (Let us recall that the first Stokes parameter 
S^ = Cos2w 0032^ corresponds to the light intensity along 
one of the reference • axes. If the amplitude of light •vib- 
ration along X-axis is and the intensity is I, then 
I = K a| = — - - ---- ) . The angle 2c< is measured clockwise 

from H. The angle 2& is measured upward from the equator,. 
The upper and lower poles and P 2 represent left and 
right circularly polarized light (2w = 90 °? tan co = 1 ) . 

Points on the equator indicate polarized light (plane, 2o = 0 
tan o)-0). The action of a linear retarder can be con- 
veniently represented on Poincare' sphere. Consider an 
elliptic ally polarized light represented by the point P, 

Fig. 8. with ellipticity tan to and azimuth °c, to be incident 
on a linear retarder with retardance A and azimuth o 0 j_ . 

Its axes are represented by line RPLj in the figure. The 
x co-ordinate of the, point P corresponds to the intensity 
of this light vibration along one of the reference axes. 

The projected length of the radius vector OP along the 
axis ROR^ corresponds to the light intensity along one of 
the retarder axes. (This can be proved as follows; X - Y 
plane is shown in Fig. 7(a). OX 1 , OY 1 are the new reference 
axes with XOX' =23: . P is the projection of P on the 

X - Y plane. The projection of "P on OX* is given by, 
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x' = x Cos2°c 4- y Sin2 5 Cos2 5 + S 2 Sin2 Sc = S-j (say). But, 
S-j is also the first Stokes parameter obtained after rota- 
tion of the reference axes through an angle °c . i.e. the 
projection of the point P (which again is the projection of 
P co X - Y plane) on any axis OX* in the X - Y plane corres- 
ponds to the light intensity (of light vibration represented 
by P) along the same axis OX 1 ) . 

Since this intensity remains constant during the 
addition of retardation, the process of addition of retar- 
dation in the retarder is represented by an arc from the 
point P to P* about the axis ROR] . This arc PP* on the sphere 
makes an angle of & an the axis ROR-j . Rotation of the 
reference axes through an angle d°< corresponds to rotation 
about the Z-axis through an angle 2d°( cn the Poincare* 
sphere. During this rotation, the z- coordinate andJ so the 
ellipticity remains unaltered. 

The process of addition of retardation along the 
reference axes can be explained as follows. For this, 
consider Fig. 7(b). Point P represents a certain ellipti- 
cally polarized light referred to axes CK, OY in the figure. 
For point P, x = OM, y = 0, z = MP. Now rotate MP about 
axis OX through an angle A so that its new position is HP' • 
Now, for point P*, the co-ordinates are x* = OM* = OM, y = 

- MM* = - MP* Sin a = - z Sin A , z - MP* Cos A = z Cos A . Now 
consider the three Stokes parameters S-j of the pola- 

rized light represented by the point P. These are S^ = OM 
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= x, S2 = 0 , S3 = MP = z. Consider the polarized light P' 
obtained by the addition of a retardation of a along the 
reference axes. The Stokes parameters for this are given 
by s'-| = x, S2 = - z Sin a and S3 = z Cos A, The same values 
are obtained by rotation about the OX -axis. Thus, the effect 
of a single retarder of retardation A can be represented on 
the Poincare' sphere by a rotation through an angle A 
(upwards) about the corresponding axis in the equatorial 
plane. 

Combined Effect of a Series of Retarder s: 

The combined effect of a series of retarders can be 
represented by successive rotations of the Poincare' sphere 
about the corresponding axes in the equitorial plane. This 
in turn can be replaced by a single rotation about some 
general axis in Poincar 4 sphere. Alternatively, this resu- 
ltant single rotation of the sphere can be resolved into 
two rotations about perpendicular axes - the first may be 
about the Z-axis and the other will then be about some axis 

in the equitorial plane which may be determined by the 
23 

construction of Fig. 8. A rotation about AO through twice 
the internal angle at A followed by a rotation about BO 
through twice the internal angle at B is equivalent to a 
rotation about CO through twice the external angle at C. 

Thus the combination of a series of linear retar- 
ders is equivalent to a system containing two elements, one 
a rotator and the other a retardation plate* Since finite 
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rotations about nonparallel axes are noncommutative opera- 
tions, it is necessary to specify the exact sequence of the 
various elements. Interchanging any two of them would in 
general lead to a change in the final state of polarization. 

The above result has also been proved by the matrix 

qcr 

method by Hurwitz and Jones . 

Effect of a Series of Retarders Using Stokes Vector : 

The above result will now be proved in detail using 
the Stokes vector ignoring the variation in the total inten- 
sity. The transformation tint the Stokes vector undergoes 
when a given elliptic ally polarized light traverses through 
a two-dimensional photoelastic model can be obtained from 
Eq. (2.20) as: 

V = M 1 V 

where 



V 


s 1 

V = 

S 2 

; v = 

s 2 


S 3 

L 



are the initial and final Stokes vectors and Rj is given by 

(Cos 2 2°^ + &Ln^2°tj CosCj) Sin2°^ Cos2°tj Sin^ Sin2° ( -j 

(1 - CosC 1 ) 

Sin 2°^ Cos2°tj (1 - CosG-j ) Cos 2 2°C 1 Cos£j - Sin£, 0032 °^ 

+ S±n 2 2^ 

- S±n2 e <j. SinG 1 Cos2^ SinG n Cose 1 
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In the above matrix ®Cj and are the azimuth and retardance 
of the retarder respectively. 

This means, that the effect of a two-dimen si onaJL 
model at a particular point can be represented by a 3 x 3 
matrix. If the light represented by the vector V is inci- 
dent on n successive models, the out coming light vector is 
given by 


V = (M n ....M 3 M 2 M 1 ) V = MY (3.1) 

where M = M n ....M 3 M 2 M 1 

This means that the effect of a series of two- 
-dimensional photoelastic models can also be represented 
by a 3 x 3 matrix. 

Now consider the relation (3-1) 


m-^ -j m-| 2 m*j 3 


m 11 m 21 m 31 


M - m 21 m 22 m ^ ; M 1 = a , 2 

m 31 m 32 m 33 “13 m 23 m 33 


Now Y - (^ 1 ...,M 3 M 2 M 1 )” 1 (?) 

- -1 -1 _ 

= (K, M 2 ....) (V) 


-1 T 


One can show from the M-j matrix that Mj = 
Similarly = mJ; ; M 3 "^ = M 3 etc. Hence 


V = (M^ M 2 ....) (V) 


(3.2) 


= (^ U n _ 1 .*.,K 2 K { ) T (V) 


i.e. V = M T ? 


(3.3) 
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(It may be noted that the property of a retarder 
is unchanged even if the direction of light is reversed. 
That is, the amount of retardation and the direction along 
which this retardation is added remain unaltered even if 
the direction of light is reversed. If the light path is 
reversed for the same series of retarders we then have 
V = (M^ M^M^. . . .M^) where is the Stokes vector of 
light incident on the retarder and V-j is the Stokes 
vector of the outcoming light from the retarder ) . 


So S. = m. . S. 

i i j j 

i, j = 1,2,3 

(3A)- (3.6) 

and S ± - m^i S-j 

i } 2 = 1,2,3 

(3.7)- (3.9) 

Also, we have the relations. 

f 


„2 2 
s 1 + s 2 + S 3 — 1 


( 3 . 10 ) 

2 2 2 
s 1 + s 2 + s 3 = 1 


(3.11 ) 

using the above Eqs. (3. 4) ■ 

- ( 3 . 11 ), we obtain 

the follow- 

ing equations in m , ±,j = 

If S 1 = 1 , s 2 = S 3 = 0 

1,2,3. 


_2 _2 2 2 2 2 

S.j + S 2 + S 3 = m-j^ + m 2 ^ + m 3 -] = 1 

( 3 . 12 ) 

2 

= 0 = S 3 | S 2 —Ij fti"| 2 

2 2 
m 22 + ^32 = 1 

(3.13) 

S^j — 0 — ^2 * ^3 z ~ ^ 5 ^13 

2 2 _ -| 
m 23 + m 33 " 1 

(3.1^) 

If ^ = 1 , s 2 = s 3 = 0 



S 1 + S§ + S3 = n£, + m ^ 2 + >»13 = 1 

(3.15) 



FIG- 9- REPRESENTATION OF THE AXES OF 
THE BIREFRINGENT PLATES IN THE 
TWO EQUIVALENT SYSTEMS. 
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= S3 = 0 , S 2 = 1 ; m|i + 11122 + m§3 = 1 ( 3 . 16 ) 

_ __ _ 2 2 2 

S-[ = S2 = 0 , S3 = 1 ; m3-] + 11132 + 11133 ~ ^ ( 3 . 17 ) 


Thus, we have six relations connecting the nine 
quantities i,j = 1,2,3 so that the matrix M has only 

three independent quantities. So, a series of two-dimen- 
sional photoelastic models (for a particular light path) or 
birefringent plates, in general is equivalent to an optica- 
lly equivalent system with three independent parameters. A 
single birefringent plate has two parameters namely, its 
azimuth and its retardation p . A pure rotator which 
rotates the incident light ellipse has one parameter 4> 1 
namely its rotatory power. The most simple optical systems 
consisting of only birefringent plates and pure rotators 
and having three parameters are (i) a system with a bire- 
fringent plate and a pure rotator in that order and (ii) a 
system with a pure rotator and a birefringent plate in that 
order. 

If the systems (i) and (ii) are to be equivalent, 
the pure rotatory power in both systems must be same, the 
angle between the axes of the birefringent plates in the 
systems (i) and (ii), Fig. 9* must be equal to this rotatory 
power and the retardations of the plates in both systems 


must be equal. This can be proved as follows. For the 
first system, let the rotatory power be , the azimuth of 


the birefringent plate be p and its birefringence be 5. 

For the second s y t" * th'eb'e"'' nuah titles be , 6* amd 

n. i, T. KANPUR I 1 ■ ' 

I CENTRAL. U BftARV 1 
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respectively. Now consider polarized light to be incident 
along one of the axes of the birefringent plate in the 
first system, so that after this plate, it still remains 
plane polarized along the same direction. The pure rotator 
rotates it through an angle of ^ so that at the exit of 
the first system we have plane polarized light at an azimuth 
of ( f ) . In the second system also, let the incident 
light be polarized at an azimuth of P . The pure rotator 
makes its azimuth ( P 4 - ^). This light should be polarized 
along one of the axes of the birefringent plate in the 
second system so that we have plane polarized light at the 
exit of the second system also with the same azimuth as in 
the first system. Therefore, P + ^ = P + ik or iK = ip* . 

That is, the rotatory powers in both the systems are same 
and the difference between the azimuths of the birefringent 
plates in the two systems is equal to the rotatory power of 
either system. Now consider circularly polarized light at 
the entrance of both the systems. For the ellipticities of 
the light ellipses at the exit of these systems to be equal, 
the birefringence of the plates in both the systems must be 
equal. That is, 5 = 5* 

Optically Equivalent System for a Cintinuous Body : 

Now, we shall discuss the existence of an optically 
equivalent system for continuous systems. Aben ^proved the 
existence of an optically equivalent system (consisting of 
a retarder and a rotator) for an optical path in a three- 



-dimensional photoelastic model along which the secondary 
principal-stresses and their directions vary continuously. 
Since the governing differential Eqs. (2.34) , (2.35) are 
linear, ordinary simultaneous differential equations? the 
light propagation is governed again by equations of the 
form (3*4) -(3 . 6 ). The rest of the proof is similar to the 
one presented by Hurwitz and Jones 35 . Aben called the two 
axes of the retarder as primary characteristic directions 
and its retardation as the characteristic retardation. He 
also called the mutually perpendicular directions which 
made an angle equal to the pure rotatory power with these 
axes as the secondary characteristic directions. These 
three quantities namely, the primary characteristic direc- 
tions, the secondary characteristic directions and the 
characteristic retardation, are referred to as the charac- 
teristic parameters. His proof for the existence of these 
parameters is based on the theory of unitary matrices. It 
may now be noted that, plane polarized light entering a 
photoelastic model along one of the primary characteristic 
directions at the entrance comes out as plane polarized 
light at the exit along one of the secondary characteristic 
directions. The corresponding primary and secondary charac- 
teristic directions are termed the conjugate characteristic 
directions. The existence of these characteristic para- 
meters in the case of continuous systems will now be proved 
in detail using the optical equations in terms of Stokes 
parameters. Consider Eqs. (2.37) and ( 2 . 38 ), we have: 
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r i 




20 ' 


= S 2 = - 2S 1 0' - S 3 p' 


The primes denote differentiation with respect to 
the distance co-ordinate 2 . 


>1 


_ 20 ! . 


+ 2S-J0' 


i 

^ P 


- S 2 p - - 20' 


using Eqs. (2*39) and (2.37)* 

On rearranging various terms and simplifying, we obtain 


S 1 - S" 


20" _p" 

r + e’ 


+ s. 


^ - p' 2 + ^ n 


0 


0 ’ 


0 ' 


+ 
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1 


0 U 0' - -0-- 2pt> 
pV. 


= 0 


(3.18) 


This is a third order ordinary linear differential 
equation in with variable coefficients. The initial 
conditions are; 

For a given initial set of (S-j, 82 ) 83 ) i.e. for a known type 
of incident light, the following quantities can be deter- 
mined at z = 0 since the values 0 and p are also assumed 
to be known; 

= 2S 2 0' and s” = 2S 2 0 n + 20,’ S 2 = 2S 2 0 n 
- 20 * ( 2 S-^ 0 ' - 83 ?’) are also known. 

The solution for this initial value problem exists 
and is unique over the closed interval in which the co- 
efficients in Eq. (3.18) are continuous. The solution will 
be of the form 


(3.19) 


S 1 = A-j f| (z) + A 2 f 2 (z) 4- A 3 13 ( 2 ) 


( 3 - 20 ) 
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where A^ } and A3 are arbitrary constants which are to be 
determined from the initial conditions given by Eqs.(3«19). 


i . e . S-) 


• z=0 


= A-j f-j (z) 


z=0 


+ A 2 f 2 ( 2 ) 


' z=0 


+ A 3 f 3 (z) i z=0 ~ S 1 


Si = A-j f i ( z ) 
z=0 


z=0 


+ A2f2(z) + A3f3(z) = 2S20 f 

2=0 2=0 z =9 


— tt it 

S-| = A'j f >[ ( 2 ) 

2=0 . 


2=0 


+ A‘2 f 2 ( 2 ) 


2=0 


+ A3^3(2) 


2=0 




2320 Lo + 


9-0 ' ( - 2S-i 0 -SoP’ )] 

•P J 7,-1 


z=0 


From these relations, the constants A^ ? A 2 and 
can be obtained as linear functions of the initial Stokes 
parameters* Aslo at exit, we have 



(2SiJZ> + S 2 ) 

and S-3 = - from Eqs. (2.81 ) and (2.82) 

p f 

where 0 1 , p* , are the values at exit. 

Hence, "s^ , "§2 and S3 are linear functions of the initial 
Stokes parameters so that we can write 

( 3.21 ) 

where M is a 3x3 matrix. 

The elements of the matrix M are not functions 
of the initial Stokes parameters S^ , S2 and S3. They are 


s n 

5 2 
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M 


>3 
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however functions of 0 ' and p’ and also of the wavelength 
of the light that is used. Since Eq. (3.21) is of the same 
form as Eqs. ( 3 ) - (3.6), all the conclusions regarding 
the existence of an optically equivalent system and the 
characteristic parameters follow as derived in the previous 
section. We also note that these parameters are dependent 
on the wavelength of light that is used since the elements 
of the matrix M depend upon the wavelength of light. 

Some Special Cases: 

For some particular cases of stress distribution 
in the photoelastic model, there exist some relations bet- 
ween the characteristic parameters. In such cases, the 
experimental determination of these parameters is direct 
and simple. 

Case (a)„ Pure Rotator: 

We shall now show that a system of two half-wave 
plates with the angle between their axes being consti- 
tutes a pure rotator. If V is the Stokes vector at the 
entrance of this system referred to the axes of the first 
half-wave plate, the transformations through this system 
are given by 


C os2 w Cos2 °< 


C os2oi Cos2o( 


~Cos2to Cos( -2°0 

Cos2<o Sin2°^ 

becomes 

- €os2oi Sin2°C 

= 

Cos2<o Sin (-2=0 

Sin2 w 


- Sin2 to 


- Sin2 to 


V 
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v 2 


_ 3 — 

Cos2<uCos2 (-°C-j -°0 


Gos2wCos2 ( oC + °(l) 

Cos2 to Sin2 (- c < 1 -°0 

— 

Cos2 u Sin2 ( °f +oC^ ) 

- Sin2w 


Sin2 to 



Cos2oj Cos2 ( °C + 2o(j) 


Cos2w Sin2 (°C + 2°(]) 

1 

Sin2u> 



Where is the Stokes vector after the first 
half-wave plate with respect to its axes, tan w is the 
elliptic it y and is the azimuth with respect to the axes 
of the first plate, of the incident light ellipse, is 
the Stokes vector of light after the first half-wave plate 
with respect to the axes of the second half -wave plate, 
is the Stokes vector of light after the second half- 
-wave plate referred to its axes and is the Stokes vector 
after the second half-wave plate referred to the axes of the 
first half-wave plate. 

This shows that the effect of the combination of 
two half-wave plates is to rotate the light ellipse through 
twice the angle between the corresponding axes of the half.- 
-wave plates and in the same direction as the angle through 
which the first half-wave plate is to be rotated so that 
its axes become parallel to the axes of the second half- 
-wave plate. 
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Case (b)« Reversal of Light Path; 


For studying the optical phenomena when the light 

32 

path is reversed, the Jones vector is more convenient 
since it contains only two elements. 


The Jones vector V is given by V 


xT 

7 




The phase angle for the Y component is taken to 
be zero since we are only concerned with the relative phase 
difference and not the absolute phase angle in our present 
investigation. The transformation of the Jones vector 
through a retarder with retardation and azimuth *-j_ is 
given by 


V 


Cos'* - Sin* 
Sin* Cos* 

i£. 


i£i 


e 

o 


0 

1 


Cosc( Sine* 

- Sin* Cos <* 


M 


e ^ Cos. + Sin^*j Sin2°^ ^i^ _ 


) 


Sin2*i 


iG 1 

(e - 1 ) 


i£i 


Sln£5_ (e 1 + 1 ) 

2 


[V] = K,V 


(Both the Jones vectors V and V are referred to 
the same reference axes). 


T 

So, E| is a symmetric matrix, (i.e. ) 

It is of the form: 
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Let the matrix Mg corresponding to another retarder 
be given by 


M, 


&2 b2 
*>2 c 2 


So, Mg M 1 = M2 Mi = (Mi M 2 ) T 

Let M3 be the matrix representing a third retarder. 

«3 = “b 


T T T T 

Nov/, M 1 Mg = M 1 Mg = (M^ Mg M 1 ) 


For n retarders, we have 


M 1 Mg M- 


3....^ = wSJ l? 2 a (M n ....M 3 Mg M| )° 


Thus, the reversal of the light path is equivalent 
to transposing the matrix of the series of retarders. 

It is proved earlier, that a system of linear 
retarders is equivalent to a system containing a single 
retarder (of retardation G and azimuth * ) and a pure rota- 
tor (of rotatory power ijq ) . The transformation matrix for 
such a system is given by 


Cos'*-] 

Sin'*-] 


- Sin'*-] 
C os'*-] 


where *^ ~ <* + % 
iG 

e Cos'* Cos<j 
+ Sin* S±n*i 
iG 

: e Cos* Sin*.j 


- Sin* Cos**] 


iG 


r -i 

e 0 


Cos* Sin* J 

0 1 


- Sin* Cos* I 

— w 


Jf 


iG 

e Sin* Cos*-] 
- Since-] C os* ■ 
iG 

e Sin* Sin*.] 
+ Cos* Cos*,] 
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The rotatory power ^ for this system is given by 
tfi-j = °C-] - °c 

For a system consisting of a pure rotator (rotatory power 
) and a linear retarder (azimuth <*1, retardation £ ) in 
that order, the transformation matrix is again M-] * That 
is, the transformation matrix can represent either one of 
the above two systems. Now, the series of retarders M-j , 

M2 , My is equivalent to an optically equivalent 

system with a retarder and a rotator. Let this be repre- 
sented by the matrix « If the light path is reversed, 
the transformation matrix of the equivalent system will be 
the transpose of the matrix M-j i.e. M-j-. Now, M^ represents 
a system with a pure rotator of rotatory power - ^ = °c - oc-j 
and a retarder of retardation and azimuth <=( . 

Case (c). Reflection Polarjscope; 

In a reflection polariscope, the characteristic 
retardation p is added along the primary characteristic 
directions and the light ellipse, is rotated through an angle 
ip for the forward ray of light. When this is reflected, 
it is rotated in the opposite direction through the same 
angle \p and a retardation of p is again added along the same 
directions. Thus the combined optical path (forward and 
backward) is equivalent to a single retarder. (In this 
case, if the forward optical path is equivalent to a system 
represented by matrix M, the backward optical path will be 
equivalent to a system represented by the matrix M^ . ) 




TO. 10- DISTRIBUTION OF SHEAR STRESS t AND 

ISOCLINIC PARAMETER d> ALONG THE CHORD 
ABC IN AN AXI- SYMMETRIC CASE . 
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Case (d). Axi -symmetric Case; 

In this case, the magnitude and directions of the 
secondary principal-stresses are symmetrical about the 
middle point of the optical path in a transverse plane* 
(example; the optical path along a chord in a plane perpen- 
dicular to the loading axis of a sphere under diametral 
compression.) If the total optical path is divided into 
two equal parts, and if the first part is equivalent to a 
retarder (azimuth °f, retardation p) and a pure rotator 
(rotatory powers), then the second half is equivalent to a 
pure rotator (rotatory power -t|>) and a retarder (retarda- 
tion p and azimuth «<) . 

Now observe Fig. 10. The optically equivalent 
systems for both the light paths AB and CB will be same. 

Let the corresponding matrix for these optically equivalent 
systems be M. Then for the optical path BC, the matrix 
representing the optically equivalent system will be . 

So, for the total optical path ABC the matrix representing 
the optically equivalent system will be mFm. So this total 
optical path ABC is equivalent to a single retarder with 
azimuth °c and retardation 2p * The quantities °c and 2p 
however depend upon the rate of rotation of the secondary 
principal-stresses and the rate of retardation along the 
optical path ABC* 




FIG -11 . DISTRIBUTION OF THE DIRECTION <t> AND THE 

DIFFERENCE (p'-q') OF THE SECONDARY PRINCIPAL- 
STRESS ALONG THE OPTICAL PATH ABC • FIGURE 
SHOWS THE CROSS SECTION OF A RECTANGULAR 
BAR UNDER COMBINED TENSION AND TORSION 
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Case (e). Anti -symmetric Case; 

In this case the secondary principal-stress direc 
tions are anti -symmetrical and the secondary principal- 
-stress difference is symmetrical about the middle point o. 
the optical path (example; the optical path in a plane perp* 
dicular to the direction of tensile loading in a rectanguls 
model under combined tension and torsion.) 

In Fig. 11. let the optical path AB be equivalent 
to a system with a retarder (of retardation p and azimuth 
°0 and a rotator (rotatory powers). Then the optical path 
CB is equivalent to a system with a retarder (retardation p 
and azimuth -°0 and a rotator (rotatory power -tJj). The 
optical path BC is equivalent to a rotator (rotatory power 
if) ) and a retarder (azimuth°(- and retardation p). So the 
total optical path ABC is equivalent to a retarder (retar- 
datlon p and azimuth^), a rotator (rotatory power - 2<P) and^ 
another retarder (retardation p and azimuth -<=(). This 
system is again equivalent to a single retarder and a pure 
rot at or. The primary and secondary characteristic direc- 
tions will be inclined equally, but in opposite directions 
to the axes of the first and second retarders or to the 
secondary principal-stress axes at the points A and C . This 
can be proved as follows; 

From Eqs. (2.15) 

tan2w = tan£ Sin2o( I 

= tan (180 -C) Sin (- 2=0 and Cos2e = Cos2 <d Cos2o( 




FIG.12. EQUIVALENT SYSTEM FOR THE OPTICAL PATH ABC. 
OA,OB ARE THE AXES OF THE RETARDER Ri . 2y 
IS THE ROTATORY POWER OF THE ROTATOR. OA',OB' 
ARE THE AXES OF THE RETARDER R 2 -0R|,0'P 2 ARE 
THE PRIMARY AND SECONDARY CHARACTERISTIC 
DIRECTIONS OF THE SYSTEM RESPECTIVELY. 
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= Cos2© Cos(- 2=0 i.e. along axes inclined equally but in 
opposite directions (=(, - <=0 , to the major axis of a par- 
ticular light ellipse of ellipticity tan©, the phase diffe- 
rences will be 0 and 1 80 - 0 respectively. Also, the ampti- 
tudes of light vibration are equal along these axes. In 
Fig. 12. OA, OB are the axes of the retarder R-j 2$ is the 
rotatory power of the rotator, OA*, OB' are the axes of the 
retarder R 2 and OP^ , 0 1 P 2 are the primary and secondary 
characteristic directions of the system respectively* Now, 
let the incident light on this system be plane polarized 
along the primary characteristic direction 0P 1 . Then at the 
exit of the system, the light will be plane polarized along 
the secondary characteristic direction o'Pg. The changes 
in the properties of the light ellipse as it passes through 
the t .0 retarders will be as follows: 

The phase difference and the- amplitude ratio after 
the first retarder along its axes will be p and tan©j 
respectively. The light ellipse is then rotated through 
an angle 2^ by the pure rotator after which a retardation 
of p is added along the axes of the second retarder. 

Also, the light is plane polarized at the exit of 
the system along 0 1 P 2 , one of the secondary characteristic 
directions. So, the phase difference and the amplitude 
ratio of the light vibration before the second retarder 
along its axes should be ( it- p) and tan( - & 2 ) respectively. 
The same light ellipse has a phase difference of p and an 
amplitude ratio of tan e -j when referred to the axes of the 
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first retarder. 

So, 

tan 6^ = tan( -8 2 ) 
or 0^ = - 6 2 

Also, the angle between the axis of the first 
retarder and the major axis of this light ellipse is 

= (ill + °0. (Because along two sets of axes equ- 
ally inclined, but in opposite directions to the major 
axis of the light ellipse, the amplitude ratios are equal 
and the phase differences differ by ir . The effective 
(taking into account the rotatory power 2$) angle between 
the axes of the two retarders is (2 ^ + 2=0 . 

If ip +=< =s =<-1 

tan 2°c^ = tan2e^ Cosp - using Eqs. (2.15) 

or tan2( if> 4°c) = tan2©^ Cos p 

This relation connects the angle 6^ (the angle 
between the primary characteristic directions and the axes 
of the first retarder) and , =< and p . 

Thus, in the anti -symmetric cases, the primary 
and secondary characteristic directions are inclined equa- 
lly but in opposite directions to the secondary principal- 
-stress axes at the entrance and exit respectively. (Becau- 
se, the axes of the two retarders are inclined equally but 
in the opposite directions to the secondary principal - 
-stress directions at the entrance and the exit respectively). 
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Theoretical Determination of Characteristic Parameters; 

The matrix M in Eqs. (3.21) in the most general 

case will have three independent quantities namely the 

three characteristic parameters. In these equations, only 

2 2 2 _2 __2 _2 

two are independent since S-j + S2 + S3 = S-j + Sp + S3 = 1 . 

These are linear equations in Sq , S2) and • 

To determine these three characteristic parameters we need 
the exit light ellipse parameters S-j , S2 and S3 for atleast 
two sets of independent initial Stokes parameters S-^, S2 
and S3 so that we have four equations in the three unknowns. 
One of these three equations can be used for checking. The 
exit light ellipse parameters for two types of incident 
light ellipse can be found theoretically by integrating 
numerically the governing differential equations when the 
functions 0* and p’are known. From these values, the three 
characteristic parameters can be calculated. 

Experimental Determination of the Characteristic Parameters: 

The exit light ellipse parameters for two types of 
incident light can be found experimentally and the characte- 
ristic parameters can be calculated from this data. If a 
reflection polarise ope is used, the primary characteristic 
direction and the characteristic retardation can be found 
directly as the combination of the forward and (reflected) 
backward optical paths is equivalent to a single retarder. 
However, this can be used only for comparatively thin models. 
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These quantities can be found directly in a transmission 
polariscope for axi -symmetric bodies for, in this case also 
the optical path is equivalent to a single retarder*. 

It will be very useful, if there exist direct and 

simple experimental methods to find these parameters in the 

general case without any calculations or iterations. Robert 
14 . 15,16 

and Guillemet * ’presented a method of finding these 
parameters experimentally in a direct way. The principle 
of this method is the same as the one suggested by Kent 
and Lawson for the ellipticity measurement. In this 
method, the ellipticity of light from the analyser is 
calibrated in terms of the ratio of the constant and alter- 
nating components of light intensity, when the analyser is 
rotated at a particular angular frequency Z . In this case, 
the intensity of light from the analyser is given by I = 

(1 + Cos2u Cos2mt) where K 2 is a constant, Eqs. (2.15)* 
This method however requires additional instrumentation 
for measuring the alternating and constant components of 
voltage. To find the primary characteristic direction 
using this method, the polarizer at the entrance of the 
model is rotated until the ellipticity at the exit is gero. 

We shall now describe two direct and simple meth- 
ods to determine these characteristic parameters experi- 
mentally, which dbnot involve the rotation of the analyser 
at a particular speed or any additional instrumentation. 
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Method 1 : 


If the light at entrance to the model is circular- 
ly polarized, the light ellipse at exit will have its major 
axis at an angle of 45° to the secondary characteristic 
directions. This is so since the amplitudes along the pri- 
mary and secondary characteristic directions will be equal 
in this case (if the absorption of light along the optical 
path is neglected). If the light path is considered to be 
equivalent to a pure rotator and a retarder, the pure 
rotator will not have any effect on the incident circularly 
polarized light while the retarder increases the phase 
difference along its axes from 90° to (90 +5) where <5 is 
the characteristic retardation. Along these axes the amp- 
litudes will be equal even after the addition of retardation. 
These axes correspond to the secondary characteristic dire- 
ctions. This angle of ^90 +5) can be measured by using 
Tardy's method. One of the primary characteristic direc- 
tions is found by rotating the polarizer until the Inten- 
sity along one of the secondary characteristic directions 
is minimum. The final position of the polarizer axes are 
the primary characteristic directions.. 

Experiment al Pr oc edure : 

Keep the first quarter -wave plate axes at 45° to 
the polarizer axis so that we have a circularly polarized 
light at entrance to the model. With the second quarter- 
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-wave plate out of the field, rotate the analyser for mini- 
mum intensity and note this direction . This direction 
corresponds to the minor axis of the light ellipse. Along 
the directions + 45° the amplitudes of the components 
of the light vector are equal and these directions corres- 
pond to the secondary characteristic directions. 

Now, introduce the second quarter -wave plate so 
that its axes are parallel to the and D-j ± 90° directions . 
Rotate the analyser for minimum intensity which will be 
absolute minimum. Let the minimum intensity position be 
D . The 2 (D 2 - ) = (90 +8). Now remove both the quar- 

ter-wave plates and keep the analyser axes along one of the 
secondary characteristic directions. Rotate the polarizer 
for absolute minimum intensity. The final orientation of 
the polarizer axis corresponds to one of the primary chara- 
cteristic directions. 

Method 2 ; 

The principle of this method is the same as that 
of the intensity calibration method described in chapter 2. 

Procedure; 


Let the optical set up consist of the polarizer 
(at arbitrary orientation), the stressed model, the quarter- 
-wave plate (at arbitrary orientation) and the analyser. 
Rotate the analyser for minimum intensity thus finding the 
orientation of the major and minor axes of the light ellipse. 
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Rotate the analyser through b^° in either direction and note 
the intensity 1^ with a photometer. All ellipses ( of 
different orientations and ellipticities ) will have an in- 
tensity equal to V along directions which make an angle 
of 45° with their major or minor axis. Now keep the angle 
between the axes of the second quarter-wave plate and analy- 
ser equal to 45° and rotate the polarizer until the intensi- 
ty from the analyser is equal to the previously found value 
of 1^* The final orientation of the axis of the polarizer 
corresponds to one of the primary characteristic directions. 
The reasoning for this is as follows. When the intensity 
of light from the analyser is equal to V it means that 
the axes of the light ellipse are at to the analyser 
axes. That is, they are along the axes of the second quar- 
ter-wave plate. Since the phase difference along the major 
and minor axes is 90 ° t the phase difference along the same 
axes but before the second quarter-wave plate must have 
been zero which corresponds to the case of plane polarized 
light. The orientations of the polarizer axes which pro- 
duces plane polarized light after the model correspond to 
those of the primary characteristic directions. 

Now remove the second quarter-wave plate and rotate 
the analyser for absolute minimum intensity. The final 
orientations of the analyser axes correspond to the secon- 
dary characteristic directions. The characteristic retar- 
dation can be found as in the first method or by any other 


method 
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With minor modifications, the first method can "be 
used in the scattered light technique where the model acts' 
as an analyser and the second method in the technique where 
the model acts as a polarizer. These will be described in 
the next chapter. 

It may be noted now that in a transmission polari- 
se ope only the fractional part of the characteristic retarda- 
tion can be found directly while in general, it is difficult 
to find the integral part of the characteristic retardation. 
The usual method of using the Babinet-Soleil compensator with 
white light, as in two-dimensional problems can not be used 
here since the characteristic retardation in general, is a 
function of not only the wavelength of light, but also the 
orientation and magnitude of the secondary principal-stresses 
along the optical path. The importance of the concept of 
optically equivalent system and the characteristic parameters 
lies in the fact that in a three-dimensional model, in gene- 
ral, the preferential directions (directions which can be 
directly found by experiment) are the characteristic direc- 
tions and not the secondary principal-stress directions* 

In this chapter, the proof for the existence of 
an optically equivalent model is presented in detail. The 
behaviour of light when the light path is reversed is dis- 
cussed. The relationship between the two optically equiva- 
lent systems one with the retarder followed by a pure rota- 
tor and the other with a pure rotator followed by a retar- 
der, is derived. Using this information, the optical 
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phenomenon for each of the following three cases are discussed. 

(a) Reflection Polarise ope 

(b ) Model hawing axi- symmetric stress distribution 

(c) Model hawing anti -symmetric stress distribution 

Two direct and simple methods to determine experi- 
mentally the characteristic parameters are presented. These 
methods can be extended for use in a scattered light pola- 
riscope also. This will be discussed in detail in 
Chapter 



CHAPTER, if 


SCATTERED LIGHT METHODS IN THREE-DIMENSIONAL 

EHOTOELASIICITY 

Unlike the stress -freezing -technique, the scattered 
light method is non -destructive for analysis of stresses in 
three-dimensional phot oelasticity. In the scattered light 
technique, the models can he subjected to live loads at 
room temperatures and analysed, and no frozen stress model 
is necessary. In this chapter, we shall first make a brief 
survey of the existing techniques in scattered light photo- 
elasticity, investigate their limitations, and then suggest 
refined technioues. 

17 

It was first suggested by Weller that the three - 
-dimensional photoelastic analysis could be done by measur- 
ing the intensities of the scattered light components. His 
suggestion was based on the theory of absorption and re- 
-emission of light by vibrating particles. • ' 

According to this theory, (i) the light that is 
scattered at right angles to the direction of propagation is 
always linearly polarized (ii) the intensity of this 
scattered light is proportional to the square of its appare- 
nt amplitude when viewed from that direction. By 'apparent 
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amplitude’ is meant,. the component of amplitude normal to 
the line of sight. His method in the three-dimensional case 
consists of computing the secondary principal-stress differe- 
nce from the minimum spacing of the scattered light fringes. 
That is, the model is rotated about the light path under 
consideration until the spacing of the fringes becomes mini- 
mum. If T d' is the fringe spacing, the secondary principal- 
-stress difference (p 1 - q 1 ) is given by 

| | Tp 

(p - q ) = — where F is the material fringe constant 

( P si .in/fringe ) 0+.1 ) 

The directions of the secondary principal-stresses 
are found by rotating the model until the fringe system 
vanishes. This method doesnot take into account the effect 

of rotation of the secondary principal-stress axes. 

18 

Drucker and Mindlin investigated the effect of 
rotation of the secondary principal-stress axes along the 
wavenormal on the retardation in a photoelastic model. They 
obtained an analytical solution for the case where the prin- 
cipal-stress magnitude remains constant and their directions 
vary linearly along the wavenormal. They concluded that 
(i) the rotation increased the relative retardation by a 
factor S = (1 + 4R 2 )£ where R is the ratio of the rate of 
rotation to the rate of retardation and that (ii) at high 
stress, the light vector rotated with the secondary principal- 
-stress directions; plane polarized light entering along a 
direction of one of the secondary principal-stresses rema- 
ined practically plane polarized in the direction of the 
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stress through the entire model. They suggested that the 
factor 'S’ be incorporated in the scattered light analysis 
so that we have 

( p 1 - q ’ ) - — — (4.2) 

Sd 

37 

Saleme suggested a method of analysis which in- 
volves the measurement of absolute intensities in two sets 
of mutually perpendicular directions, the angle between 
which is 4-5°. His analysis is as follows: 

Let = (Si, S2, S3) be the Stokes vector of light referred 
to the secondary principal-stress axes at a particular point 
in the model. 

Si Cos2© 

V 1 = S 2 = Sin2© CosG (4.3)(a) 

S^ Sin26 SinG When the same light ellipse is 

referred to OXY making an angle °( with Opq, we get 

S-j Cos2© Cos2c< + Sin2S Cose SinSbC Cos2e^ 

V 2 = S 2 = Sin20 Cose Cos2=( - Cos26 Sin2=c = Sinl^ Cose-] 

S 3 Sin2© Sine Sin2e.j SinCi 


(4.3Kb) 


Let the amplitudes of light vector components along X,Y-axes 

a x ‘ 


be a x and a • tan©^ = ~t ^1 is the phase . difference corres- 


p 2 . — 2 

ponding to Stokes vector a^ - a^ = S^ = 2a x - 1 


1 - 2a: 


Y 


q 2 , a 2 _ 1 

Ciy^ *T* — * 


The intensities and along these .axes are given by 

. 0 1 + S-. 

I x = KaJ = K( 2 ~ i ) 


2 


(4*4) 
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1 4* S 

Iy = Ka| = K ( ^ — -) (4.5) 

where K is a constant. 

Now consider the intensities Iy' and Iy along the 
axes which are inclined at an angle of 45° to the X,Y-axes, 
Fig. 13 . These are obtained by substituting (°f + 45°) for of 
in the expression for in the above equations. 

So, Iy f = Jj (1 - Cos20Sin2 c < + Sin2© Cos£ Cos2=<) (4.6) 

I y ' = f + Cos2e Sin 2 of - Sin2e Cose Cos2=f) (4.7) 

Let z denote the distance along the wavenormal. 
Neglecting the rotation of the secondary principal-stress 
axes, i.e. 9-j and °f are independent of z, we have 

(4.8) 

(4.9) 

(4.10) 

(4.11 ) 

(4.12) 

from Eqs. (4.4) and (4.5) 

dl x dly 

d? Hz - 

aI X ' <% * 

dz dz 


dl 


X = 


dz 


- % Sin 2© Sin2of Sin€ ~ 
? dz 


dly 

~dz 


■5 Sin2 e Sin2°< SinG 
^ dz 


jS' = - ? Sin2e Cos2* SinG & 
dz 2 dz 


^ — Sin2© Cos2c< SinG ~ 
2 dz 


Hz 


Iy + Iy 


K 


tan2x 


(4.13) 
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from Eqs. (4.8) - (4.11 ) 


I x - I Y = KS 1 = K0os26 1 (4.14) 

X X' “ ly' = Sin20 1 Cos€ 1 (4.15) 


From the above Eqs. (4.12) - (4.15), the quantities 
K, =< , 9-j and £-| can be computed from the values of the 


intensities I-. r I v I v t and I v t . 

a, Y, X x 


Since (p 


T Hg 

q ; = F ■— 
dz 


( 4 . 16 ) 


The quantity (p - q ) can be calculated if the rotational 

effect is neglected. 
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Jessop attempted to arrive at a qualitative 
estimate of the magnitude of the rotational effect and 
deduced certain conditions under which this effect will be 
small enough to be neglected. He considered the Maxwell- 
-Neumann equations 


de _ d£ 
dz dz 


+ 


C 


fX - fY 
a Y a X 


) Sin 8 



(4.17) 


da X 

dp 


= a v Cos € 
x do 


and 


da Y 

dT* 


a Y Cos ^ 
x dp 


(4.18) 


Now, a relation connecting the quantity ~ (to be 
determined) and the quantity € (that can be observed experi- 
mentally) is to be arrived. 

The quantity --is equal to JJr at points where 
dz dz 

su. = a v or SinC = 0. The first condition can be 


either 
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satisfied at the point of entry of the beam into the model, 

by orienting the model so that the secondary principal- 

-stresses at this point are inclined at 45° to the direction 

of vibration of the incident light, but it is not in general 

satisfied elsewhere. The second condition is satisfied at 

points where whole or half-fringes appear but the value of 

at these points cannot be obtained from a smooth curve 
c 2 

plotted through the observed values. Also, in the absence 

of rotation, the quantity il£ will be equal to To find 

dz dz 

the effect of neglecting the second term in Eq. (4.17) which 
represents the rotational effect, we consider, 


d_ 

do 


/fX _ fYs 
'■ay ay' 


(2 + — 2 + — 2 ) Cos£4p 


a. 


*Y 

from Eq. (4.18) 


(4.19) 


^ 4Cos eff, if ay and ay are initially equal. 

If both 4^ and its derivatives are small 
dP 


I +CosC 


SinG 

d0dp 


4 (SinG $&) 
dP dG 


so that (If) _ ! % ^ 
Y a x 

fx_fy 
a y 


SinC 


Sin 0-^ — Sin 2 e 

dp d o dP 


where y 


d0 . d0 
d G dp 


frp, Eq. (4.19) 
„ . 2 

= y Sm £ 
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So, u is a small variable coefficient of the order 

An estimate of the magnitude of in any particular 
110 dp 

case can be made by observing the dimunition of the sharpness 

of the fringes in the model after rotating about the light 
path. If the total change in the orientation of the secon- 
dary principal-stress axes over the light path in the model 
is 30° and the relative retardation introduced (due to stre- 
sses only) is n wavelengths, the value of iHL can not exceed 

dp 

(1/12n). If the fringes are sharply', defined over a range 

of two or more wavelengths, the error introduced in by 

dz 

drawing the , z graph as a smooth curve through the points 
obtained by measurements of the positions of the fringes 
will not exceed something of the order of 1 % and this will 
generally be less than the possible error of measurement. 

If) 1] Ip 

Based on the above findings, Frocht and Srinath > } 

developed methods for measuring the secondary principal - 
-stress difference and their orientation. .An account of 
these methods will now be presented. 

Consider Eqs . (2.1) - (2.5) of chapter 2. Let 
X-and Y-axes represent the directions of the secondary 
principal-stresses. The intensity T* for a direction of 
observation making an angle °f with the X-axis is given by 
the expression 

I ^ = ay Sind* + ay Cos 2c * - 2ay ay Cose* Sin<* Cose (2.5) 
if the proportionality constant is assumed to be unity. 

This light intensity will be maximum when the 
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direction of observation is along the minor axis of the light 
ellipse and is a minimum when the direction of observation 
is along the major axis. That is, the intensity is a maxi- 
mum or minimum when the angle* satisfies Eq. (2.3) 

2 Uy ay 

tan2 0< = Cos C = tan2 0Cos £ (2.3) 

a X _a Y 

From Eq. (2.5), we note that for any fixed direc- 
tion of observation °(and for fixed values of ay and ay? the 
intensity of scattered light varies from a maximum to a 
minimum depending upon the value of £ . When °c is less than 
V2, the intensity is a minimum (although not necessarily 
zero) when £= 2mTr and becomes a maximum when £ = (2m - 1 )ir 
, When * is greater than ir/2, the condition for maximum 
and minimum will be reversed. Except for special cases, E 
will vary from point to point along the axis of propagation 
and in general, will assume the values necessary for a maxi- 
mum and a minimum light intensity. This gives rise to a 
sequence of bright and dark regions which forms the scattered 
light pattern. Also for approximately equal amplitudes 
ay and a.y, the axes of the light ellipse are at ^5° to the 
directions of the secondary principal-stresses and therefore 
this direction gives a pattern of sharpest contrast. 

If the directions of the secondary principal-stress 
axes donot vary along the optical path and the amplitudes 
of the light components along these axes at the entrance to 
the model are equal, (i) these remain equal all along the 
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optical path and. (ii) whatever retardation is added along 
these axes at the entrance, the same will appear as an addi- 
tional retardation at all points on the optical path in the 
model. We then have, 

(p' - q* ) = F — (4.20) 

dz 

A compensator can be placed at the entrance to the 
model with its axes parallel to the secondary principal-stress 
axes. If the axes of the polarizer, which is kept before the 
compensator make an angle of 45 ° with the axes of the compen- 
sator, the amplitudes along the secondary principal -stress 
axes in the model will be equal. The compensator then acts 
as a continuation of the model. Along directions making an 
angle of 45° with the secondary principal-stress axes, the 
scattered light fringes will have the sharpest contrast. 

This is so since, the major and minor axes of the light 
ellipse at any point on the optical path in the model will 
be along these axes. From these fringes, the n-z curve can 

be plotted from which the slope — can be computed. 

dz 

Even in the presence of rotation of the secondary 
principal-stress axes, the above method of plotting the 
n-z curve and thus calculating the secondary principal -stress 
can be used provided that the conditions mentioned by 
Jessop ^ 8 are satisfied. 

Directions of the Secondary Principal -Stresses : 

'■■■■■ . , ■ • ■ . ■ ■%. : 

Approximate Methods; 

The approximate locations of the secondary 
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principal -stress axes can be readily found by visual observa- 
tion without a compensator or a photometer. 

(a) Displacement Pattern Method: 

When looking along either one of the two secondary 
principal-stress axes, the center of the local scattered 
light pattern is of uniform brightness, and the fringes on 
opposite sides of this region become displaced relative to 
each other by half a fringe. The discontinuity in the frin- 
ges arises because, for a direction of observation on one 
side of the uniform region =< is less than y/2 and on the 
other side it is greater than ir /2. Therefore, the same 
retardation 2mir will produce dark fringes on one side and 
bright fringes on the other as can be seen from Eq. (2.5 )• 

The stress pattern resembles the picture of an X-ray of the 
spiral column with emanating ribs. 

(b) Sharpest Contrast Method: 

The model is rotated until the scattered light 
stress pattern appears of sharpest contrast to the naked 
eye when viewed along the line of vision of the photometer. 

In this position, the secondary principal axes are approxi- 
mately at + *+5° to the line of vision of the photometer. 

(c) Minimum Intensity Method: 

If the light is circularly polarized at the entrance, 
the axes of the light ellipse at any point are independent 
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of the birefringence and are always inclined at .4-5°, approxi- 
mately to the directions of the secondary principal-stresses. 
Hence, the determination of one of the axes of the ellipse 
also determines the directions of the secondary principal- 
-stresses. The direction of the major axis of the ellipse 
can be readily and accurately located with a photometer 
since, the intensity of the light is a minimum for observa- 
tions along the major axis. 

(d) Exact Methods: 

(i) The Constant Intensity Method: 

When is small, the changes in the amplitudes 
dz 

ay and a^-, when retardation in the compensator (whose axes 
are set parallel to the directions of the secondary princi- 
pal-stresses) is varied, are negligible. So, from Eq. (2.5), 
when ** = 0, the intensity 1^ is independent of the bire- 
fringence and it remains constant even though the retarda- 
tion at that point is changed. 


( ii) The Zero Intensity or Line Ellipse Method: 


From Eq. (2.5), it can be observed that for a fixed 

value of* ( <ir /2) , the intensity is a maximum when e= (2m - 1 ) 

* and a minimum when £ = 2mTT . In both these cases, the 

light ellipse degenerates into a straight line. When 

-1 . a X\ 

6 = 2mir , the line is inclined at an angle ^ = tan (~) 

and when' £ = (2m - 1 ) T , the line is inclined at an angle 
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at -1 7 a Y\ 

= tan (g^) to the X-axis. The intensities of scattered-* 
light for observations along these lines will be zero* To 
determine the directions of the secondary principal-stresses, 
the retardation in the compensator is adjusted to make the 
intensity at a point minimum for an arbitrary direction of 
observation say OM, Fig. 1*f. Assuming * to be less than 
* rr /2, £ will then have assumed a value of 2mir and the light 
ellipse will have become a straight line CO, Fig. 14. making 
an angle of = tan~^~) with the X-axis, which can be located 
with a photometer. Next, keeping approximately the same 
direction of observation OM, the retardation in the compen- 
sator is adjusted to give the maximum light intensity. When 
this occurs, the light ellipse will have become once again 
a straight line, only this time it will be DO inclined at 

an angle = tan (~0 to the X-axis. This line is also 

a X 

located with a photometer. The directions CO and DO, are 
equally inclined with respect to the direction OX which is 
one of the secondary principal-stress axes. 

(iii) Absolute Intensity Method; 

This method which involves the measurement of the 

39 

absolute intensity has been suggested by Menges and 
Salame 37 , This has been described earlier. The main dis- 
advantage with this method is that it involves the measure- 
ment of absolute intensity which will depend upon the absorp- 
tion along the varying optical paths. 




FIG. 15. SCATTERED LIGHT TECHNIQUE USING THE 
MODEL AS A POLARIZER, 
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40 

Cheng described a new method which records simul- 
taneously the intensities of scattered light along two direc- 
tions of observation 1 +5° apart in a plane normal to the beam 

of light. The theory of this method is the same as the one 

37 

described by Salame except that the effect of varying 
lengths of light path is eliminated by recording the inten- 
sities before the loading of the model. He made use of 
Drue her and Mindlin’s conclusion that for cases where the 
parameter R is small, plane polarized light entering the 
model along one of the secondary principal-stress axes at 
the entrance follows that axis throughout the model. The 
polarizer axis is set along this direction and both the pola- 
rizer and the model are rotated in unision until the light 
intensity is minimum at the particular point under conside- 
ration. When the minimum occurs, the direction of observa- 
tion is along one of the principal-stress axes. 

Robert and Guillemet siHbpted an altogether 

different approach towards scattered light analysis in the 
three-dimensional cases. Since any length of light path 
in a photoelastic model under load can be replaced by the 
combination of a linear retarder and a pure rotator, the 
characteristic parameters of this equivalent system can be 
found at each point on the optical path with respect to the 
point at the entrance to the model, Fig. 15* After knowing 
these parameters, the secondary principal-stress difference 
and orientations can be found by using the governing optical 
equations. To find these parameters, they adopted the method 



that is popular in crystal physics, of rotating the analyser 
at a particular speed, and. then measuring the ratio of alter- 
nating and constant components of light from the analyser. 
This ratio is dependent only upon the ellipticity of light 
from the analyser (see Chapter 3) • The light incident on the 
model in this method is unpolarized and the model acts as a 
polarizer. The intensity from the analyser is given by 

I = K (.1 1 51 ) 

2 

where is the first Stokes parameter of the light 
ellipse and K is a constant, 
or I = | (1 - Cos2w Cos2=0 

Tan o) is the ellipticity and is the azimuth (with 
respect to the analyser axes) of the light ellipse incident 
on the analyser. If the analyser is rotated at an angular 
velocity of u , = co t (t is the time.) The ratio of the 

alternating to the constant components of this intensity is 
obviously 

1 - e 2 

C os2 w = J 

1 + e^ where e = tan w 

Also Cos2w = \/\ - S 2 = - Sin 2 2 © Sin^£ 

= Y ( say) 



then Cos 6 . 
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To find the primary characteristic directions the 
model is rotated with respect to the incident light -vibra- 
tion -until the value of Y observed is zero, which corres- 
ponds to plane polarized light (e = 0) before the analyser. 
The position of the incident light beam with respect to the 
model axes will be the primary characteristic direction. 

Now this beam is rotated through an angle of ^5° and the 
value of Y is noted. Then Y = Cos£ where £ corresponds to 
the characteristic retardation. The direction of the major 
axis of the light ellipse before the analyser when the inci- 
dent light beam is along one of the primary characteristic 
directions corresponds to cne of the secondary characteris- 
tic directions. For finding these parameters the model is 
rotated about the axis 00* , Fig. 15« The§e parameters are 

found for the optical paths OMLj , 0M 2 , .... 0^,.*.., 00* • 

XS 

Aderholdt, Mckinney, Ranson and Swinscn in their 
scattered light method assumed that for problems where the 
rotation to retardation ratio is small, the light waves 
follow the secondary principal-stress axes. That is, the 
amplitudes of the components of light along these axes re- 
main constant along the light path. This assumption is 

lo 

again based on the findings of Drucker and Mindlin . Keep- 
ing the light at the entrance to the model circularly pola- 
rized, the. directions at each point of the optical path in 
the model along which the fringes appear to be fading (or 
the intensity is constant) are noted which give the directions 
of secondary principal-stresses. The fringes are counted by 
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moving along the light path and at the same time, rotating 
the observation angle such that the angle of observation at 
any point is 2 +5° from one of the secondary principal-stress 


axes . 
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Srinath made an analysis of some of the previous 
scattered light techniques. He suggested a method called 
mini -max method, which is a combination of the zero and 
absolute intensity methods, to find the directions of the 
secondary principal-stresses. In this method, the orienta- 
tion of these directions is found from the formula 

,J 1 + y 2 


tan * 


v- 


l 3 + Vi}. 


where 


2 2 
v-j = K (a Sin 1 * + b Cos*) = (I*) 


y 2 = K (a Sin* - b Cos*)*" = (I*) 

^ ) • 

v 2 = K (a Cos* + b Sin*) 2 = (Ip) 


= K (a Cos* - b Sin*) 2 = (Ip) 


. rr 
*+ -jy 


and the retardation p is obtained from the equation 


max 


min 


max 


mm 


Cos p = 


-2L + (V-j +v 2 ) 


(vi - V2) 


where 


I* = K (a 2 Sin 2 * + b 2 Cos 2 * - 2ab Sin* Cos* Cos p) 
- obtained from Eq. (2*5) 
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Critical Analysis ; 


A critical analysis of the above described scattered 
light methods will now be made. Except the method described 
by Robert and Guillemet, all the methods are. connected with 
the findings of Drue her and Mindlin ^ . These findings in 
turn are based on the analytical solution for the case where 
the secondary principal-stress difference is constant and 
their orientation varies linearly along the optical path. 
These results are confined to the cases where the parameter 
R is small. 

We shall now examine the validity of the assumptions 
involved in the previous scattered light methods for the case 
where the rate of rotation 0' can be any arbitrary function, 

but the ratio R = ^ is a constant. 

p* 

The analytical solution for this case is given by Eqs. (2.46) 

- (2.48). These are: 

A 1 1 . 




2RK‘ 


~ + \ (1 - -4;)^ Sin ( 2K0 + Ag) 


A 2 a 

S 2 = O - -p)^ Cos(2K0 + Ag) 


and 


_ Si 

S 3 = 2R + A 1 


Al 

K2 


(1 -4 

__Llf 

2RK 


Sin (2K0 + & q ) 


(4.21 ) 

(4.22) 

(4.23) 


where A-j are arbitrary constants which are to 

be determined from the initial conditions. 

Let us now examine this solution for two types of 
incident light, (i) plane polarized light along cne of the 
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secondary principal-stress axes at the entrance and (ii) 
circularly polarized light. 


C ase (i) ; Incident Light is Plane Polarized Along one of 
the Secondary Principal -Stress axes At the Entrance 


At Z = 0, 0 = 0 (say), S-j = 1 , S 2 = S3 ± 0 


From Eq. (4.23) 

A 1 = - 1/2R 
and from Eq. (4.21 ) 

1 = — 1 +1 0 1 )i 

4« 2 + 1 K 4H a + 1 


SinAg 


X.e SinAg = 1 
or A 2 = it / 2 . 

So, at any point on the optical path we have 

% = (1 '^r—/ Sln (ZKI> + 2> 

= — ~ (Cos2K0 - 1 ) 

2RK 2 

The light is plane polarized at points where 

a 3 = 0 

i.e. when 2K0 = 2nir 
At these points 

g _ - - + i (1 - - — i— 75)^ 

4r 2 K 2 K 4r 2 K 2 

= 1 


( 


(4.24) 
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i.e. the intensity attains its absolute maximum ’value. So, 
the light vector will be along the direction of one of the 


secondary principal-stresses. 


The change in the orientation of the secondary 


principal-stress axes between two fringes is given by 
Eg. (4 . 24- ) as 

d0 = tt/K = Ido 


(4.25) 


i.e. dp = "ff/m 

i.e. between two fringes, the change in the retardation dp 
introduced due to the stresses only is given by Eq. (4.25). 

If the effect of the rotation of the secondary principal- 
-stress axes is neglected, this is given by 
dp = 2 TT 

So, the effect of rotation is to decrease the fringe 
1 1 

spacing by a factor — — or , 

2KR (4r 2 + i )2 


Gase (ii) ; Incident Light is Circularly Polarized. 


At z = oj s 1 = s 2 = o; S3 1 , ft = 0 

From Eq. (4.23 ) 

M = 1 

and from Eq. (4.21 ) 

0 = ' SF + i (1 ' F )1? SinA s 

i.e. SinA 2 = 1 or = ir / 2 
At any point on the optical path, 

4R 2 Cos2K# 

3 4h 2 + 1 + 4r 2 + 1 

= f 
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if Cos2K0 = - 4R 2 

If4R 2 «1 

i.e. 2K0 = (2n + 1 ) 
At these points, 


S. 


2RK* 


+ ~ (K 2 -1)V 2 Gos2K0 


1 1 

■ 9 + o Cos2K0 

2RK?' 2RK 2 


1 


(1 - Cos2K0) 


2RK^ 

2R (1 + 4R 2 ) 
(4R 2 + 1 ) 




0 . 


S o ~ -'(1 


- Sin2K0 


(4.26) 

(4.27) 


1 

(4r 2 + i )i 


Sin 2K0 


- Sin -^4-±— 1 ) . 7r 
2 


(4.28) 


The second Stokes parameter gives us the inten- 
sity along axes making an angle of 45° with the reference 
axes (in this case, the secondary principal-stress axes). 

From Eqs. (4,27) and (4.28) we observe that between 
two successive fringes, 


d0 = Rdp = ^ 


(4.29) 


So, in this case also, the effect of the rotation 
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of the secondary principal-stress axes is to decrease the 

A 

fringe spacing by a factor — — x . The same conclusions 

(4R 2 + 1 )? 

to 

were dram by Drucker and Mindlin for the case where both 

of the functions 0' and p'are constants. 

Now, we shall examine the above analytical solution 

in light of the assumptions of the previous scattered light 

methods. In the scattered light method developed by Frocht 

10 , 11.12 

and Srinath , the following assumptions are made: With 

the addition.' of a certain amount of retardation A£q at the 
entrance along the secondary principal-stress axes when the 
amplitudes of light vector components along these directions 
are equal, (1) a retardation of aG q appears as an additional 
retardation at every point on the optical path along the 
secondary principal-stress axes and (2) the amplitudes (or 
the intensities) along the secondary principal-stress axes 
at any point on the optical path do not change with the 
change in the retardation G^ added at the entrance. At 
Z = 0 let S-| =0 (i.e. the amplitudes are equal along the 
seccodary principal-stress axes at the entrance). 

= Sin£ 0 , 0 = 0 ± , S 2 = Cos^q 


The light intensity I = K a- 


'X 

_ 1 + S, 

K( — - — - 


) 


where K is a constant. 


dl _ K 
dG^ 2 dG 0 


0+.3O) 
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Let us evaluate the constants A^ and A 2 of the 
solution for these initial conditions. 

From Eq. (4.23) 

S 3 = Sin€ Q = a^ (4- *31 ) 

From Eq. (4-. 22) 


Sin^S 


0 . 


S 2 = Cose 0 = (1 - Cos(2K0 ± + A 2 ) 

K 


or C os ( 2KjZ5^ + Ag) = Cos £ 0 


(4.32) 


(1 - 


Sin^S^ 


1/2 


and from Eq. (4.21 ) 


Sin(2K0j_ + A 2 ) 


SinG, 


D_ 


Sin^ 0 , 
2HK(1 - — — )i 
K 2 


(4.33) 


Substituting these expressions in Eqs. (4.21) - (4.23) so 
that the constants A^ and A 2 are eliminated, we obtain the 
following expressions for S-j , S 2 and S 3 at any point. 


Sine 


Si = - 


0 1 


2 R.r 


s 2 = CosK 2 Cose Q _ 


~ + ^ (SinK 2 CosCq + 
SinK 2 Sin0Q 


CosKL, Sin£ 


0 


2R.K 


) 


2RK 


and 


SinE, 


S 3 = ^ ° + (SinK 2 Cose 0 + 

where Kg = 2K(j3 - 0 ± ) 


CosKg SincQ 
’ 2RK 


(4.34) 

( 4 . 35 ) 

(4.36) 

(4.37) 
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dCE _ _ "&R 
& 0 (1 + 4R^) L 

- CosK 2 Cose^j 

= 0 

for all ■values of 0 


[cose 0 + SinK 2 Sin0 o (1 + 4R 2 )i 


■0 


when 


K 2 = 2m = 2K(0 


0±) 


(4.38) 


Also 


ClI <rsJ 


dS, 


0 


0 


for all values of £ 


0 


when R«1 (4.39) 

The phase difference 0 along the secondary principal-stress 
axes is given by 

Sin0Q SinK 2 Cos£q CosK 2 Sin£o 


K c 


2RK 


4R 2 K 2 


tan€ 


h_ 

S 2 SinK 2 Sin£ 0 

CosK 0 Cos€ n 

2 0 2RK 

4r 2 Sin0 o + (4R 2 + 1)£ SinK 2 Cos% + CosK 2 Sin% 
(1 + 4R 2 ) CosK 2 Cos^q - (1 + 4R 2 )'S' SinK 2 Sin€^ 


tan€, 


0 


if = 2nir 


(4 .40 ) 


tan£ tan(0Q + K^) 


(4.41 ) 

^ £ 

if 4R^«1 

The quantity 4R may be neglected when compared to 
unity if it is less than 0*01 i.e. for R<0.05» 

So, the assumptions in the above scattered ligljt 
method may be considered valid for the cases where the 
parameter R is less than 0.05* 
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13 

Aderholdt, Mckinney, Ranson and Swinson , in their 
scattered light method assume that the intensity along the 
secondary principal-stress axes remains almost constant (for 
certain class of problems where the rate of rotation is not 
high) along the optical path. 

Also from the above solution, i.e. Eq. (4.3*0, 

'0 


daf = . = 0 , ( flnK, Sin£, 


dz 


2RK 


j - CosK^ CosCq) 


<20 t 

i.e. da? <2d0 (4.42) 

2 2 2 2 
(where - 3 % ~ a Y an 4 a X + a y = ^ •) 

i.e. the intensity variation along the secondary principal - 
-stress axes depends upon the rate of rotation and not on 
the parameter R. The change in the square of the amplitude 
(and so the intensity) will be less than 0.1 if the rate of 
rotation is less than 0.05 radians per one tenth of a milli- 
meter. The above scattered light method may be taken to be 
valid for rates of rotation lesser than this value. Even if 
these assumptions are satisfied at each point, the errors 
involved may be cumitlative and prove to be disastrous over 
long light paths. This sets a serious restriction on these 
methods. Although the above remarks are based on the solution 
for a special case where the ratio R is a constant, this 
solution holds over a length of light path even in a very 
general case. This makes the above remarks relevant. 
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38 

Jessop from his analysis concludes that the error 
introduced by drawing theg , z graph as a smooth curve through 
the points obtained by measurements of the positions of the 
fringes will be y Sin 2 g where y is a small variable coeffi- 
cient of the order In his derivation he assumes that 

dP 

both the quantities and its derivatives are small so 

dp 


that 


Cos 


dp 


— (Sine -£) 

dp d£ 


Sin£ 


d 2 0 

d€dp 


(SinC M) 
dp dp 


Analysis as Applied to a Model Under Combined 
Tension and Torsion : 

Figs. 16-25 show the plots of various quantities 

for the case of a model under combined tension and torsion. 

The model material is paraplex. The material fringe constant 

is 139 Psi. in. /fringe. The values of shear stress are com- 

42 

puted from the series solution, Timoshenko . The integra- 
ted optical effects are computed by numerically integrating 
the governing optical equations in terms of Stokes parameters. 
Runge-Kutta-Gill method of numerical Integration is used. 

It was found that the results didnot vary when the number 
of steps is increased from 50 to 100. At each step the 

values of the three Stokes parameters S-| , S 2 and S 3 were 

2 2 2 

found to satisfy the equation S>| + S 2 + So = 1 . 
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From Figs. 16 and 17 ve observe that the quantities 
R and^-S are of the same order of magnitude along most of the 

d£ OS 

length of the light path in the model. Hence Jessop’s 

assumption of neglecting the term Sin£ compared to the 

at 

term RoCQ5£ is not justified in these cases. It can be said 
■that in general, the quantity need not be small and negli- 
gib le. 


The conclusions of Drucker and Mindlin 18 are based 
on the solution for the case where 0' and p' are constants 
and their ratio is a small quantity. However, these con- 
clusions cannot be extended to the general case as can be 
seen from Figs. 19 -21. Fig. 19* shows that the phase di- 
fference along the secondary principal-stress axes calculat- 
ing without considering the rotation is greater than that cal 
culated with rotational effect, ’along some part of the opti- 
cal path. This is contrary to the prediction by Drucker 
and Mindlinl 8 

Now consider the assumptions in the methods descri- 
bed by Frocht and Srinath 10,11, . 12 These are not satisfied for 
the particular loadings considered as can be seen from Figs. 
22 and 23. The variation of the amplitudes along the secon- 
dary principal-stress axes with -distance z along the light 

path is quite marked for these cases, Fig. 18. and the method 

13 

described by Swinson et . al is not applicable here. So, 
before using either of these methods, one has to check for 
the problem under consideration whether the rate of rota- 
tion 0* or the parameter R is within the limits suggested in 



24. VARIATION OF THE SEC. PR. STR. AXES ORIENTATDN 
4 ALONG THE OPTICAL PATH IN THE MODEL. 









115 


the previous section. The plot of 0 Vs. z and (p * — q’) 

Vs. z are presented in Figs. ^+. and 25. for reference. 

-|Lj. 1 5 -[ 5 

The method described by Robert and Guillemet- ’ y 1 
doesnot involve any of these above mentioned assumptions. 

But the analysis here is through the characteristic para- 
meters. This method requires the rotation of the analyser 
at a particular speed and some additional instrumentation. 


New Methods; 

It is highly desirable and useful to have methods 
which donot have any of these drawbacks. These methods should 
essentially be based on the Maxwell -Neumann equations or 
their equivalent forms since these are the equations availa- 
ble at present which describe the optical phenomena in the 
presence of the rotation of the secondary principal-stress 
axes, accurately as sho-wn by various dynamical derivations 
(Chapter 2). It is now clear that at any point on the optical 
path in a photoelastic model under three-dimensional state 
of stress in general, the preferential directions (the dire- 
ctions which can be determined directly by experiment) are 
the characteristic directions and not the secondary princi- 
pal-stress directions. So, in general, the secondary princi- 
pal-stress axes cannot be found directly from the experimen- 
tal at a single point unless it is coupled with the data 

at the neighbouring points. 

The scattered light technique involves the follow- 
ing steps: (1) measurement of scattered light intensity 
along different directions taking into account, the variation 
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of the lengths of the light path involved, (ii) finding the 
parameters of the light ellipse and the characteristic para- 
meters at the point under consideration and (iii) evaluating 
the functions 0 and p using these values and the governing 
optical equations. 

One of the methods of avoiding the error due to the 
variations of the length of optical path is to obtain absolu- 
te minimum (zero) intensity at the point under consideration 

in 11 ip 

along any particular direction (Frocht and Srinath ’ f ). 
Another method is to note the intensities along the required 
directions before and after the loading of the model (Cheng ^) . 
The former method has the drawback in that, it would not be 
possible to obtain the zero intensity in general, using the 
methods described previously because of the rotational 

14. TK 

effect. In the method described by Robert and Guillemet ’ ’ 

the variation in the length of light path doesnot enter into 
picture. 

Now, let us recall that (i) the amplitudes along 
directions at 45° inclination to the major axis of a light 
ellipse on either side are equal and (ii) any part of the 
light path in the model is equivalent to a system consisting 
of a linear retarder and a rotator. Now consider the passage 
of light through a system consisting of a linear retarder 
(retar dance p-j , azimuth of ^ ) and a pure rotator of rotatory 
power . The transformations through this system are 
given by 
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f -t 




j S 1 


S 1 


8iCos2^ + Sin 2 ^(S 2 Cos Pi 





- S 2 Sinpi ) 

OJ 

m 


S^Cos?^ - S^Sinoi 


(S 2 CosPi - S^Sinpi ) Cos 2 b 





- S.,Sin2ip 

k S 3 


S-,Cospi + S p Sin pi 


S 3 C 0 S pij +S 0 Sinpi 

„ _ 

i 

r__ n 

- 

- J 




=3 


(4.lf2) 


V = (S-| , Sg , S 3 ) is the Stokes vector of the incident 
light referred to the axes of the linear retarder V = (1L 3 
S ' 2 } S 3 ) is the Stokes vector of the light at the exit of 
the system referred to some arbitrary set of axes $ = (of - ^ ) 

! j 

whereof is the angle between the set of arbitrary axes X, Y 
chosen and the axes of the linear retarder. 

Also , 

- 2 2 

Si = - ayi where , a y i are the amplitudes of 

the components of light along the X f , Y ( axes. 

So, 

_ o 2 

- aj 1 = SiCosS^ + (S 2 CosPi - S^Sinp^ ) Sin2 ^ 

2 


5 1 


Consider now, the following cases; 


1 ) 

(a) 

3 1 

= S 2 = 0, S 3 = 1, S-] = 

a X ! 

- ay’ 

= - Sin pi Sin2 i|» 


(b) 

S 1 

= S 2 = 0, S 3 = -1, Si = 

2 

a X’ 

2 

' a Y 5 

= Sin Pi Sin2 ij> 

2 ) 

(a ) 

S 1 

= Cos2g, S 2 = Sin2p, S 3 : 

= 0 , Si = 

(Cos2 p Cos2 v 


+ Sin2P Cos Pi Sin2\0 
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(b) S 1 = Cos2( p + J), S 2 = 3in2( p + |), S 3 =0, 

S-] = - Cos2P Cos 2'p - Sin2B Cos?-j S in 2 s 

3) (a) S>j = 0, S 2 = CosG, S3 = SinG, = Sin2<KCos£>[ Cosp-j 

- SinG Sinp-j ) 

= Cos(G + P-j ) Sin2'4' 

(b) S 1 =0, S 2 =Cos(G + 7 t), = Sin(e + 1 r) 

S>] = - Sin2q Cos( £+ p-j) 

Vie note that by changing the incident Stokes para- 
meters from the case (a) to those of case (b), the sign of 
the first Stokes parameter of the light at exit is reversed. 
This means that, if the amplitudes along the X 1 axis is ax’ 
in case (a), it will become ay' in case (b). This is true 
for any direction of observation since the angle of is taken 
to be arbitrary. We are thus able to obtain the component 
ayf along the axis X'. This can be achieved in many more 
ways. But the above mentioned ways can be said to be repre- 
sentative of the rest. Of these again, the first two cases 
(1) and (2) are of inport ance since these donot require the 
knowledge of the primary characteristic directions. In case 
(1), right circularly polarized light is converted into left 
circularly polarized light at the entrance. This can be 
achieved by introducing a half-wave plate or by rotating 
the quarter -wave Plate after the polarizer through 90°. In 
case (2), the polarizer is rotated through 90°. Using these 
respite,, we hj.ve 0 the 0 following . sqqtjt ered... light techniques^ 



S - Light source( laser) 

P - Polarizer 
(X/4)- Quarter - wave plate 
EM - Direction of observation 


f E 


FIG. 26- SCATTERED LIGHT TECHNIQUE WITH THE MODEL 
ACTING AS AN ANALYSER. 
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(1 ) The Model Acts As an Analyser 

A schematic diagram of the set up is shown in 

Fig. 26. 

Method (a) 

Let the incident light be right circularly polari- 
sed. Note the intensity 1^ along the direction of observa- 
tion EM at the point under consideration M. Now make the 
incident light left circularly polarized and note the inten- 
sity I 2 * Now change the direction of observation until 

II = I 2‘ 

If Lj corresponds to an amplitude of a then I 

2 2 

corresponds to where ay + ay = 1 . When, the condition 
a X ~ as stained, it means that the direction of obser- 
vation is making an angle of 45° with the major axis of the 
light ellipse at the point under consideration. So the 
final position of the direction of observation will then, 
correspond to one of the secondary characteristic directions 
(for the light path.OM). Now keep the direction of observa- 
tion in this position, rotate the polarizer (after removing 
the quarter-wave plate) until the intensity of light becomes 
zero. The final positions of the axes of the polarizer are . 
the primary characteristic directions. Now, keep the axes 
of the quarter -wave plate 45° to the primary characteristic 
directions and the direction of observation 45° to one of 
the secondary characteristic directions. Rotate the 
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polarizer until the intensity observed is zero. Then the 

characteristic retardation will be equal to 2°C where is 

the angle between the axes of the polarizer and the quarter- 

43 

-wave plate (inverse Tardy method, Cheng ). This is so 
since the zero intensity (corresponding to plane polarized 
light) will be obtained at the point under consideration 
only if the retardation in the linear retarder of the equi- 
valent system is nullified at the entrance. 

The total optical path in the model can be broken 
into small segments so that each segment can be represented 
by a single linear retarder. The axes of these retarders 
can be taken to be the directions of the secondary principal- 
-stress axes at the middle points of the corresponding seg- 
ments. The secondary principal-stress difference can be cal- 
culated from the retardance. of this linear retarder. Consi- 
der the segment AB of the optical path, Fig. 26. First find 
■the characteristic parameters for the optical path OA as 
described above. Keep the axes of the quarter-wave plate 
at ^+5° to the primary characteristic directions. Keep the 
angle between the axes of the polarizer and quarter-wave 
plate as (? r A -°</2) so that we have circularly polarized 
light (say right handed) at the point A. Note the intensity 
1^ along any direction at the point B. Now change the light 
at the point A from right handed to left handed circularly 
polarized light by rotating the polarizer through 90° a nd 
then note the intensity Ig. Change the direction of obser- 
vation until I 1 = Ig. The final direction of observation 



122 


corresponds to one of the secondary principal-stresses in 
the region AS. To find the retardation corresponding to the 
segment AS, nullify the effect of the equivalent system for 
the optical path OA. Let the characteristic parameters of 
this equivalent system be P-j (primary characteristic direc- 
tion), p (characteristic retardation) and ^ (rotatory power). 
Now keep a compensator at the entrance with its axes along 
these primary characteristic directions and with a retarda- 
tion of -o . Keep the axes of the quarter-wave plate along 
a direction making an angle of (v/4 - $ ) with secondary 
principal-stress directions of the segment AE. Fix the dire- 
ction of observation along a direction making an angle of 
v/lf with these secondary principal-stress axes. Now rotate 
the polarizer until the intensity is minimum. If the angle 
between the axes of the polarizer (in the final position) 
and the quarter -wave plate is o (. , then the retardation corr- 
esponding to the segment AB is - 2=< . 

The orientation and difference of the secondary 

principal-stresses can also be found from the characteristic 

parameters using the governing optical equations referred 

to some fixed axes, Eqs. (2.42) - (2-44). For a particular 

incident light the three Stokes parameters Sq , Sp and S 3 

with respect to some fixed axes are found at each point 

using the characteristic parameters determined experimentally. 

The derivatives S^ , Sg and S 3 (with respect to the distance 

z along the optical path) are also found. Then from Eqs. 

(2.42) - (2.44) si 

tan2jZ5 = - -7 
b 2 




FIG. 27. SCATTERED LIGHT OBSERVATION ALONG 
TWO SETS OF AXES 45° APART. 
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and 

S 1 

D * = SrjSin20 

(b ) Determination of Characteristic Parameters From the 
Measurement of Intensities Along Two Directions 
Apart 

Consider any point £ on the optical path in the 
model. Choose two directions which are ^5° apart and denote 
these by the numbers 1 and 2, Fig. 27» let the incident 
light be circularly polarized. Note the intensities I 1 and 
I 2 along these directions. 

We then have 

1 1 = K-| af 

1 2 = K 2 a| 

where K, K are constants which take into account 
1 » 2 

the absorption through the corresponding lengths of scattered 
light paths. Now make the incident light left circularly 
polarized and note the intensities 1^ and 1 2 along the dire- 
ctions 1 and 2. 

I.j = K,(1 - af ) 

= K, aj 

l ' 2 = Kgd - a 2 ) 

= 


and 
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So I 1 + I-J = It, 


and 



Now consider the Stokes vector transformation when the 
reference axes are rotated through 4-5°. 



i.e. in this transformation the first and second Stokes 
vectors are interchanged with a sign change. 

Also 

2 2 2 
S 1 + S 2 = Cos 2 to 



Cos^2o> 



tan2°c 


d 2 - ig) di + ij ) 
(I 2 + I 2 ) (Lj - l\) 


The characteristic retardation p of the equivalent 
system (for the optical path from the point of entrance of 
the beam into the model to the point P) is given by 
p = ( 2<c) - tt/2 ) 



i 



FIG. 28. SCATTERED LIGHT TECHNIQUE WHERE THE MODEL 
ACTS AS A POLARIZER. 
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The orientation of the secondary 

= =< + t r /+• 

characteristic direction 

To find the primary characteristic directions rotate 
the polarizer (with the quarter-wave plate out of field) at 
the entrance to the model until the condition 


„2 a 2 i£i 

D-j 4- S2 = i 


II 


r2 


: < 1 O 

il 2 - 12! 


I 1 + I 




T 2 + l2j 


= 1 


is satisfied. 

The final positions of the axes of the polarizer correspond 

to the primary characteristic directions. Alternatively, 

keep the angle between the axes of the polarizer and the 

quarter -wave plate equal to (v /b - p/ 2) and rotate them in 
I 1 I 2 

unisicn until — = (£, and K 0 are already known). When 

Ki K2 

this condition is satisfied, the light ellipse at the point 
P is circularly polarized and so the final directions of 
the quarter-wave plate axes will be at *+5° to the primary 
characteristic directions. 

In Appendix A are described two more scattered 
light techniques, (i) determination of characteristic para- 
meters when the direction of observation is fixed and (ii) 
a whole field method. 

(2) Scattered Light Technique Where the Model Acts 
As a Polarizer: 


A schematic diagram of the set up is shown in 


Fig . 28* 
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The technique to be described is essentially an 
extension of the intensity calibration method for the deter- 
mination of characteristic parameters (Chapter 3)- Let the 
angle between the axes of the quarter -wave plate and analyser 
he 45°, Fig» 28. Note the intensity 1^ . Rotate the analyser 
through 90° and again note the intensity I 2 . Rotate the 
model or the incident light beam about the direction of 
observation D ; ij until I-j = I 2 • The final position of the 
incident beam corresponds to one of the primary cnaracteris 
tic directions. The reasoning is as follows. 

The intensities I-j and I 2 are the intensities of 
the components of light vector along two mutually perpendi- 
cular directions. If these are equal, this means that the 
light ellipse before the analyser has its major and minor 
axes along the axes of the quarter-wave plate so that before 
this plate the light is plane polarized. The corresponding 
■position of the light beam is one of the primary characteris- 
tic directions. To find the secondary charact exist ic dire 
ctions remove the quarter-wave plate and rotate the analyser 
for minimum intensity. In the final position, the axes of 
the analyser are along the secondary characteristic direc- 
tions. Now keep the light beam and the quarter-wave plate 
axes at 4-5° to the primary and secondary characteristic 
directions respectively. In this case the amplitudes of 
light components along the primary characteristic directions 
are equal and so the amplitudes of light components along the 
secondary characteristic directions are equal at the exit of 
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the scattered light path in the model. 'The light ellipse 
at the exit has its major and minor axes along the axes of 
the quarter -wave plate. Now rotate the analyser for minimum 
intensity and let the axes of the analyser in this position 
make an angle oc with the axes of the quarter -nave plate. 

The elliptic it y of the light ellipse before the quarter-wave 
plate is tanef . The phase angle along the secondary charac- 
teristic directions is 2X so that p = 2x where p is the 
characteristic retardation of the optical path 3^ Fig. 28. 

The characteristic parameters of the equivalent systems for 
the optical paths Mq M-j , M 3 Mq , .... etc. are thus 

found. Using this data and the governing optical equations, 
the difference and the directions of the secondary principal- 
's tresses are found as described earlier. 

In this chapter, we have discussed the assumptions 
involved in some of the previous scattered light methods. 

This is done using the analytical solution obtained for the 
case where the rate of rotation is arbitrary, but the ratio 
R = ,0*/ p ? is a constant and numerical results for optical 
phenomena for some loadings of a rectangular model under 
combined tension and torsion. For the cases shown, the 
scattered light methods developed by Frocht and Srinath and 
Swinson et . al. are found to be inadequate. Five methods of 
scattered light experimentation and analysis are described 
which neither involve any assumptions nor the rotation of 
analyser at a particular speed as in the method described 
by Robert and Guillemet. Of these, method 1(b) which 
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involves the measurement of intensities along axes V? 
apart is a very fast method of determining the characteristic 
Parameters at various points, if the computation of one pro- 
ducts and quotients of the measured intensities is made auto 
matic. As nondestructive techniques these methods hold high 
potentiality for three-dimensional stress analysis. 



chapter 5 


MISCELLANEOUS TOPICS 


Introduction: 

In the previous chapters methods have been descri- 
bed to determine the characteristic parameters for an opti- 
cally equivalent system either for a series of retarders 
or for any optical path in a photoelastic model under a 
three-dimensional state of stress. Using the concept of 
light ellipse and Stokes parameters, we shall investigate 
three additional topics in photoelasticity, namely: A. The 
Graphical Methods, B. The Theory of Compensation and C . 
Determination of Integral Fringe Orders In Photoelasticity. 
These topics are useful in the general three-dimensional 
photoelastic analysis. 

A. Graphical Methods in Photoelasticit y 

Graphical methods in photoelasticity are essentia- 
lly pictorial representations of the changes in the para- 
meters of the light ellipse as it passes through a photo- 
elastic model, or the polariscope elements. These are 
useful (i) in finding the changes in the light ellipse 
parameters as the polarized light passes through various 
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polarise ope elements, without cumbersome calculations, 

(ii) in finding the characteristic parameters of a parti- 
cular light path in a three-dimensional photoelastic model 
or for a series of linear retarders and (iii) in estimating 
the effects of various errors in the elements of the polari 
scope (exs error in the orientation of the axes of the pola 
rizer, deviation of the retardation of the quarter-wave 

plate from ir / 2 etc.)- Most of the existing graphical 

23 

methods are derived from the Poincare' sphere (see 
Chapter 3) representation of polarized forms. After a 

brief description of the graphical methods described by 

04 . 14,15,16 25 26 

Jerrard ^ , Robert and Guillemet , Cuske , Schwieger 

27 

and Cernosek , we shall describe a new and simpler appro- 
ach to graphical methods. 

Jerrard 24 has given a detailed account of Poincare 
sphere representation, its derivation and plane representa- 
tion of the polarized light by using conformal transforma- 

14,15,16 

tion. Robert and Guillemet also have discussed these. 

In the derivation of the Poincare' sphere representation 
the light vibration is first represented in a complex 
plane. The most general type of polarized light can be 
represented by a combination of two simple harmonic motions 
along the axes OX and OY as; 

x = 3 .% Cos (<5t +£ x ) = Real Ae^^ (5*^) 

with A = a x e i, " x 




FIG- 29. REPRESENTATION OF POLARIZED LIGHT FORMS IN 
A COMPLEX PLANE. 
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Y — Uy CoS ( wfc -f ) ■— F.eal Be * (5*2) 

• 4-i — i£y 

wit a r5 = a y e 1 

‘•/here as is the angular frequency, t is the time and 0^- 
£y are the phase angles. 

The ratio B/A is in general a complex quantity and 
can he written as; ' 


VJ = B/A = u + iv 


a.j ~ 1(0 x ~ ^ 
a X ~ 


where 


t an e 



= tan 0 (Cos £ + i Sin£) 

( 5 - 3 ) 


tan 6 is the ratio of the amplitudes of the components of 
•vibrations along the reference axes OX and OY, Fig. 29, and 
g is the difference between the phase angles of these 
components. So, any general type of polarized light can - 
be represented by a complex quantity 3/A = u + iv as des- 
cribed above. This is represented by point M in the com- 
plex plane u - v, Fig. 29(a). 

OM = - i Amplitude ratio 

a x 

The inclination £ of OM with the u axis is the 
relative phase difference. If this phase difference is 
increased by £ , the point M moves to- M where OM = OM* 

and MOM’ = £ . In this mapping, all points on the real 
axis represent linear vibrations since in this case, the 
phase angle £ (inclination of OM with the u axis when M is 
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on the u axis) is zero. The length of OM again represents 
the ratio of the amplitudes along the reference axes OX and 
OY. If B/A is purely imaginary, the representative point M 
lies on the imaginary axis, the phase difference £ is 7/2 


and so the major and minor axes coincide with the reference 
axes OX and OY. (Because, with respect to the major and minor 
axes the phase difference is ir/2.) Right hand vibrations 
in which the OY vibration leads the OX vibration, are repre- 
sented by points above and left handed vibrations by points 
below the real axis. Thus the amplitude ratio and the phase 
difference of a particular light ellipse referred to a par- 
ticular set of axes, can be represented in the complex plane. 
It will however be useful if the ellipticity and the azimuth 
of the light ellipse (with respect to a particular set of 
axes) are also represented in this mapping. For this consi- 
der an elliptic vibration with ellipticity b/a and azimuth cc 
(with reference to the axes-X,Y), Fig. 29(b). This can be 
defined in terms of the vibration along the axes u,v ( i.e. 
the major and minor axes are along the u,v axes), so 
u = a Cos iiit 
v = b Sin tiTt 
b/a = e = ellipticity 

Resolving along the axes OX and OY, we have, component alcng 
the OX axis = x 


(a Co s «t Cos=< - b Sina t Since) 


= Real 


(a Cos=( + i b Since) e 


icut"} 


(5 ^ ) 
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C omponenu along^ _ y ( a q os a > t gi_nc< + b Sin ^ t Cos =0 
uhe ui axis J 


= Real I (a Sin=c - i b Co 3=<) e 


i^t 


(5*5) 


sparing these equations with Eos. (5*1) - (5* 3), we obtain, 


a Sin- - i b Cos=< 
a Cos =c + i b Sin=< 


tan~< - 1 e 
1 + i e tan** 

This equation is a special form of the bilinear transformation 

w = a-L-b t (5.6) 

c + dt 

with ad f be 

(condition for conformality) 

In this conformal transformation circles and stra- 
ights lines are mapped into circles and straight lines. 

(a) Mapping of Light Ellipses of Constant Orientation 

( ** = Constant) and Varying Elliptic it y (- 00 < e < <*) 


w _ tan=< - i e 
1 + i e tan<* 

= u + 1 v 

2 2 
- tan°c (1 - e ) - i e (1 + tan op 

(1 + e 2 tan 2 *) 

so that u = C -e 2 ) t|n £. 

1 + tan^** 

u 0 as o( 0 or tr/2 (5*7) 



IV 



FIG. 30. CONSTANT ELLIPTICITY (e) AND CONSTANT 

AMPLITUDE RATIO CIRCLES IN THE CONFORMAL 
MAPPING. 


138 


2 

and v = ~ e < 1 + tan * ) 

(1 + e 2 tan 2 °0 

v -* - e or-1/e as =< -* 0 or ir/2 

Now consider two values of e i.e. e^ = e and 
= -- e, then for these values (with sane •=( ), the values 
of u are the sane --here as the values of v are equal in mag- 
nitude but opposite in sign. So, the napping is symmetrical 
about the u - axis, Also u varies with e so that it can 
take nonzero values also (if 1 ^ jt 0, ir/2). Therefore, the 

mappings of constant azimuth flight vibrations are circles 
with centers on the u - axis, (-'hen <=c = 0 or v/2, this dege- 
nerates into a straight lines, the v - axis). Also u = 0 
and v = + 1 if e = t 1 irrespective of the values of c< . 
Therefore, all the circles pass through the points P(0,1 ) 
and P' (0, - 1 ), Fig. 30. 

(b) Mapping of light Vibrations of Constant Ellipticity e 
and Varying azimuths* 

Consider two values for tan * in Eqs. ( 5 * 7 ) i.e. 
for tan * 1 = f and tan=< 2 = - f (e is same for both). For 
mappings of both these vibrations the values of v are the 
same, but the values of u are equal in magnitude but oppo- 
site in sign. Hence these mappings are symmetrical about 
the v - axis. Also v varies with ^ if e y 0 so that it 
takes nonzero values also. Therefore, these mappings are 
circles with their centers on the v - axis. When e = 0, 
this degenerates into a straight line, the u axis. 




'These two families of circles are 'orthogonal 
since OF. OP* = OH. Oil' = 1. Thus, if a vibration is represen 
ted by a point M the orientation of the ellipse is obtained 
from the circle through P, P' and M with OQ = tan*, and the 
eliipticity is obtained from the orthogonal circle through 
K with e = ON. 

(ii) Husks' s^J Circle Methods 


Consider a plane polarized light of unit amplitude 
represented by the straight line AOA as shown in Fig- 31* 

The components along the principal-stress axes p, q are Cos 0* 
and Sine’. If a phase difference oft' is added between 
these two components in the model, 

Vibration component along) 
p axis j 


Cos 0 f Cos u t 


Vibration component along 3 _ 

j) = Sin 9 Cos(uit +0 ) 
q axis ) 


Now with respect to tile axes OB and OC which are 
inclined at an angle y to the p, q axes, the vibration com- 
ponents are given by 

Component along OB = Cos^t (Cos 8 Cos y+ Sin© CosCSinY) 
- Sin©' SinuTt Sin£ Sin y 

= a f Cos ( w t + 0^ ) 

Component along OC = Cos “t (Sin 9* Cos £Cos Y 

- Cos© Sin Y ) - Sin co t Sine r Sin£ Cos y 
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and 


Sin 2 ©' Sin 2 Y + Cos 2 ©* Cos 2 v + Sin2e* Cos£ Cosy Siny 

( 5 * 8 ) 

Sinks' Cos 2 v Costs' Sin 2 y - Sin2e' CosC CosY Sin Y 


( 5 . 9 ) 


o find tie maximum and minimum values of 


and b 1 (corres- 


ponding to the lengths of semi -major and. minor axes a and 
b), we differentiate the above expressions with respect to 
y and equate to zero 


,2 ,2 
da db 
dY " dY ~ 


giving tan2y = Cos£tan2§’ (5*10) 

Substituting this expression in terms of y in Sqs. ( 5 * 8 ) 
and (5*9)? we obtain 


a 2 = -k j 1 + (Cos^e 1 + Cos 2 £ Sin 2 2 ©* )*J (5*11 ) 

= 2" (1 + J ) 

where J = (Cos 2 2e’ + Cos 2 £ Sin 2 2e' )£ 

and b 2 = i (1 - J) (5*12) 

The quantity (Cos 2 29 r + Sin 2 e* Cos 2 e is called 
J. The elliptic ity and the orientation of the major axis 
can be represented by an arrow (vector) having length equal 
to J and oriented at a chosen angle. The angle between 
this direction and a chosen reference direction are doubled 
and so are all the other angles involved. 

t 

Wow, <1 = 2a^‘ -1=1- 2b 2 = 

1 + e 1 

• where e = § = elliptic ity 

SI 


( 5 * 13 ) 
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So, the magnitude of J depends on the ellipticity 
only. Further, 


e 2 = L- . r - 4 (5-1^) 

1 + J 

So, J = 1 represents plane polarized light and 
J = 0 represents circularly polarized light. That is in 
Fig. 32, a tip of the J "vector on the circumference of the 
circle means plane polarized light, a tip in the center 
corresponds to circularly polarized light while all other 
positions represent elliptically polarized light. Since 
all the angles are doubled, the mutually perpendicular 
principal-stress axes are represented by the straight line 
qOp, since qOp =v. If plane polarized light is at an angle 
6* to one of the principal-stress axes, the components along 
the principal-stress axes are Gos 0* and Sin 6* . Consider 
the plane polarized light represented by the point A on the 
circumference of the circle. The line AC is perpendicular 
to line -qCp. 

Now, 

q B = 1 - Cos2 0 * = 2Sin 2 ©* 


and BP = 2Cos 2 e* 

i.e. q B and BP are proportional to the squares of the amp- 
litudes of components along the principal-stress axes. In 
order to find the intensities of light components along two 
mutually perpendicular directions', a perpendicular is to 
be drawn from the tip of the J vector on the straight line 
which corresponds to these mutually perpendicular directions* 
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The intercepts correspond to the required intensities. Since, 
during the addition of retardation along the principal-stress 
axes, the intensities of components of light along these 
axes remain constant, the representative point moves along 
the straight line AC which is perpendicular to the line qOp. 

At the point A, 

J = (Cos 2 2 9 5 + Sin 2 2©' Cos 2 G)^ = 1 

or Cos 2 ^ = 1 

i.e. G = 0, it, , nir 
At the point B, 

J = OB 

= - Cos2©' 

= (Cos 2 2 ©* + Cos 2 C Sin 2 2©* )i 
or Cos 2 G = 0 

, (2n 4- 1 k 

i.e. G = tt/2, 3V2, 2 

This can be interpreted as follows; At the point 
A, the phase difference is zero. As the phase difference 
is increased toir/2, the representative point moves to B 
along line AB. As the phase difference is further increased 
t s ir it moves to the point C along line BC . When it is 
3ir/2, the representative point comes back to B and so on. 
(iii) Schwieger 's 26 Analysis; 

Schwieger presented some rotation rules to be used 
with J circle method. He introduced the Wulff’s grid method 
of crystallography into photoelasticity and described methods 
to find the characteristic parameters for a combination of 




parallel projection of the sphere of 
polarization (J CIRCLE ) 
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linear retarders- A "brief suinary of his analysis is pre- 
sented below: 

The four parameters of light ellipse, namely, , 
the azimuth, the angle 0 (or tan 0, tne amplitude ratio), 
the angle w (or tan a> , the eilipticity) and 0 the phase 
difference. From the equations of Chapter 2 we have the 
following relations connecting these four parameters. 

SinC Sin2 © = Sin2e (5^5) 

C os £ tan2 © = tan2c< (5-16) 

Cos2 c < Cos2 a = Cos2 6 (5*^7) 

The angles of a spherical triangle satisfy these 
relations and so these parameters can be represented on 
Poincare’ sphere. 

(a) Parallel Projection: 

In Fig. 33. all points of the upper hemi-sphere 
represent left elliptically polarized light while points in 
the lower hemi-sphere represent right elliptically polarized 
light. The area of the J circle can be obtained by parallel 
projection on a plane having a normal which is parallel to 
the axis connecting the poles of the sphere. 

The chords drawn in the J circle, Fig. 3*+ * are 
lines of constant amplitude ratio, the ellipses are loci of 
points of constant phase difference. The radial lines are 
lines of constant azimuth. The actual azimuth will be half 
of the angular inclination of these radial lines with the 
OK axis. The reference axes are the lines OX and OY. The 
scale. on the right hand side of this plane representation 
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RG.3 5- STEREOGRAPHIC 
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is calibrated in terns of elliptic ity. This plane represen- 
tation can be readily used to find the resulting photoelastic 
effect after the passage through a set of linear retarders. 
Addition of phase difference corresponds to moving along the 
chords of constant amplitude ratio, the length of movement 
corresponding to the value of retardation added. Rotation 
of reference axes corresponds to moving along a constant 
elliptic. ity circle until the angular inclination to OX is 
increased by double the amount of rotation of reference axes, 
(iv) b.-.Iff 's Grid- 

This is the stereographic projection of Poincare 1 
sphere of polarization. 

Tho stereographic map is the projection of a sphere 

on tho plane of a great circle, Fig. 35 , in such a way that 

tho projection of a point A on tho sphere is the point in 

which che line through A and a pole of that great circle 

intersects the piano of the groat circle, a hi s projection 
44 

is conformal 

The chords of constant amplitude ratio in the para- 
llel projection correspond to the curved lines surrounding 
the YOX axis in the stereographic projection. The radial 
lines here again represent lines of constant azimuth. The 
loci of points of constant phase difference are arcs of 
circles in this projection. 
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Determination of Parameters of the Equivalent Optical 
System For A Combination of Two Linear Retarders-: 

A = Retardation introduced in the first, plate 

= Retardation introduced in the second plate 
a 

The optically equivalent system consists of a linear 
retarder and a pure rotator. If the incident light is plane 
polarized along the axes of this linear retarder it remains 
plane polarised after the retarder but is rotated by the 
pure rotator. That is, it is plane polarized at the exit. 

To find the rotatory power of this pure rotator and the orie- 
ntation of this linear retarder draw two lines CK^ and CK^ 
corresponding to the axes of the two plane models. Cons- 
truct the circular arcs on these lines corresponding to 

- constant and 180 - A ^ = constant. From this point of 
intersection P of these two arcs, proceed to the point Q and 
R on the circumference of the circle along the constant 
amplitude ratio curves corresponding to the axes CSC-] and 
0X 2 « The angle ROQ will be twice the rotatory power of the 
system. Also, OQ = 2P re pp, where P re pp is l* 1 * 3 angle bet- 
ween the axes of the first plate and the axes of the linear 
retarder in the optically equivalent system. The procedure 
can be easily understood from Fig. 3 7* In this figure 
?i P 2'^2 re P resGn ts the passage of light through the two 
retarders.. 

P^P 2 = 2G 

The actual problem is to locate the points P^ and 




FIG. 36- DETERMINATION OF THE ANGLE <j> ACTIVE 

OF THE ACTIVE PLATE OF AN EQUIVALENT 
SYSTEM FOR TWO PLANE MODELS. 




FIG- 38. REPRESENTATION OF AXES OF 
THE TWO RETARDERS. 
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P_ which lie on the constant phase difference curves A = 

and A = 180 - A^ respectively such that 
A_ 

P 1 0P 2 = 20 

This is done by rotating one curve with respect to 
the other by an angle 20 as in Fig. 36. so that their 
point of intersection P corresponds to the points P^ and ? 2 
of Fig. 37. Fig. 3 8 . shows the orientation of the axes of 
the recorders. 

In Fig . 36 • , 

r$q = 2 (e 2 - e 1 ) 

A 

ROQ f 2 p ac+ . ive as mentioned £y Schwieger. 

To find the retardance A re pi of the optically 
equivalent system consider plane polarized light along the 
axes of the first plate represented by the point P in Fig. 
39. Retardation A in this plate brings it to P^ . Rotation 
of the reference axes to make them coincide with the axes 
of the second plate brings it to the point P 2 * Addition of 
a retardation A 2 in the second plate is represented by the 
path P 2 P 3 so that the point P^ represents the state of pola- 
rization at the exit of the two plates when referred to the 
axes of the second plate. We should arrive at the same 
point P_ even if we use the 'optically equivalent system. 

3 ' A 

Now proceed back to the point where P^Oi^ = 2p ac tive 
and Oil = 0P 3 . Find the point P^ on. the circumference such 
that PGP^ = 2©^ . That is, the point, P^ represents the state 
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of polarization at the entrance -hen referred to the axes 

of the equivalent retarder. The points and P, lie on 

5 ©- 

the same curve of constant amplitude ratio. (This fact can 


be used to check the values of ©^ and c ac ^j-,-o* > The Path 
P^P^ represents the equivalent retardation A rnp ^. In the 
case of systems with many number of retarders, this method 
can be used successively. Schwieger, in his example with 
ten plates, mentions an iterative procedure where the azimuth 
of the incident plane polarized light is varied until the 
light at the exit is plane polarized. 

(vi) Graphical Method Using Quaternions 

* 27 

Cernosek devised a geometrical method based on 
the properties of quaternions to find the characteristic 
parameters of a system of linear retarders. He made use of 
Jones vector in this derivation. The principle of this 
method is described below. 

Let the Jones vector for a certain elliptic ally 
polarized light referred to the X,Y axes be given by 


i(£y + to t) 


a e 

Y 


l(Gy +u t) 


where £y £y are the phase angles and A^ Ay are 


the amplitudes of the components along the X, Y .axes 


respect iyely . r 

,1 0 

Let I = 

10 1 


i 0 
0 -i 


; o = 


o 1 

-1 0 


k = 


0 i 

1 o 


( 5 . 18 ) 




FIG-40- CO-ORDINATE SET USED FOR 

THE DESCRIPTION OF THE GENERAL 
QUATERNION Mi,M 2 ,M 3 ,M 4 ARE 
THE COEFFICIENTS OF THE QUATERNION. 
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The general unitary 2x2 matrix, which represents 
the transf or mat ion of polarized light from its entrance th- 
rough a system of linear retarders to its emergence can be 
uniquely expressed as a quaternion 

M = Mj I + MU, i + M. j + k ( 5 * 19 ) 

The coefficients of the quaternion M are connected 
by the unitary norm relation 

hf + + M3 + '1^=1 (5*20) 

The quaternion can be described graphically as 
shown in Fig. 4-0. The co-ordinates of the four points 
K-j I <2 ana Kq. in Fig. 40. are the coefficients of the 
quaternion M defined by Eq. (5*19), ^ 

and M 2 % respectively in four different co-ordinate systems 
(m m ), (m m ), (m m ) and (m m ) which are in the same 
plane but have different orientations. The values of M-j ^ M 2 ^ 
M 3 and ^ for a rotator and a retarder placed before or a 
after the system can be easily found from the transformation 
relations. Bach of these transformations corresponds to a 
particular type of rotation of some of the four radius 
vectors OK^ 0 K Q OK^ and OK^. Thus, the final positions 
of the points K-j K 3 and can be found graphically very 
easily for either a retarder or a rotator when it is placed 
either before the system or after the system by simple rota- 
tion of the radius vectors and or and K^. Two 
other points in each step will then be found by drawing 
lines parallel to the coordinate systems. 
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The Characteristic Quantities of a General System of 
ilet ardors? 

Before studying a system of retarders, ve shall 
presort the quaternions for simple operations. 

Quat err ion of rot at i on : 

S 0*.,) = Cos* • I - Sin* 2 . J 

=*2 is the rotation angle. 

Quaternion TW = S (* ) M when S ) is placed after M. 

2 a ^ 

then 

U — A-j 2 ^- A22^ A^2 j -^. 2 ^ (5*21 ) 

where 


A. 

= M. 

CosF< 

+ 1L 

Sir?* 

(5*22) 

12 

1 

2 

3 

2 


OJ 

= m 2 

Cos* 2 

- \ 

Sinp<2 

(5*23) 

A 

= M 

Cos'* 

- M 

Sin<* 

(5*24) 

32 

3 

2 

1 

2 


\ 2 


Cos °( 2 

+ M 2 

Sin°( 2 

(5.25) 


Quaternion for a linear retarder of retardance A 2 is given 
by 

G(A 2 )=Cos a 2 * i+ Sin & 2 

If this retarder is placed after M, the quaternion 
for the combined system is given by 

U ( A 2 ) = G ( A 2 ) M — 2 ^ + (5*26) 

B = M Cos A 
12 1 2 


where 


- M Sin A 
2 2 


(5*27 ) 
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— x CosA^ i a Lj oinA 2 


B^p = CosA 2 - ^ SinA 2 
CosA 2 + Mg Sin A 2 


( 5 * 28 ) 

(5-29) 

(5-30) 


If the rotator and. the retarder are placed before the system 
M, the corresponding quaternions for the combined systems 
will be as given below? 


U (°f ^ ) — ^ 0 (°^q ) = A-j -j I + A 2 "ji + A^-j 1 + -'^L}_q 1 c 
whore 


and 

U* ( = MG(4 1 )=B 11 I+ B| 1 i + B 31 j + B^k 

where 


(5-31 ) 


A 11 


o 

o 

+ M 3 

Sincr 

1 

(5-32) 

A 21 

= M 2 

Cos^ 

1 

+ 

Sine* 

1 

(5-33) 

A 31 

= M 3 

Cosoc^ 

- M, 

SinoOj 

(5.3k) 


= \ 

Cos^ 

1 

- M 2 

Sirf* 

1 

(5.35) 


(5-36) 


B n 

M 1 

C os A 

1 

-m 2 

SinA„ 

1 

(5*37) 

B 21 

II 

C OS A 

1 

+ M 

1 

Sin A 

1 

(5-38) 

B 

31 

ii 

COS A 

1 

+ X 

Sin A 

1 

(5*39) 

\l 

= % 

Cos A 1 

-m 3 

SinA^ 

(5*ko) 


The coefficients of the resulting quaternion U c 
after a rotation through an angle =< 2 can be foimd by 



m] 



FIG .41. TRANSFORMATION OF A CO-ORDINATE 
SET AT THE POINT OF EMERGENCE «( 2 
AZIMUTH OF THE OLD AND NEW 
REFERENCE CO-ORDINATE SYSTEM. 


FIG. 42 THE ADDITION OF THE RETARDER AT 
THE POINT OF EMERGENCE .COORDINATE 
SETS m 1 -m 2 AND m 3 -rri 4 ROTATED IN 
OPPOSITE DIRECTIONS, 2A?:PHASE RE- 
TARDATION OF THE ADDED RETARDER ' 
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FIG. 43. 


THE CHARACTERISTIC QUANTITIES OF A 


GENERAL SYSTEM OF RETARDERS, ARE 
THE ORIENTATIONS OF THE PRIMARY AND 
SECONDARY CHARACTERISTIC DIRECTIONS 


2p IS THE CHARACTERISTIC PHASE 
RETARDAim , . v . . 

•' tan, 

4 , , 

' M/ !, * v ' ' 

tan ftsl- , . 
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rotating the radius vectors of the points and in the 
clockwise direction, Fig. 4-1. as can be soon from Eqs.(5*22)- 
(5*25). In the same way the coefficients of the quaternion 
for a retarder whose axes are inclined at an angle oc, with 
the reference axes can be found by rotating the radius vec- 
tors of and (multiplication by the quaternion of 
rotation) through an angle ^ in the opposite directions 
and then rotating the radius vectors of and by an 
angle A-j , both in the same direction (for the addition of 
a retaraance of A-| in the linear retarder), Fig. 4-2* 

Similarly, the geometrical procedures correspond- 
ing to the equations ( 5 * 32 ) - ( 5 * 35 ) and ( 5 » 37 ) - ( 5 » ! + 0 ) 
can be established. 

Using .the above relations, ve have the following 
representation of the quaternion for the combination of a 
linear retarder and a pure rotator (of the optically equiva- 
lent system), Fig. 3 . Following the procedures described 
above, the quaternion for a system of linear retarders can 
be found graphically. Comparing this with the quaternion 
for the optically equivalent system. Fig. b-3* the characte- 
ristic parameters of the system of linear retarders can be 
found. In Fig. 4-3. the characteristic quantities of a gene- 
ral system of retarders are given by 

e x = °c + f? 

^2 = °f-P 
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cCj £ ar ,3 the orientations of the primary and secondary cha- 
racteristic directions. 2p is the characteristic phase 
retardation . 


ban 0^ 

= tan 

(«f + P) 

Mo 

= 

tan £ 0 

= tan 

(* - B) 

= % 

2 



Mo 






v 2 

m2 p- 


tan p 

x 2 

- ^ ) 2 



?i 2 

+ 



'1 

3 



The New Approach 

The light ellipse can be characterised by any two 
of its four parameters which are related by the following 
equations (from Chapter 2) 

S = ai -a 2 = Cos29 = Cos2w Cos2c< (.5^ ) 

1 A I 

S = 2a v a v Cos£ = Sin2e CosC = Cos2w Sin2°<(5 ) 

cL Js. X 

= 2a x ay Sin £ = Sin2© Sin€ = Sin2 u 0^3) 


Also 





ay,a x are the amplitudes of the components of light along 

the X, Y axes. £ is the relative phase difference between 

these components. is the azimuth and tan wis the ellip- 

ticity of the light ellipse tan© = — * 

a x 

So, we have actually two independent equations 
connecting the four parameters (the azimuth, amplitude 
ratio, ellipticity and phase difference) of the light ellipse* 
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So, the variation of any two of these parameters with respect 
to one another can be plotted keeping one of the other two 
parameters constant. For example, the azimuth o( and the 
phase difference 0 can be plotted along X and Y axes so that 
we have constant amplitude ratio curves and constant ellip- 
tic ity curves. We will have other systems of graphs if 
other combinations of the parameters are taken with different 


coordinate systems. 

From the above equations, we have 


and 


tan€ 


_ tan2 M 

“ Sin2* ” (1 _ e 


2 av Uy 

tan 2 c < = — 2 2 CosC 

a X " a Y 


2e 

2 ) Sin2 c * 


, 2e X 

0 


cose 


where 

e = b/a = ellipticity = tan w 
e X = = am Pli tude ratio = tqne 


(5.M+) 

(5.M) 


From these relations, it is clear that for given 
values of e and °f , the values of C is fixed and for given 
values of e x and 0 , the value of * is fixed. The values of 
these parameters can be tabulated and this table can be used 
for computing the effect of a series of retarders.' 

From Eqs. (5. Mi-) and (5*^5), we note that for a 
given value of e the ellipticity (not equal to unity), the 
phase difference 6 has a minimum value which occurs at 
of = 1+5° . For a given value of e-^ the amplitude ratio (not 
equal to unity), the azimuth °< has a maximum (or min imum) 




FIG- 44- THE FIRST GRAPHICAL METHOD 
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value which occurs at £ =0. When the amplit-uio ratio 
equals unity and the phase difference £ f 90° ? the azimuth 
°C equals ^+5°» When e^- = 1 and£ = 9^°? it is circularly 
polarized light and °f takes all values. 

One can plot Eqs. (p.4V) and (5 *~5 ) 3-s shown in 
Fig. 44.. jn this figure, the horizontal axis represents 
the azimuth °c.and the vertical axes the phase differences . 
Use of Eq. (5**+^+) gives ’’Constant Elliptic ity” curves. To 
get these, one plots vs £ 5 for fixed values of e getting 
curves like APB, AEB, etc. Use of Eq. (5.^5 ) gives "Constant 
Amplitude Ratio” curves. These are obtained by plotting 
c( vs £ for fixed value of ey . At points A and B the phase 
difference £ is f/2 and the azimuths °fare 0 and v/2 respec- 
tively. So points A and B represent states of polarization 
when the reference axes coincide with the major and minor 
axes of the light ellipse. The line AB represents circu- 
larly polarized light and the lines AO, 0C and CB represent 
plane polarized light. This form of graphical representa- 
tion of light ellipse is very convenient, since in photo- 
elastic analysis, we add a certain amount of retardation 
along a particular set of axes (during this addition, the 
amplitude ratio along these axes remains constant and so 
this process is represented by movement along the constant 
amplitude ratio curves) and rotation of the reference axes 
(during this process, the ellipticity remains constant and 
so it is represented by movement along the constant exlip- 
ticity curves). The same graph can be used for ranges 




OR- Axis of the first retarder 
OF^-Axis of the second retarder 
OP -Reference axis 


FIG- 45- ORIENTATION OF THE AXES OF 
THE TWO RETARDERS. 
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beyond ir /2 radians for the azimuth and the phase difference 
by taking into account the sign variations and using Eqs. 

(5 .44) and (5 *45). This is demonstrated at the end of this 
section. 

Procedure to Use the Gr a phical Method: 

Consider two retarders, R-j and whose axes are 
inclined with respect to each other by an angle © 2 and - let 
the axes of *make angle 9^ and 6 - | + 7/2 with respect to a 
given fixed axis, Fig. 45. Let the initial light be plane 
polarized along the refering axis. This initial light elli- 
pse is represented by point 0 (with £ = 0 and 1 * =0), Fig. 

46. When this is referred to the axes of iL , the initial 
light ellipse will be given by point (with e = 0 and amp- 
litude ratio e-£ given by point ). Addition of retardation 
'--I by R^ is represented by movement along constant amplitude 
ratio curve P-j P 2 up to the point P^, such that the ordinate 
(phase difference £ ) of P 2 is greater than that of the point 
P^ by an amount e^* Now, we refer this light ellipse to the 
axes of the second retarder by moving to the point P^ from 
the point P^ along the constant ellipticity curve such that 
the azimuth of the point P^ is greater than that of the 
point P by an amount © P , where © is the angle between 
the axes of the two retarders. Addition of a phase differe- 
nce of £ 2 in the second retarder is represented by movement 
along the constant amplitude ratio curve P^ from the point 



FIG- 47. VARIATION OF LIGHT ELLIPSE PARAMETERS 
IN A PLANE POLARISCOPE . 


i 

i 
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the Point P^ where the ordinate £ of the point 3^ 
is greater than that of the point by an amount £ 2 ? where 
is the retardation introduced by the second retarder. 
Finally, the point P^ in the chart represents the light 
vibration at the exit of the second retarder with respect 


to its axes. So, during the representation of the varia- 
tion of light ellipse as it passes through a system of li- 
near retarders on this graph, the point under consideration 
always represents light vibration with respect to the axes 
along which the retardation has just been added or about to 
be added. 


Passage of Light through Elements of a Polariscope: 


(i) Plane Polariscope: 

In Fig. 47. point 0 represents the plane polarized 
light along the axis of the polarizer; point P^ , the plane 
polarized light referred to the principal-stress axes of the 
model at its entrance; point "P 3 the light at- its exit after 
an addition of 0-j retardation referred to the axes of the 
model, and P , the same light vibration but referred to the 
axes of the analyser which is kept crossed with respect to 
the polarizer. If one of the axes of the retarder (i.e. 
model) coincides with the direction of the initially pola- 
rized light, then we move from 0 to P 2 *, where OP^ = ♦ 

When the analyzer is at ir/2 to polarizer axis, we move along 
the constant ellipt.icity line P 2 OEC (e = 0) until we come 
to the orientation of the analyzer axes, that is, point G, 
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The amplitude ratio is zero and no light emerges tHe 

analyzer . 

(ii) Circular Polariscope: 

Fig . 4-8 shows the elements of the circular pola- 

riscope. The taints ?., and P 0 in Fig-. 4-9. represent light 
vibrations before .rad after the first quarter -wave plate 
referred to its .axes, she points P^_ , P^ •> the light vibra- 
tions before and after the model referred to its principal- 
-stress axes at the particular point under consideration. 

The points 1^, represent the light vibrations before and 
after the second quarter -wave plate with respect to its axes 

and the point P , the light vibration before the analyser 
o 

referred to its axes. Now, we have to prove that Gy^ , the 

amplitude ratio corresponding to the point P depends oily 

6 

on the retardation introduced in the model and not on the 
angle 0 , the orientation of the principal-stress axes with 
respect to the axes of the quarter-wave plate. In a graph 
of Fig. 49,. which is drawn to scale, point will always 
lie on the e., = e v line, irrespective of the value of 0 

JL 

used. This can also be shown by Eqs. (5.44) and (5* 5 + 5) 
that are used to construct Fig. 49. To show this we proceed 
as follows: 

Fran Eqs. (5*4-4) and (5*4-5)? we obtain 

tanC = tan3i ( 5 . 46) 

Sin2c( 

and tan2°< = tan2 8Cos£ (5*4-7) 
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where tan© is the amplitude ratio. 

Let £ 3,0 k, £ 5 , £6 he the phase differences at the 
points respectively. 

From Fig . 4-9 

e 3 = ir /2 - e* 


_ TT/ Lj. 
3 - /H- 


tan-, -3 = -fcan2 . 6 i3 = tan 2 w 
J Sin2 c 3 


2ui 3 = = tt/2 - £ 


For the point P, 


4 


tan^ = 


t an 2 o. 


3 


TT 


Sin 2(__ + 0) 

4 


= Cotf 
Cos20 


e = 0 - tt /2 

5 4 

. . -can £> = - Cott j, = - * 

5 v Cote* 

from Eq. (5*5^ ) 
Also, for points 3+ and- P 5 


tan 2 (j5 + 0 ) = - Cot 20 = tan 204 . CosC^ 
and tan 2 = tan 29^ Cos £^ = tan 20^ Sin£ 4 

. . tan 2°^ = - Cot 20 tan £4 

_ Cot 20 Cote* 

Cos 20 

from Eq. (5.49) 

_ _ Cot £* 

Sin 20 


( 5 . 4 - 8 ) 


(5-4-9) 


(5-50) 


(5-52) 

(5*53) 


(5-54) 
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ITow, for the points and we have 


and 


tan 2<v = tan 6 - Sin 2=c 

5 5 5 

tan 2^ = tan 


(5-55) 


= tang 6 Sin 2*^ 


tanG.^ Cos 2*^- 


(5.56) 


; an £ 

o 


tan Gc- tarr<2 ^ 

5 5 


°C _ of _ 7T Al 

6 - 5 /H ^ 


= - Cot 20 


from Eqs. (5 *5©) and (5*5^) 


= ir/2 + 20 


(5-57) 


[■low for point P 


tan 20/ 


tan 2< , Cotp<2 


5 


* Cos G 6 Sin 20 

from Eq. (5 *57) 

Sin 20 


Cot £* Sin 20 

from Eq. (5 .5*+ ) 
- tan ( it- G*) 


or 


2©^ = ( it - G ) 

„* 


G 6 


jr _? 
2 " 2 


(5-58) 
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So, the amplitude ratio tan© g = 
retardation £ introduced in the model and not on the orien- 
tation 0 of the principal-stress axes with respect to the 
axes of the first quarter -wave plate. In the case of light 
field (the quarter-wave plates are crossed but the polarizer 
and analyser are parallel) polarise ope, 


depends only on the 


©v = ror 

whenever £* 1 = (2 n + 1 ) ir 
then 


tan © 


6 


a Y/ a X 


0 


i.e. the extinction of light occurs. 

In the dark field set up, the analyser and polari- 
zer are crossed and the two quarter -wave plates are crossed. 
In this case 

*6 = *5 +7r/4 
Here Eq. (5*56) becomes 

tan 2^ = tan Cos 

and tan G^ = Cot 20 

= tan (1 - 20) 

2 

or e 6 = ir/2 - 20 (5.57) 

becomes 

© 6 = £*/2 




FIG. 50- OPTICAL SET-UP IN TARDY'S METHOD. 




wKKm 



wKKm 


182 


so tli at 

t an , = tan£*/2 
6 

rz 0 


i.e. extinction of light occurs. 

•fieri e n = 2 a>r 

These results cm bo proved in exactly the same way 
for other alternative set ups of the dark and light field 
arrangements. 


(iii) Tardy’s Method: 

The set up is shown in Fig. 50. Point 0 in Fig. 51 


represents plane polarized light referred to the polarizer 
axis. Points 1 and 2 represent light before and after the 
first quarter -wave plate referred to its axes. Points 3 
and 4 (path 3-2-4) represent light before and after the 
model with respect to the principal-stress axes. Points 5 


and 6 represent light before and after the second quarter- 
-wave plate. In order to cut off light, the amplitude ratio 
along the analyser axis must be 00 . This can be achieved 
by making point 6 come to point 7 along the constant elli- 
pticity line. The orientation of the analyser axis is then 
- c<^) W hose value can be found as follows: 

Fran Fig. 51 we note that 

^ = it /2 - e* 




tan 2 


\ 


tan€ ^ Sin 3=< 4 
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= tan £ , 

M- 

= tan (!L - £*) 

or = ir/2 - £ 

Als Oj 

£5 = tt /2 

From 3q. (5-^3) 

Sin 2© SinG = Sin 2o 
when £ = ir/2 

© = to 

. . ©^ = (o^ = © 6 

tan 2°(g = tan 2©^ Cos = tan 3^ 

* 

i.e. c( 6=( o 5 =^=|-| 

The angle between the axes of the second quart c. 
-wave plate and analyser is given by 

*7 -«6 = -"6 

. . The orientation of analyser axes is given by 
* = (c < 7 - c< 6 ) + ir/f 

= 3 ir A - * 6 



or £*= 2< - * 

From this graphical representation the essential 
features of light ellipse can be observed clearly. For any 
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light ellipse, for' directions nearly at b^° to its major or 
minor axis, the variation of phase difference with direction 
is very small. (Because, the constant ellipticity curves 
are almost flat for values of =C around *+5° ) . This variation 
is high for directions close to the major or minor axis. 

Also, we see that as the ellipticity nears unity, the varia- 
tion of amplitude with direction becomes small. This chart 
will be very useful in finding the effects of small errors 
in the retar dances of quarter-wave plates readily. 

Determination of Characteristic Parameters: 

To find the characteristic parameters of a system 
of linear retarders Schwiegar suggested a geometrical cons- 
truction and an iterative method, and cernosek presented 
a method based on quaternions. A much simpler and direct 
method to find these characteristic parameters will now be 
described. If the incident light on the equivalent optical 
system is circularly polarized, the amplitudes of the compo- 
nents of light vibration before and after the linear retar- 
der along its axes (primary characteristic directions) are 
equal, and these components are just rotated in space by the 
pure rotator to coincide with the secondary characteristic 
directions (Because, the angle between the primary and sec- 
ondary characteristic directions is equal to the rotatory 
power). So, the mutually perpendicular directions at the 
exit of the system along which the light, components have 
equal amplitudes are the secondary characteristic directions. 
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This method can be used with an y of the above described 
graphical methods (except the one using quaternions) and 
with the second graphical method to be described now. Once 
the secondary characteristic directions are found, the pri- 
mary characteristic directions and the characteristic retar 
dation can be found directly by using a different type of 
incident polarized light other than circularly polarized 
light. An example of finding the characteristic parameters 
for a combination of two retarders using this chart will be 
presented at the end of this section. 

The Second Graphical Method: 

Squaring .and adding Eqs . (5.4t)and (5.4-2), 

we obtain 

Sin 2 2© Cos 2 e + Cos 2 2© = Cos 2 2« 

Cos 2 2©Sin 2 e + Cos 2 C = Cos 2 2 oj 
C os 2 2co = x and Cos 2 2© = y 
y = - Cot 2 C (5-59) 

Sin 2 0 

For a given value ofC , this equation represents 
a straight line. 

Since the set of values x = 1 , y = 1 satisfies 
Eq. ( 5 * 59 ) irrespective of the value of £ , Eq. (5*59) repre 
sents a family of straight lines all passing through the 
point x =1, y = 1 . Also, 

if 0= 0, x = 1 ; and if € =ir/ 2, then x = y 


or 

Let 

then 
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if j = 0 

x = Cos 2 e = Cos 2 2 05 
or e= 2 05 

Squaring and adding Eqs. (59+2) and (59+3), we obtain 
Sin 2 26 = Cos 2 2 o) Sin 2 2=f + Sin 2 2e 
or j = x C &s 2 2c< (5 .60 ) 

If x = 1, y = Cos 2 2=<= Cos 2 2© 

or c< = e 

The families of straight lines represented by Eqs. 
(5»5 9) and (5.60) can be plotted as shown in Fig. 52. 

Eq. (5*59) represents straight lines of constant 
phase difference AC, AG, AF, AB,.... etc. Eq. (5.60) re- 
presents straight lines of constant azimuth BA, BE, BD, BC, 
etc. The ordinates are lines of constant ellipticity 
and the abscissae are lines of constant amplitude ratio. 

Geometrical Construction of x = Cos 2 oj And 

y = Cos 2 2 © Scales (Figs. 53 And 5*+ ) 

A n 

In' Fig. 53. BC- = AC = 1 and BCA = 90°. Draw a 
semi circle on BC with 0, the mid-point of BC as the cent er. 
To construct the constant phase dif ference line ( £ = £ ^ )*► 
take point D on the semicircle such that d$D = 2£<| . Draw 
the perpendicular DD^ to BC . D|A is the required constant 
phase difference line with € = £ 1 . The proof for this is 
as follows: 
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In the graph point B is the origin and line BC is 
along the x-axis. Also x = Cos 2 2m and y = Ccs 2 0 . 

On x-axis, 

0 = 2a 

as sho’-.n above 

Nov in Fig. 53, 

BD 1 = BO + 0D 1 

= i + i Cos 20 1 

= Cos 2 ^. 

If the ellipticity at the point is tan o>, the 
phase difference € is given by the relation 

€— 2oj 

Also, the value of x at is given by 
X = Cos 2 2a 
= Cos 2 C 
or2to = 0^ 



Vie know that all constant phase difference lines 
pass through, the point A (1,1). Therefore the line AHj 
is the constant phase difference line \*rith 0=0^. 

The procedure for construction of constant azimuth 
lines is exactly similar. In Fig. ♦ ABC is the circle 
on AC with Op as center (AO-j = O^C ). To construct the 

constant azimuth line << =<< -j , locate the point E on this 

A 

circle such that A (Xj E = 4°^ . Draw the perpendicular EEj 
on AC from the point E to meet AC in ij. 
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For point E. , 


x = 1 


Cos^ 2 9= C 0 1 + O-jE 
Cos 


Cos 2 2*, 


. . 0 = °f ^ 

Also = 9 when x = 1 (shown above) 

. . c( = o( ^ 

So, the azimuth at the point isc*^. 

We know that all the constant azimuth lines pass 
through the origin B. So is the constant azimuth line 

with ** = ** ^ . 

To increase the phase difference from C-j "to e 2 
we move along the constant amplitude ratio (i.e. parallel 
to X-axis) and for rotation of the reference axes from azi- 
muth *** to °f 2 we move along the- constant ellipticity line 
(i.e. parallel to Y-axis). 

The same chart can be used for other ranges also 
i.e. beyond v / 2 but the proper signs are to be taken into 
account. This is demonstrated at the end of this section. 

The striking advantages of this method are that 
(i) it is very simple and easy to draw since it involves 
only straight lines and (ii) the nonlinear scales can be 
easily constructed by using the method described above. 
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Tbs analysis for plane and circular polarise ope us- 
ing this chart is exactly similar as that described for the 
first graphical method. 

Tardy T s method : 

In Fig. 55* the notations are same as those used 
in Fig. 59. Point 0 is plane polarised light referred to 
the axes of polarizer. Point 1 represents light before the 
first quart or -wave plate referred to its own axes. Addition 
of a retardation of ir/2 is represented by movement along the 
constant amplitude ratio line (parallel to X-axis) from G = 0 
line (A C ) to G = ir/2 line (A B), that is to the point 2, 
which represents circularly polarized light. The same vib- 
ration, but referred to the principal-stress (inclined at 
45° to the first quarter-wave plate axes) is represented by 
the point 3 (which coincides with point 2). This rotation 
of the reference axes through 45° is represented on the 
chart by movement along the constant ellipticity line (i.e. 
Parallel to Y-axis) from = 45° line (B C ) to°c = 90° line 
(A B ). Addition of a retardation of €* in the model is 
represented by movement from point 3 to point 4 along cons- 
tant amplitude ratio line till the representative point is 
on the constant relative phase line with G = T/£ir+ G . 

Light vibration after the model but referred to the axes of 
the second quarter -wave plate is obtained by moving the, re- 
presentative point from 4 to 5 along the constant ellipticity 
line. Subtraction of a retardation of ir /2 in the quart er- 
-wave plate is represented by movement along the constant 
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amplitude ratio line from point 5 to point 6. In order to 
cut off light, the amplitude ratio along the analyser axis 
must be 00 . This can be achieved by making point 6 come to 
point 7 along the constant ellipticity line. The orienta- 
tion of the analyser axis is then -of^) whose value 
can be found as follows: 

From Fig. 55 . 

x c - *4. = x 6 - x^ = y ? - y 6 
But X Q = 7 r, = 1 

i *e . Co 3 ^ 2w^ n C os2 2©^ oi* toLj. = ©^ 

Also Gq - 2w r - 2<^ 

( • • when j - 0,0 = w ) 

And C c = =0 


Also 


But 


£!+ - 2^ 

*7 -*6 = e 7 " e 6 

( 

°c 7 = e 7 = 0 

Of _ 0 

6 6 

” 4 - = 2 ^ = 20 £= 2=( 6 


when x = 1 , e =oc ) 


If the retardation added in the model is P then p = 
£4 - ir/ 2 . if the angle between the analyser axis and one 
of the principal-stress axes is p , then 




09P 



FIG- 57. FIRST GRAPHICAL METHOD ( o<<=« 2TT, o :<■ € < 2T). 
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:< c = v/k- - P 

e = G 4 " 7r / 2 
= 2*6 - u /2 

= - 2 P 

Derivation of J Circle Method Without Referring To The 
Poincare' Sphere Representation of Polarized Forms 

Consider the polar coordinate system (r, 5), Fig. 56. 
Let r = J = (Cos 2 2© + Sin^G Cos^G)^ 

then r = Cos 2<a (5.62) 

- using Eqs. (5*41) and (5.4-2) 

Therefore, the magnitude of the radius vector r is 
a measure of the ellipticity. So, the curves of constant 
elliptic it y are concentric circles with center at 0, the 
origin. 

Let =< = 2=< 

where °c is the azimuth 

then r Cos°f = J Cos 2<=c (5*63) 

But from Eqs. (5*4-1 ) and (5 -^3 ) 

Cos 2* = Cos 2.yj (5.64) 

s o that 

rCos^ =' C os 26 (5.65 ) 

i.e. the x component (= rCos^) depends on e only. So,, the 
loci of points of constant amplitude ratio are straight 
lines perpendicular to the X-axis. 
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Mow, 


x = r Cos * = Cos 2 6 ■ 

7 = r Sin oC = Cos 2w Sin 2* 


• , • 


x 2 + ,_zL 
Cos 2 £ 


Cos 2 2 

1 


Cos^2tt Sin 3- 2<* 

C os 2 £ 

- using Sq. C5* i +2) 


( 5 * 66 ) 


Eq. (5*66) represents a family of ellipses with 
their centers coinciding with the origin and with the 
lengths of semi -minor and semi -major axes as Cosg and unity 
respectively. So the locus of points of constant phase 
difference is an ellipse belonging to this family of 
ellipses. The Wulff's grid (stereographic projection) 
method can also be derived in a similar way without refer- 
ring to Poincare' sphere by choosing a proper expression 
for r, the radius vector. 

The graphical methods described can be used for 
values of phase difference and azimuth beyond the range 
0 - /2. The first graphical method is plotted in Fig. 57. 

for tho ranges HtQr, 0<£<2ir 

Some Examples of Using the Graphical Method . “ 

We shall now take up the two examples of Schwiegef 
for finding the characteristic parameters. 
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3x-npl3 (a) 

Consider the combination of two models (linear 
retarders). Model (1) produces a retardation of A 1 =45 a e r 
Mode a (2) produces a retardation of & 2 = 60 dog. The angle 
t -t “-s.i '-ne axes of these two models is craven bT 


P 12 = ^5 deg. 

ihe characteristic parameters of this system are 
found by one of the methods described earlier, i.e. the 
secondary characteristic direction p 2 and the characteristic 
retardation S are found from the exit light ellipse parame- 
ters when the incident light ellipse is circularly polarized. 
The primary characteristic direction P 1 is found by deter- 
mining the azimuth of light at the entrance which produces 
plane polarized light along one of the secondary principal- 
-stress axes at the exit of the system. This is found by 
tracing back from the exit point to entrance point of the 
system. This is demonstrated in Figs. 58 (a) and 58 (b). 

From these graphical methods, 

5 = 69 . 4 - deg . 

*1 (rotatory power) = 13.8 deg. 

Orientation of primary ) 

characteristic directions with ) 

respect to the axes of the ) ~ deg. 

) 

first model ) 
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Schwi eg er 1 s value s ; 

+1 = ^3*5 leg. 

€ =68 deg . 

r =33.5 deg. 

Example (b ) 

Consider the case of constant principal-stresses 
rotating uniformly along the wave normal. 

Data: 


The phase angle introduced due to the difference 
of principal-stresses without rotation is e = iOO deg. 

Rotation angle of the principal-stresses betweetj 
th*se at the entrance and at the exit 0 = *+5 deg. Let us 
divide the continuous system into ten plates. Then each 
plate has the following data: 

The phase angle of cne sub plate = = 10 deg. 

Rotation angle between two sub plates = A 0 = 5 deg. 

Schwi eger's solution: 

Rotatory power = = 

Primary characteristic direction 
with respect to the principal- 
-stress axes at the entrance 

Characteristic retardation =5= 86 deg. 


-13 deg. 

) 

) 

) = P 1 = 29 deg. 
) 
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Results from the second graphical method: 

(The procedure is demonstrated in Fig. 59.) 

= - 10*8 deg. 

f -| = 26.7 deg. 

S =84- deg. 

iheoretical Values using the analytical solution 
(Chapter 2) : 

= - 9*8 deg. 

P 1 = 27. k deg. 

5 = 78.6 deg. 

In this section r -re have developed two new graphical 
methods to trace the passage of polarized light through a 
series of optical elements. These two methods are compared 
with a few of the existing graphical methods. Procedures 
for the use of the new graphical methods are described with 
application to various polariscope arrangements. The first 
method traces the transformation of polarized light through 
the same constant ellipticity and constant amplitude ratio 
curves. The second method also makes use of the same 
principles, but involves only straight lines in the graph 
which are convenient from the construction point of view. 
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33 * THEORY of compensation 

In this section, we shall first present a brief 
review of the existing methods of compensation and the methods 
of measuring the integral fringe orders. Then, we shall des- 
cribe, the development of a compensator which has certain 
advantages over the previous ones, and a general method of 
measuring the integral fringe orders. The theory for these 
methods makes use of the concept of optically equivalent sys- 
tem and various methods of determining this system for a 
series of retarders which were described in Chapter 3* 

The compensator can be defined as an optical device 
which introduces a required amount of retardation along a 
particular set of axes. This is equivalent to a two-dimensio- 
nal photoelastic model in which the principal-stress dif fere- 

45 

nee can be varied at will. Jerrard has given an account of 
different types of optical compensators in vogue. 

One of the commonly used compensators in photo- 
elastic investigations is the Babinet-Soleil compensator. 

This consists of two quartz wedges of same inclination with 
their fast axes parallel and a rectangular block of quartz 
in series with these wedges with its fast axis perpendicular 
to that of the wedges. The retardation is varied by the 
lateral movemoit of the wedges and thus their combined thick- 
ness in the optical path. Other types of compensators inclu- 
20 

de Coker compensator which is a tensile strip (preferably 
made of the same material as that of the model) in which the . 
tension can be varied by a loading frame, a model beam 
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compensator, (Jellyman and Milne 46 ), in which the fringes 
resulting from lending are frozen into the material; mica 
plate compensator, (Phenister in which the angle of in- 

cidence of light is varied to affect a variation of the re- 
tardation introduced; etc . The compensation methods which 
donot involve the use of any additional optical elements are 
the Tardy compensation, De Senarmont compensation and the 
inverse Tardy method (Chapter 2). The underlying principle 
in all these methods of compensation is to introduce a reta- 
rdation equal and opposite to that introduced along the prin- 
cipal-stress direction (or the characteristic direction in 
a three-dimensional case), so that total extinction of light 
is obtained when observed in the polarise ope. The retarda- 
tion introduced can be related to the (secondary) principal- 
-stress difference in the model. 

Some of the desirable characteristics of a compen- 
sator can bo listed as follows: 

(1 ) the principle of operation should be simple and the re- 
lation between the movement (angular or linear) in the com- 
pensator and the retardation that it introduces should be 
linear, (2) it should be inexpensive, (3) be accurate, (4) 
be useful for all wavelengths of light in the visible region 
(5) have a large field so that it can be used directly in 
the diffused light polarise ope and (6) its retardation should 
not vary with time. 

The Babinet -Soleil compensator involves a very 
careful preparation of the wedges and so is expensive. 
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But, it can be used with whits light to find the integral 
fringe order upto about six fringes. But it does not give 
correct results when the variation of the material fringe • 
constant of the model with wavelength is not linear (assu- 
ming that it is so for the quartz wedges and the rectangular 
quartz block in the compensator). Its field of view is very 
small and is about 1 cm. in diameter. Coker's compensator 

If 7 

is cumbersome because of the loading frame. Phemister’s 
mica plate compensator has a nonlinear relationship between 
the angle of light incidence and the retardation introduced 
in the model. It cannot be used to find the integral fringe 
order directly. The compensators which involve freezing of 
stresses into a piece of photoelastic material usually exhi- 
bit time dependence of the retardations. The Tardy, Be 

f 

Senarmont and inverse Tardy methods can not be used with 
white light to find the total fringe order directly. They 
involve some errors if the wavelength of light is changed 
since the quarter-wave plates that are used show some depen- 
dence of their retardation with the wavelength of light. 

It is the objective of this section to develop methods of 
compensation which involve minimum of the above drawbacks 
using the theory of characteristic parameters discussed in 
Chapter 3* 

28 , . 

Pancharatnam , in his analysis on achromatic cir- 
cular polarizer and achromatic quarter-wave plate suggested 
that two quarter -wave plates with a half-wave plate in bet- 
ween will act as a compensator. Tuzi^9 suggested that a 
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combination of two retarder plates of equal retardance can 
be used as a quart e r -wav a plate. But this combination will 

40 

have a rotatory power. Childs ~ ust»d the combination of 
3 retarder plates to act as a quarter-wave plate with no 
rot at or v- power. The present compensator to be described is 
essentially based on these principles. The analysis for 
this will be made using the Stoke's vector. 

Combination of Two Linear Retarders of Baual R etard ance p 
Consider two linear retarders each of retardance p 
(with the angle between their axes being ©]. The transfor- 
mation of the Stokes vector through this system is given by.* 
(bot V Q and v be the Stokes vectors of light at 
the entrance and exit of the first retarder, with respect to 
its axes} -VI-/.., 



- 


w 


!**“ ~ ^ 


s 1i 


S 1i 


; Si | 

v o = 

S 2i 


S 2i Cos p - S^ Sin p 


s 2 


_ S 3±_ 


3^ ^ C os p + S Sin p 


_ S 3_ 


Let and V 2 be the Stokes vectors at the exit 

of the first and second retarders referred to the axes of 
the second retarder. We then have. 


V 




s 1 Cos 2 0 + S 2 Sin 29 


S^Cos 2 e + S 2 Sin 2© 

S 2 Cos 29 - S 1 Sin 29 


(S 2 Cos 2© . 


- 

- S-| Sin 2 9 X 5 os P — S^Sin p 

_ S 3 


Cos p+ Sin p 

= 


(S^ Cos 2© - S^ Sin 2©) 
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(3 1 C os 26 + s 


+ S 2 Gin 26) Cos 29 


)in 


GGS p 


(S 2 Cos 26 - 3 1 Sin 26) - S 3 Sin p) 
[ S 2 c °s 28 - S 1 Sin 20 ) Cos p - 80 5 


- Sin p |Cos 2© 


+ Sin 26 (s 1 Cos 29 + S 2 Sin 29) 

S 3 Cos p 4- Sin p (S 2 Cos 2 © - S-| Sin 2©) 


s 2 = V 


(5.6?) 


where V is the Stokes vector of light at the exit of the 

second, retarder with respect to the axes of the first 
retarder. 

If V 0 and V 1 are the Stokes vectors of incident light 
referred to the axes of the first retarder (reference axes) 
and the retarder of the optically equivalent system, x;e have 
S 1 Si C os 2 pj + S 2 Sin 2p 1 " 

V 0 = S 2 ■ - s 2 Cos 2 Pi - S 1 Sin ' = (Vj) 


S 1 60s 2 + S^ Sin 26^ 

( S 2 Cos 2? i " s i Sin 2^ ) Cos 5 - S^ Sind 
! s 3 Cos5 + Sin« (S 2 Cos 2P 1 - s 1 Sin 2p 1 ) 
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where 7„ 
the retarder 
its axes. 


is the 
of the 


Stakes vector of light at the exit of 
optically equivalent system referred to 


sysa 


r 


S’- SCO 


-i CD 


z n 


system is equivalent to an optically equivalent 
single retarder (retardations, aximuthp 1 with 
axe£ the first plate) and a rotator with 


a rotatory p.-^er of . The 
vector through the optically 


transformation of the Stokes 
equivalent system is given by 


Also, 



Cos 2v (S 1 Cos 2p 11+ S 2 Sin 2*-, ) - Sin 2i> £ Cos 8 
(S $ 2P, - Sin 2 Bq) - S Sin 5 j 

Cos 2# [ (S 2 Cos 2P 1 - S 1 Sin 2^ ) Cos 5 - S 3 Sind ] 

+ Sin 2t (S 1 Cos 2p 1 + S 2 Sin 2$ 1 ) 

S' 3 Cos <5 + Sin 5 (S 2 Cos 26-, - S 1 Sin 26 1 ) 



( 5 . 68 ) 



where ip = P, + £-j and V is the Stokes vector of the light 
at the exit of the optically equivalent system (after the 
pure rotator) with respect 'to the reference axes. 

We shall now find the parameters f , 5 and * of 
the optically equivalent system in terms of p and 0 using 
Eqs. (5-67) and (5.68). 
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.bet the incident light he circularly polarized for 
the original system and the equivalent system. The Stokes 
vector of this incident light will be (0,0,1 ). Equating the 
Stokes vectors at the exit of both the systems and substi- 
tuting for the incident Stokes vector, the vector (0,0,1 ), 
we obtain 

S^ = Cos 2 P - Sin 2 p Cos 2© = Cos 5 (5*69) 

S'] = - Sin pSin 2© Cos 2© + Sine Cos p Cos 2© Sin 2© 

+ C os p Sin p Sin 2 © 


= Sin 6 Sin 27 


(5.70) 


How, let the incident light be plane polarized along 
one of the reference axes (i.e. the axes of the first retar- 
der), such that V Q is (1,0,0) for both the original and 
optically equivalent systems. Then, by equating the Stokes 
parameters of light at the exit of the two systems, we obtain 


83 = - Sin 2© Sin P = - Sin 2p 1 Sind (5.71 ) 


If 5 is it /2 and the value of p is given, the value 
of© can be calculated, i.e., the value of ©for ^iieh the 
system acts as a quarter -wave plate. Pr® Eq. (5.69), for 
d =ir / 2, we have 

tan 2 p = Sec 20 (5*72) 

Since Sec 2 > 1 > tt/4 

Tan^ p M or Pl’/t (5.73) 


Also, from Eq. (5*71), the value of 
Sin 2© Sin p = Sin 2^ 


ft| can be foundas 
( 5 . 74 ) 



From 3 q. (5-70) 


Sin 2 y = Sin p Sin 2 0 ( - Cos 2 e + Co.sp Cos 2e + Cos ?) 

( 5 - 75 ) 

from which the values of *p and ^ c an be calculated. 

Combi nation of Two Quart or -wave Plat eg 

ce shall now investigate the behaviour of a combi- 
nation jf two quarter-wave plates and its use as a compensa- 
tor. Substituting p =7r /2 in Sqs. (5.69) and (5*71), we 
obtain 

- Cos 2 8 = C cs5 

or 0 =ir - 20 (5.76) 

i.e. the characteristic retardation is linearly related to 
Q , the angle betv/een the -axes of the quarter-wave plates. 

From 3q. (5.71 ) 

2 p 1 = p = ir/2 

or P 1 = v / 4 (5-77) 

From Eq. (5*70) 

- Cos 2© = Sin 2 ip 
or ip = e - v /If 

= P 1 + *1 
= TrA + ^-j 

or *1 = © - ir/2 (5-78) 

i.e. the rotatory power is also related to© , the angle bet- 
ween the axes of the two quarter -wave plates. 

Hence, two quarter -wave plates with angle © between 
their axes are equivalent to a retarder ( a xim uth p = ir /if 
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with respect- to the axes of the first quarter -wave plat e , 
retardation 5 =r r - 2 8) and a pure rotator of rotatory power 

$ -j = (8 - ir/2). 

So this combination can be used as a compensator 
if this rotatory power is nullified. This can be done in the 
following way. Consider the system of a light source (unpo- 
larized), a polarizer and two quarter-wave plates, Fig. 60. 

In Fig. 60(a) - 

OP - Axis of the polarizer 

OF-j - Fast axis of the first quarter -wave plate CR - 
one of the axes of the retarder of the optically equivalent 
system (one of the primary characteristic directions) 

0F 2 - Fast axis of the second quarter -wave plate 
OM - One of the secondary characteristic directions 
Also, 

— IT — 7T — 7T IT 

Y= e- £ - = 9'q: - e + 2 = IJ 

i.e. the angle between the secondary characteristic direc- 
tions and the axes of the second quarter-wave plate is v /+ 
which is a constant. So, with this combination of polarizer 
and quarter -wave plate , we can produce the required ampli- 
tude ratio and phase difference along, the secondary charact- 
eristic directions (say OM and ON). Along these directions, 
the amplitude ratio will be tan © and the phase difference 
will be ( ir- 2 8 ). If the amplitude ratio is to remain cans-' 
tant, but the phase difference is to vary, we rotate the po- 
larizer and the first quarter -wave plate together by the 



V F~, 

\ _ / ‘ 

/ C* ,.x 

\ f ' 'V x'‘ 

\ i / V 



FIG. 61(a) 


Orientations of the axes of the four 
retarders 

F li F 2> F 3’ F 4 are the fast axes of th<2 
first ^second, third and fourth retarders 

respectively . 

R-|,R7 are the primary characteristic 
directions of the first and second systems* 
Mi, M2 are the corresponding secondary 
characteristic directions. 




FIG. 61(b) 


FIG. 61 (c) 
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required amount. If tho phase difference is to be increased 
by p we rotate the polarizer first quarter-wave plate combi- 
nation through an angle - c / 2 , keeping the second quart er- 

-v;ave plate fixed. 

Similarly, the combination of two quarter-wave 
plates and analyser can be used as a compensator, Fig. 60(b). 
In this case, the second quarter-wave plate and analyser are 
to be rotated together so that retardation is varied along 
the fixed direction OR (the primary characteristic directions 
remain fixed in space if the first quart er -wave plate .axes 
are fixed). The retardation added along the primary charac- 
teristic directions will again be (v - 2e). By altering the 
value ofe , we can obtain the required amount of retardation. 
Rotation of the analyser is necessary to compensate for the 
rotation of the secondary characteristic directions with the 


second quarter-wave plate. A quarter-wave plate will have 
the same effect when it is rotated through ir radians. Its 
effect will continuously vary as it is rotated through 0 - ir 
radians. So, the angle e can take values from 0 to w i.e. 

We can introduce retardations from v to - v (Eq, (5*76)) 
radians by these compensation techniques. 

Now, consider two pairs of retarders, each pair 
consisting of two retarders of equal retardation p. The 
axes of these retarders and the primary and secondary chara- 
cteristic directions of the corresponding optically equiva- 
lent systems are shown in Figs. 61(a) and 61(b). 
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F 1 J F 2 ’ ? 3 » F l+ are the fast axGS 01 ” the first, second, 
third ana fourth retarders respectively. 

R 1 , 1 2 are the primary characteristic directions of the 
i ir st and second systems 

-'*1 j ^ are corresponding secondary characteristic 
direct i ons . 

How examine Eqs. (5-72), (5.?3) and (5.7*0. In 
these equations if is replaced by - 0^ we obtain -? 1 
for P -j , ~'-P for ^ and - for (*.* ^ = i> - F ^ from Eq. 

(5.68). i.e. if the angle' between F-j and F 2 is 6^ and the 
angle between F^ and F^ is -©^, the rotatory powers of the 
equivalent systems for the pairs F 1 , F 0 and F^, will be 
equal but in opposite directions. Also, the optically equi- 
valent system can be either a retarder and a rotator or a 
rotator and a retarder. In both cases, the rotatory powers 
and characteristic direct ions will be equal. But the retar- 
der axes will be along the primary characteristic directions 
in the first case and along the secondary characteristic 
directions in the second case (Chapter 3). Let the axes of 
the four retarders be oriented as shown in Fig. 61(c). 

Let the combination of the first two retarders be 
equivalent to a retarder (azimuth *f-j with respect to the 
axes of the first retarder of retardation v/2) and a pure 
rotator of rotatory power ^ . Then the combination of the 
third and fourth retarders will be equivalent to a system 
with a retarder (azimuth - with respect to the axes of the 
third retarder of retardation ir/2) and pure rotator of 
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rotatory power - ^ y or a system with a pure rot at or of 
rotatory power - and a retarder (azimuth with 

respect to the axes of the third retarder of retardation 
7r /2). So, when the four retarders are arranged in sequence 
as shown m Fig. 61 . the rotatory power of the first two 
retarders will be cancelled by the rotatory po-er of the 
combination of the third and fourth retarders. And the com- 
bination of the four retarders acts lie e two quarter-wave 
plates. The orientation of the axes of these quart or -wave 
plates will be along the primary characteristic directions 
for the first two retarders and the secondary characteristic 
directions of the combination of third and fourth retarders 
respectively. The orientations of these directions are 
found from Eq. (5 .A) and these are fixed relative to the 
axes of the retarders once the values of 0 and p are fixed 
(pis the retardation of each of the retarders). Thus the 
combination of four retarders each of retardation (p" > 
from Eq. (5 - 73 ) ) can be used as a combination of two quarter- 
-wave plates. This in turn can be used as a compensator as 
described above. By varying the value of 8^ the angle 
between the retarders, this combination can be used as a 
compensator for any wavelength of light. The range of wave- 
length is however limited by the value of pT and the dispers- 
ion characteristics, (i.e. the variation of the retardation 
with wavelength) of the material of the retarders. 

In the above discussion, we have shown that the 
combination of two quarter-wave plates can be used as a 
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compensator provided the rotatory power introduced by the - 
system is nullified by rotating the polarizer or analyser 
as the case may be. To develop a compensator which does not 
introduce any rotatory power i.e. to make it act as a single 
retarder, let us consider the case of a half-wave plate placed 
between two quarter -wave plates. 

A Half-Wave Plate Betwe e n Tw o Quart er -Wave Plates 


It will now be shown that a system consisting of a 
quarter-wave plate, a half-wave plate and a quarter-wave 
plate in that order, can be used as a compensator. The axes 
of the two quarter-wave, plates are parallel. The angle bet- 
ween these axes and the axes of the half-wave plate is © • 
The transformations of Stokes vector through such a system 
are given below 


Let Vq and be the Stokes vectors of light at the 

entrance and exit of the first quarter -wave plate referred 

\ 

to its axes. Then, we have, 





S 1 " 

v o * 


. — > 





s 2 


If V 2 and V 3 are the Stokes vectors of light at the 
entrance and exit of the half-wave plate referred to its 


axes 
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S„ Cos 26 - Sin 26 

1 3 


S Cos 26 - S_ Sin 2 e 
1 3 


-S 3 Cos 26 - Si Sin 28 S 3 Cos 26 + Si Sin 2© 


Also, 

Cos 26 (S 1 Cos 2 © - S 3 Sin 2 © ) - Sin 2 © (S 3 Cos 2 © 

+ S-j Sin 2© 

= Cos 26 (S 3 Cos 2 © + S 1 Sin 2 ©) + Sin 2 © <S 1 Cos 2© 

- S 3 Sin 26) 


Cos 26 (S Cos 2© - S Sin 26) - Sin 2© (S_ Cos 2 © 
13 3 

+ S 1 Sin 26 ) 


Cos 2 © (S^ Cos 2© + S^ Sin 2© ) + Sin 2©(S^ Cos 2© 


- S 3 Sin 2 6) 



( 5 . 79 ) 
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sphere , V^: Stokss vector of light incident on the second 
quarter -wave plate with respect to its axes. 

V 5 : Stokes vector of lj -g^ emerging from the second quart er- 
-wave plate with respect to its axes. ' 

The Stokes vector at the exit of the optically equivalent 
system is given by aq. (5*68). Comparing this equation with 
(5*79) for S.j = S 2 = 0, = 1 , we obtain 

S 3 = Cos5 = Cos 4-6 or 5 = If 9 (5.80) 

S-j = Sin5 Sin 2\p = - Sin 4-0 


or Sin = - 1 


5 = 4- e 



II 

• 

CD 

• 

r A 

( 5 * 81 ) 

for S 

= 1 , s = s 

2 3 

= 0 


s 3 = - 

Sin 2p^ Sins 

= Sin 4-e 

( . . q = 4 - 0 ) 

• 



0 

hi 

*0 

II 

1 

^A 

(5.82) 


* * The rotatory power of the) 

_ ) = ♦ = t - P . = 0. 

combination ) 1 1 

So, this system acts as a single retarder with retardation 
4-0 and azimuth - ^A with respect to the axes of the first 
quarter-wave plate. The retardation introduced can be 
varied by varying the angle © . Thus, this system can be 
used as a compensator. 

The whole system, i.e. a quarter-wave p la te, a 
half-wave plate and a quarter-wave plate, can be constructed 
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with eight retarders each of retardation p, when tan p>1 
(for tij.e Particular wavelength of light under consideration ). 
This can be done as follows. 

Consider two types of quarter-wave plates ( a ) 
and (b) each of which is a combination of two retarders. 

Ine orientations of the axes of these retarders are shown 
in Fig. 62. In the (a ) -type combination, the rotatory power 
is clockwise where as in the (b ) type it is anticlockwise. 

The dilference between these two types of combinations is 
in that the angles between the axes of the two retarders in 
types- (a) and -(b) is 0 and - ©respectively. If the type-(b) 
combination is placed after the type-(a) combination, the 
rotatory powers of the combinations cancel each other so 
that the two combinations behave as two quart er-wave plates. 
So, to construct a compensator we can use four such combi- 
nations (combination of two retarders each of retardation p ) 
with tne first and third combinations of the type (a) and the 
second and fourth combinations of the type (b). Then the 
rotatory powers of the first and third combinations are 
nullified by the rotatory powers of the second and fourth 
combinations respectively. Thus, the whole system of these 
eight retarders acts as a series of four quart er -wave plates 
which in turn can be used as a compensator as described 
abdve. By varying the angle e, the angle between the axes 
of the two retarders in a two retarder combination, accord- 
ing to Eq. (5*72) i.e. tan 2 p =. Sec 2 6 , we can use the com- 
pensator for any wavelength of light for which the condition 



220 


tan p >1 is satisfied. (This, in practice is not a severe 
restriction since the retardation P does not vary appreciably 
within the visible region of the electromagnetic spectrum 
for most of the photoelastic materials). Using this method 
large field compensators can be made using large sheet type 
retarders (ex. frozen-stress sheets). The angle 0 can be 
found experimentally for different wavelengths of light with- 
out any calculations. It is advisable to cut all the eight 
retarders frbm a single large retarder of uniform retarda- 
tion. The use of a large field compensator in a lens polari- 
se ope eliminates refocussing of the lenses when the -compensa- 
tor is introduced into or removed out of the field of view 
of polariscope, which is necessary with a Babinet-Soleil 
compensator because of its small field. The compensator 
described in this section thus has some additional advantages 
like large field, cheapness etc. However* the integral fringe 
order cannot be deter min ed directly with this compensator 
as can be done with Babinet-Soleil compensator (upto a 
certain number of fringes, usually about 6). However, the 
integral fringe orders upto virtually any timber can be 
determined by an altogether different method. We shall new 
develop this method in the following section. 
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C . Dig IRM STATION OF INTBORAL FRINGE 0HD_JRS III PHOT OELA.ST IC XT Y 

In analysing a photoelastic coc-cimen having no 
stress-free corner or, in which the formation of fringes 
during the loading process cannot be observed, it is desi- 
rable to have some means of determining the integral fringe 
order at a point when the fractional fringe order at that 
point has been determined by the Tardy or other method. A 
Babin et-Soleil or a similar compensator can be used to dete- 
rmine the integral fringe order upto a certain value. How- 
ever, there is a serious .point which needs attention when a 
Babinet -S oleil compensator is used. In practice, to deter- 
mine the integral fringe order at a point by the help of a 
Babinet -S oleil c ompensat or , a white light source is used and 
retardations are added or subtracted until a zero-order 
fringe (a black fringe) appears on the screen. This process 
assumes that the compensator acts as a continuation of the 
model and its behaviour is similar to that of the model. 

This implies that the model and the compensator introduce 
retardations which are equal in magnitude but opposite in 
sign for all colours contained in white light. This is true 
only if the material of the model and that of the compensator 
exhibit same type of dispersion characteristics (variation 
of retardation with wavelength of light). This is not the 
case in general. Most of the compensators are supplied 
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with a calibration chart which sho^s the dependence of the 
compensator on the wavelength of light. Further, experience 
has shown that using a Sabinet -Soleil compensator and a 
photometer, instances occur where the reading on the photo- 
meter is almost zero for two consecutive fringe orders, mak- 
ing it difficult to identify which one corresponds to the 
zero fringe order. This experience is quite noticeable with 
a urethane rubber model. 

In 1963, Pant ^ showed that by determining the 
fractional fringe orders with two wavelengths A-j and A£, 
one can determine the fractional fringe order obtained by 
using a third wavelength A , His analysis ignores the de- 
pendence of the stress -optic coefficient A and also is res- 
tricted to the case where the difference between the fringe 
orders associated with A^ and .A^ is only a fraction, i.e., 
the integral fringe orders associated with a.j and A 2 are 
equal. 

We shall now show, how by using two wavelengths, 
one can determine the integral values (and their multiples) 
for both wavelengths, associated with the fractional fringe 
orders observed at any point. The solution can be obtained 
by using a nomogram or constructing a ready reckoner for any 
two wavelengths A^and A 2 as shown in this analysis. Fur- 
ther, the usual assumption that the stress -optic coefficient 
is independent of wavelength, is not used. Its dependence 
on A is explicitly shown hy experimental results. 



Stress -Op tic Law 
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xhe stress-optic law can be expressed a 

T 1 ^ - q) = (p - q) 

wiie re p • retardation 


( 5 . 83 ) 


h 


(P - q) 


thic m-jss of tne model or length of light path 
inside the model 

principal or secondary (for frozen stress slices: 
principal-stress difference 


>. : wavelength of light used 
C : Stress -optic coefficient 
f : Model fringe constant in tension 


m practice, the stress-optic coefficient C is 
considered to be independent of wavelength. However more 

careful observations have shown 20 . 49 that 0 is influenced 

to some extent by the wavelength. For our present analysis, 
we shall assume that C depends an X and hence, f is a func- 
tion of x . if, ln any experiment, C is found to be inde- 
pendent of x , the analysis remains unchanged. 

Let x, and x £ be the wavelengths of the light used 
P 1 , P 2 the corresponding retardations. In general, 
these retardations will consist of an integral value and a 
fractional value. Let ,, = (n, + ) and P 2 = („ 2 + 

«here ^ , n 2 are positive integers and S 1 , S 2 are positive 
.£>*. « . 


fractions less than unity. Without loss of generality, 
let f, be greater than f 2 , so that their ratio a - 5. 

~ f 2 


>1 



From Eq. (5.83) 

(n-j +3-]) = ^ ( p “ q) 

( n 2 + 02) = 4 ( p - q) ' 

” L 2 

n 2 + 02 _ £l 
n 1 + 51 f 2 
or 

n 2 * 2 n 1 = 2 5i - 5 2 = 0 ( 5 . 8 +) 

If the stress-optic coefficient happens to be inde- 
pendent of x , then, 

z = z*= ^/X 

The experimental investigations will show that using 

Jft 

z instead of z will lead to erroneous results. In practice 
the fractional fringe orders 6 1 and 5 2 are determined and 
hence £ becomes an experimentally determined factor. Eq. 

(5 •8b-), will- then have to be solved for n 1 and n 2 subject 
to the condition that these are whole numbers. 

General Observations : 

a) Since z > 1 (as experimentally observed), P 2 > pj . This 
is obvious when C is independent of X , 

Hence, 

^2 ^2 ^ ^1 + 5 -j 

and n 2 > ix, (5*85) 

b) Since 0 S-j < 1 and 0<5 2 <1, for a given z, 
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and 


e 

max 

= cz s i 

) 

max 

<S 2 ) . 

^ mm 

0 . 

mm 

11 

tQ 

Op 

^min 

v ‘^2^max 

i.e. 

z > e 

> - 1 



( 5 . 86 ) 


c) In Eq. (5.84), if we add N to n 2 and N/z to where 
N and N/z are integers, the equation is still satisfied. 
Let N be the smallest integer such that N/z is also an 
integer. Then the solutions for n 2 and are 


(n^ + SN) and (n^ + SN/z), where S = 0,1 ,2, 

(5-87) 

In other words, if we find two basic solutions n 2 
and n^ , satisfying Eq. (5*8+ ), then the multiple solutions 
are given by the above equation. 


Construction of Hemogram : 

The wavelength in visible spectrum varies from 7 700° A 
(extreme red) to 3900°A (extreme violet) and for this range 
z = -j~ < 2. In practice, the usual range is from red (about 
7500 °a) to blue green (about 5000 °a) giving for z a mlue 
less than 1.5* 

In order to solve Eq. (5*84), one can construct a 
nomogram as shown in Fig. 63. line AF represents z, line 
AD, the integral values of n^ and line BE, the integral 
values of n^ . The scales for the values of z, n 2 and 
are the same. From A, rays A^ , A 2 , A^ , etc are 
drawn passing through the integral values of n-j = 1 ,2,3, . • .etc. 
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The slopes of these rays are 1 ,2,3, . ..etc., by construction. 

Let us say that in a particular example, the ratio 
of the two model fringe constants corresponding to , x 2 , 
is f-j /?2 = Z'l * Fig, 63. A G = z-j . A horisontal line G H 
is drawn cut sing the rays AA 1 , AA 2 , AA-,....etc. at g,N,G... 
etc. as shown. As the slopes of these rays are 1,2,3,... -etc 

G K = z 1 ; G N = 2 z 1 , G 0 = 3 z 1 , etc. Line G H cuts 

the vertical lines representing the integral values of n 2 
at L, II, P, 0,...., etc. Distances L X, M K, P N, T II, etc 
represent the values of £ that one determines experimentally. 
The procedure for using the nomogram is as follows; 

For a given problem, let the measured values of fra- 
ctional fringe orders be 5| , ^ so that 2 5-] - s 2 = Gl • 
Identify a point which is the intersection between the line 
G H and one of the vertical lines representing the values of 
n 2 (like L, M, P, T,...) such that the distance between this 
point and the slant ray immediately to the left of it is 
equal to G-j . Let us say in our case, M is such a point 
where M K =G -j where K is a point on the ray immediately to 
the left of M, i.e. ray A^ . The vertical line on which point 
M lies gives the value of n 2 (in this case n 2 = 3 ) and the 
slant ray immediately to the left gives the value of n 1 
associated with it (here has n-j = 1 ). The proof for this 
is evident from the figure, since MG -KG = M K, i.e., 
n 2 ~ n 1 2 1 = ^1 f (5«84-), where n 2 = M G = 3 and z-j n-j 
= K G = a-j with n-j = 1 . 
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Generally for a given value of £ 1 , there will be 
more than one point ( lying on one side of the vertical lines), 
such that the distance between this point and the slant ray 
immediately to its left is equal to G 1 * These give multi- 
ple values according to Eq. (5.87). This aspect will be dis- 
cussed later. > 

One observes from Fig. 6 ^ 3 . the values of £ can only 
be discrete. This is evident by the nature of Eq. (5.85). 

For, we can have a problem where n g* and z are at our 

command, but n 1 and ^ have to satisfy Eq. (5.8f) and also 
that n^ should be an integer. 

In the above procedure, we have assumed that € is 
positive. If £ turns out to be negative say , € = - G 2 ? then 
we measure the distance from the point such as M to the slant 
ray immediately to the right ; i.e, if M is the point such 
that M N (where N lies cn the ray immediately to the right) 
is equal to f - 62 ! » then ^ = 3 and n-j = 2 , since the ray 
on which N lies is with associated =2. For negative 
values of £ , the point can also lie cn the segment 0 G* 
Negative values cannot be numerically greater than 1 as shown 
by our previous analysis . . 

To nake the above procedure clearer, let us take 

a specific example. Let z = 1.25. A horizontal line pr 

corresponding to this value is drawn as shown in Fig. 63, 

cutting the vertical lines at a , a , a etc. Let 

* 2 3 

us assume that we find experimentally z S-j - 5 2 = 0 , 25 t 
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Now we look for a point lying on one of the vertical lines 
( such as a.j or or or...) such that the distance bet- 
ween this Point and the slant ray immediately to its left is 
0 . 25 * In our case is such a point since b^ = 0 * 25 . 
Hence n 2 = ^ ( point 34. lies on this) and = 3 (since b^ 
lies on ray A3). We can see that these values of Oj and n 2 
satisfy Er. (5*8+) with z =1.25 and J = O.25. On the other 
hand i_ c happens to U - 0.5, we select a point such that 
tne distance between this point and the slant ray imnedia- 
tely to the right is equal to 0.5. Hence, we get a point 

a 2 s -*- nce a 2 ^2 = ^*5 and b 2 lies to the right of a 2 * 

Then n 2 = 2 and n-j =2. 

A value of G= - 0.25 will give us points a^ and 
b 1 with n 2 = 1 and n 1 =1; a value of £ = - 0.75 gives us 
points a3 and b3 with n 2 = 3 and n, = 3; a value of q = 0 
gives us point a^ (b^ coincides with it) with n 2 = 5 and 

n.j = 

Multiple Values for rtj and n 2 : 

From Fig. 63 . we notice that for e = - 0 . 25 , we get 
points a-j and b-j as noted above with n 2 = n-j =1, and also 
points a^ and bg such that b^ = |- 0.251 with associated 
values n 2 = 6 and ri| = 5 * These also satisfy Eq. ( 5 « 8 N-) f 
since, 6 - 1 .25 x 5 = - 0 . 25 » This is in agreement with 
Sq» ( 5 * 8 7 ) which gives all multiple values. For our case, 
the lowest integer N such that N/z is an integer is 5. The 
solution for n 2 and n-j according to Eq. ( 5 . 87 ) are: 
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n 2 + S N = 1 + 5 S = 1 ,6,11,16, , etc . 

and n-j + $-1 _ -] + 5+ s= 1,5, 9, 13 ,...., etc. 

2 

all satisfying Eq. (5 .84- ) . This should not cause much of a 
problem in practice, since the difference between one set 
of values for n-| ,np and the subsequent set is so large that 
an analysis of the problem can give a clue as to the appro- 
priate value. Further, as will be shown later, one can 
always choose two- appropriate wavelengths such that the 
difference between two multiple values is so large that 
there will be little doubt as to the correct value. If 
ambiguity still exists, one can use the franctional fringe 
order 5^ for a third wavelength of light That is, one 

set of multiple solutions are found using the -values of 5-j 
and 5 2 anc * a second set is found using the- values of 3 2 
and 3^. The value of common to both these sets gives 
us the correct solution. For example let these two sets of 
multiple solutions be given by: 

(n-j , n 2 ) = (1,2), (6,6), (11,10), (16, 14), (21 ,18),.. .etc. 
and 

(n 2 , n 3 ) = (0,1), (3,3), (6,5), (9,7), (12,9), ....etc. 
then the correct values of the integral fringe orders are 
given by: 

(n^ , n 2 , n^) = (6,6,5). 

All Possible Values for £ And Associated n 2 , Values: 

In practice it is not necessary to use a nomogram. 
For a given ratio z of -the wavelengths, one can write down 
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all b he possible values of G that one can determine experi- 
mentally and write down the corresponding values of n-j and 
n 2* Tilis P roc edure has several advantages, since one has 
a sort of a ready -reckoner using which, the associated 
values of n 2 and a, with any given G can be read out. Also, 
if duo to some experimental error, the value of G is not de- 
termined exactly, the ready reckoner can immediately tell 
what should be the correct value of G for any z and a given 
range of n 2 , n 1 . This is so because, for any z, the values 
of G can only be discrete as was pointed out earlier and as 


seen in Fig. 63 . and also for any range of np , 
between any two successive rays JL A or A A. 
etc.) the values of € are fixed. 


n-| ( lying 

or ^3 > * * * 


In Fig. 64. , a portion of Fig. 63 ., is sho^m vjLth 
G H, the horizontal line cutting the rays A 1 , A 2 , A 3 , . . . 
etc. at distance of z, 2z, 3a,..., etc. from G. Let y 
be the integer next to z, i.e. if z = a 1 + P-, where a 1 is 
the integer Part and p-j the fractional part less than 1, 
then y 1 = a 1 + 1. If ^ - 0, y 1 = a 1 . Let y 2 be the integer 
ne^t to 2z; i.e. if 2z = a 2 + v 2 where v 2 is the fractional 
part less than 1, then r 2 = a 2 + 1 if » 2 f 0, otherwise, 
y 2 = a 2* ^&is is continued until Y m , where Y ffl is the 
integer such that y m = mz. The points corresponding to 

V y 2» y 3’ ***’ Y m are shcy/m in Fig, 64 . Note that the 

slant ray A^ in Fig. 6s. passing through the vertical line 

n 2 = m will pass through the point y . 

m ■ 
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From Fig, 6 3, the values of € and the ass -Delated 
■values of n 2 and n-j can be written down. This is shown in 

Table 5-1. The procedure for constructing this table is as 
follows; 

£ = n 2 - n-j z 
and z > G >- 1 

Case (i) £ is Positive: 

Write down the values off , n £ for = 1 , then 
for % = 2 etc. For n^= 1, the first possible value for £ 
would ^ - z with ng = y-j * To obtain the other values, 

the value of n 2 should be increased in steps of 1 S o that 
we obtain 1 + Yj - z, 2 + r 1 - z, 3 + r-j - z, ....etc. for 
G and 1 + r 1f 2 + y 1 , 3 + Y 1 , . . . .etc. for n 2 . This is 
continued until the value of £ exceeds z. Let the maximum 
value of G that does not exceed the value of z with nq= 1 
bee . 

i.e. £ < z 

= n 2 - n-j z 
= n 2 - z 

Cn 2 , n-j are the corresponding values of 
n 2 and n-j respectively so that £ = n 2 - n-j z) 
i.e. n 2 - z <z 
n 2 < 2z 

The maximum value of n 2 that satisfies this inequ- 
ality is (y 2 - 1 ). Only the maximum value is to be taken 
since ^ is to be maximum. 


( 5 . 84 -) 

( 5 . 86 ) 
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m 

* * e= r 2 - 1 - Z 

n 2 = Y 2 “ 1 and rLj = 1 

Similarly, 

for n.] = 2 
£ — np - 2s 

n 2 - 2z ■£ z 
or n 2 C 3z 

"2 = r 3 - 1 

for n 1 = 1 

^ — Hp 1 Z S 

n 2 < (1 + 1 ) 2 
H 2= t 1+1 - 1 

Case (ii) £ is Negative 

€ = n 2 - n-j z, z <e< -1 

n 1 can not be zero ( If it is zero n 2 should be tie 
gative to make £ negative) 

if n-j = 1 a positive integer 
£ = n^ - n^ z = n 2 - 1 z > - 1 
0 > n 2 - 1 z > -1 
i.e. n 2 <1 z 
and n 2 > 1 z - 1 
i.e. Iz > ru > 1 z - 1 
t * t n 2 = y-j_ - 1 

The table for negative values of 0 is constructed 
using this equation. 
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table 5.1 



Pos 

sable 

Values r 

of € 

and 

Associated 

and n. 



Pos 

itivo 0 : 








e 

n 1 

n 2 


0 

n 1 

n 2 


Y 1 

- Z 

1 

Y 1 


Y 2 - 2z 

2 

Y 2 

1 + 

Y 1 

- z 

1 1 

+ Y 1 

1 *-+ 

y 2 - 2z 

2 

1 + y 2 

2 4 - 

Y 1 

- 2 

1 2 

+ Y 1 

2 + 

Y 2 - 2 * 

2 

2 + y 2 

3 + 

Y 1 

• 

- Z 

1 3 

+ Y 1 

• 

3 + 

Y 2 ' 2z 

# 

2 

3 4- y 2 


• 

• 


• 

• 

» 

• 
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• 

• 

• 

• 

» 


* 


• 

• 


• 

• 

• 

f 

eg 

- 1 

- z 

1 Y 2 

- 1 

Y 3 " 

1 - 2z 

2 

V 3 - 1 






0 


n 1 


n 2 


continue upto 

~ 1} 2 


m - 1 

Y c 

1 ~ 1 

n 1 = m . 

1 + Y m-1 - 1 ) 

z 

(m - 1 ) 

1 

+ Y m-1 

where m z is 

2 + -£m - 1 ) 

z 

(m - 1 ) 

2 

+ Y m-1 

an integer 

3 + Y m _i -(m - 1 ) 

z 

(m - 1 ) 

3 

+ Y m-1 

- 7™ 

m 

* 


• 


m' 


Cm -1 ) V, ' 1 


Y m-1 -1 " 2 




WAVELENGTH Vs. ABSORPTION 
FILTER NO- 2938 

MAXIMUM TRANSMISSION AT 5950 A 



WAVELENGTH 





WAVELENGTH Vs. ABSORPTION 

FILTER NO- 3214 

MAXIMUM TRANSMISSION AT 5004 


o< 






Negative £ : 


e 

n 1 


n 2 

1 

— x 

1 

£3 

1 

Y 1 

- 1 

To - 1 - 2z 

d 

2 

Y 2 

- 1 

Y 3 - 1 " 32 

3 

y 3 

- 1 

Y^ - 1 - 4-3 

4 

y 4 

- 1 

continue till 

• 


* 

Y - 1 - ms 

• 

• 


• 

• 

m 

III 

Y 

- 1 

= - 1 

m 


. .Y = nz 

. m 


Example And. Experimental Results; 

To show the dependence of the stress -optic coe- 
fficient on the wavelength, three filters with wavelengths 
5950 °a, 5520 o A and 50$+°A were used to measure the fringe 
constant. The model was a circular disc of diameter (d) 
4.73 cm. with a thickness (h) of 0.648' cm. Figs. 65(a), 

(b) and (c) show the characteristics of the three filters. 
Fig. 66. shows the plots , of load in gms. vs. fringes at 
the center of the disc. The geometry of loading is also 
shown in the figure* Table 5.2 gives the values of the 
material fringe constant F and the stress-optic coefficient 
C for each of the filters. 

At the center of the disc (p - q) = 8w/4rt d 
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TABLE 5.2 


Wavelength 

Material fringe 
constant (F) 
gm./cm. fringe 

C 

°A fringe, cm/gm 

595 0°A 

210.1 

28.28 

55 20° A 

192.0 

28.79 

5cc4°a 

172.0 

29.il 


Though the values of C donot differ very much, 
stiU neglecting them can give erroneous results as will 
be shown subsequently. For the determination of the inte- 
gral fringe orders from the measured fractional fringe 
orders, two filters with 5950°A and 5520° A wavelengths were 
used. The fractional fringe orders were measured by the 
Tardy method. Table 5»3 gives the values of the experi- 
mentally observed fringe orders (integral plus fractional 
values) and the theoretically calculated values based cn 
the model fringe constant. This table reveals that the 
fractional fringe orders measured by the Tardy method are 
quite accurate*'' 1 ■ 
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5.3 


L °ad P(gni) at 
end of lever 

59 50 °A f i It or 

Sxptal. The or. 

55 20 ° A filter 

Exptal. The or. 

575 

5.25 

5.23 

5-73 

5-73 

625 

5.70 

5-68 

6.23 

6.23 

675 

6.14- 

6 .14- 

6.73 

6.73 

725 

6.61 

6.59 

7.23 

7-22 

775 

7.06 

7.0+ 

7-73 

7.72 

825 

7.50 

7-50 

S .20 

8.22 

In Fig. 

66 w 

= % = 39- 7P 



(arb) 11.2 


= 3.54-7P 

Frcsn Table 5*2, we have 

*1 = ^950OA, F 1 = 210.1 gm/em. fringe 
*2 - 5520°a, F 2 = 192.0 gm/em. fringe 



Tables 5.4- and 
with § and z* , 


= - 1 # °5^ and z* = il - 1 .078 

5 .'5 give the values of 6 that one can expect 
These are the "ready -reckoner » tables. 


table 5.5 which 


independence of 0 on » ) will give erroneous 


( assuming 
as shown 


in Table 5*6* 
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TABL5 5.4 
(Using z) 


lil s 

£( Positive ) 

ni 

n 2 

G (Negative ) 

n 1 

n 2 

s = 1 

0.906 

1 

2 

- 0.094 

1 

.1 

2z =2.1 33 

0.812 

2 

3 

- 0.1 °8 

O 

tL 

2 

3 3 =3 

CO 

• 

C3 

3 

4 

- 0.282 

3 

3 

4z = 4 .376 

0.624 

4 

5 

- 0.376 

4 

4 

5 2 = 5 .4 70 

0.530 

5 

6 

- 0.470 

5 

5 

6z = 6.564 

0.436 

6 

7 

- 0.564 

6 

6 

72 = 7.658 

0.342 

7 

8 

- 0.658 

7 

7 

8z = 8.752 

0.248 

8 

9 

- 0.752 

8 

8 

92 = 9.846 

0.154 

9 

10 

- 0.8+6 

9 

9 

10z =10,940 

0.060 

10 

11 

- 0.940 

10 

10 


1 .060 

10 

12 



• * 


etc . 


etc 


etc 
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TABLE 5.5 
(using z* ) 


a|t 

m z £ (Positive) 

z = 1 .078 0.922 

2z = 2.156 0.3*4 

3z = 3 -234 0.766 

4z= 4.312 0.688 

5z = 5.390 0.610 

6z = 6.468 0.532 

7z = 7-546 0.454 

8z = 8.624 0.376 

9z = 9.702 0.298 

lOz =10.780 


n 2 €*( Negative) n^ n 2 


1 

2 

- 0.078 

1 1 

1 ! 

2 

3 

- 0.156 

2 

2 

"3 

4 

- 0.234 

3 

3 

4 

5 

- 0.312 

4 

4 

5 

6 

- 0.390 

5 

5 

6 

7 

- 0.468 

6 

6 

7 

8 

- 0.546 

7 

7 

8 

9 

- 0.624 

8 

8 

9 

10 

- 0.702 

9 

9 

10 

11 

- O.78O 

10 

10 


0.220 


10 10 
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taels 5.6 


load F 
(g~s. ) 

Tract . 

Ti 

Fr. Ord. 

5 

2 


e* 

z * 3 1 -»2 

Sxptal. 
n 1 n 2 

fromG from £* 
n 1 n 2 n 1 n 2 

575 

0.25 

0.73 

- 0.46 

- 0.46 

5 

ET 

J 

5 

5 

6 6 

625 

0 . 70 

0.23 

+ 0.54 

+ 0.53 

5 

6 

5 

6 

6 7 

675 

0.14 

0.73 

- O .58 

- 0.61 

6 

6 

6 

6 

8 8 

725 

0.61 

0.23 

+ 0.44 

+ 0 .4 3 

6 

7 

6 

7 

7 8 

775 

0.06 

0.73 

- 0.66 

- 0.66 

7 

7 

7 

7 

8 8 

825 

0.50 

0.22 

+ 0.33 

+ 0.33 

7 

8 

7 

8 

9 10 


For the wavelengths \ and * 2 , the fractional 
fringe orders ^ and ^ are experimentally determined. 
Using z and z* the values of e and £* are then calculated. 

FrCSa Tables and 5.5, for the values of £ and e* , the 
corresponding values of ^ and n 2 are noted. These values 
are compared with the experimentally observed values. As 
pointed out, values based on z*, i.e., neglecting the depe- 
ndence of C on X can lead to wrong values. It is more 
appropriate to use z (i.e. %/Fg ) . Also, assuming that the 
values in Table 5.4 are calculated more accurately, the 
experimentally determined values of £ may differ slightly 
as given in Table 5.6 depending on the method employed to 
determine the fractional fringe orders. 



Z4€ 


Conclusions : 

Using two wavelengths and knowing the material 
or model fringe constants for these two wavelengths, one 
can determine the integral fringe order existing at any 
p oint in a plane model by measuring the associated fractional 
fringe orders. It is shown, that the usual assumption of the 
stress-optic coefficient being independent of X may lead to 
erroneous results. For any two given wavelengths and a 
model, a 'ready -reckoner ' table can be formulated to deter- 
mine the integral fringe orders based on the measured frac- 
tional fringe orders. Alternatively , a nomogram can be 
constructed. 

In this chapter, we have discussed various methods 
of finding the integrated optical effect in a photoelastic 
model under three-dimensional state of stress graphically. 

Two new graphical methods which are derived without referr- 
ing to Poincare sphere representation are presented. One of 
these methods has the advantage that contains only straight 
lines and so it is easy to construct, in the second section 
we have presented the construction of a large field compen- 
sator using 8 retarders of equal retar dance. This is esse- 
ntially an extension of the works of Pancharat nam and Tuzi. 
Also presented are a nomogram and a ‘ready reckoner table* 
with which one can find the integral fringe orders u sing the 
fractional fringe orders for two wavelengths . This is essen- 
tially an extension of Pant's work. Where as with Pant’s 



method or the use of Babinet -Soleil compensator one can find 
the integral fringe order upto few fringes (about 6), with 
the nomogram or the ready -reckoner table one can determine 
the integral fringe order virtually upto any number. 



CHAPTER 6 


experimental RESULTS 


the previous chapters, different experimental 
techniques for the determinate of various characteristic 
parameters were discussed. These were essentially built 
on the concept of light ellipse. These experimental tech- 


niques are very general in nature. In the present chapter 
we shall present some experimental results in support of 
these techniques for completeness. In addition to these, 
we have already described some experimental results in the 
previous chapter, connected with the determination of inte- 
gral fringe orders in photoelasticity. 


Test Mode l: 

Tests were conducted on a rectangular bar under 
combined torsion and tension. Fig. 67 shows the model and 
the light path under consideration Fig. 68. is the sketch 
Of the loading device. Figs. 69. and 70. show the experi- 
mental set up and the loading arrangements, respectively. 
The model material was paraplex. The experiments were con- 
ducted with a 8" - lens type transmission polariscope shown 

in Fig. 69. Two models with the following dimensions were 
used: 
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UJ a = 1B ; b = 0.528"; c = 6}" 

(ii) a = 1 "; b = 1/2 ". c = ^ 

The light sources used uere mercury and sodium vapour 


lamps* 

5oo4°a 


Three filters of wavelengths r?50°A, 5520 °a and 
were also used with white light. The model is checked 


for residual stresses with a laser (Helium-Neon; 6328°A, 5mv). 


The scattered light fringe pattern is observed when (i) the 
laser beam is along the optical path under consideration and 


(ii) the laser beam is perpendicular to this optical path. 

No fringes were observed in case CD . The fringe pattern 
observed in ca'se (ii) is shown in Fig. 71 . This shows that 

for the optical path under consideration no residual stresses 
exist. Also the residual stresses are normal to the surface 
of the model (i.e. parallel)to the optical path. 

The objects of the -various experimental tests are as follows: 
O ) Experimental determination of the characteristic 


parameters for different loadings and checking tfe 
validity of the governing optical equations; 
Experimental determination of the characteristic 
parameters for different wavelengths of light and 
checking the validity of the governing optical 
equations; 

Tests in support of the existence of an optically 
equivalent system; 

Checking the validity of the assumptions involved 
in the previous scattered light methods^* ^ ^ *^^or 
different loadings. 
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-he details of these tests will be now presented 

(1 ) Experimental Determination of Characteristic -Parameters 
Eor various Loadings (Refer Chapter 3) 

Ex per intent 3,1 Ppo c e dur e : 

Method 1 : Load the model and keep the first quarter-wave 
Plate axes at 45° to the polarizer axis so that we have a 
circularly polarized light at entrance to the model, with 
the second quarter-wave plate out of the field, rotate the 
analyser for minimum intensity and note this direction a, . 
This direction corresponds to the minor axis of the light 
ellipse. Along the directions tt, ±45°, the amplitudes of 
the components of the light vector are equal and these dire- 
ctions correspond to the secondary characteristic directions. 

Now, introduce the second quarter-wave plate so 
that its axes are Parallel to n, and n, + 9 0° directions. 
Rotate the analyser for minimum intensity which will be 
absolute minimum. Let the minimum intensity position be Dg. 
Then a (Dg - D 1 ) = (90 + S) where S is the characteristic 
retardation. Now remove both quarter-wave plates and keep 
analyser axes along one of the secondary characteristic di- 

rect ions. R ot at e the polarii^pr fr>T* ( ~ --o ±. \ , 

jjuxj.rizer lor Cab solute) minimum inten- 
sity. The final orientation of the polarizer axis corres- 
ponds to one of the primary characteristic directions. 

Method 2: Let the optical set up consist of the polarizer 
(at arbitrary orientation), the stressed model, the quart er- 
-wave plate (at arbitrary orientation) and the analyser. 
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Rotate the analyser for minimum intensity thus finding the 
orientation of the major and minor axes of the light ellipse* 
Rotate the analyser through 4-5° in either direction and note 
the intensity with a photometer. All ellipses (of 

different orientations and ellipticities ) will have an inte- 
nsity equal to along directions which make an angle of 
+ 45° with their major axis. Now keep the angle between 
the .axes of the second quarter -wave plate and analyser 
equal to 45° and rotate the polarizer until the intensity 
from the analyser is equal to the previously found value of 
I^cj. The final orientation of the axis of the polarizer 
corresponds to one of the primary characteristic directions. 
Now remove the second quarter-wave plate and rotate the 
analyser for absolute minimum intensity. The final orienta- 
tions of the analyser axes correspond to the secondary cha- 
racteristic directions. The characteristic retardation can 
be found as in the first method or by any other method. 

The characteristic parameters for various loadings 
are presented in Table 6.1. Also presented in this table 
are the theoretical values. All the theoretical values 
presented in this chapter are found by integrating numeri- 
cally the governing optical equations in terms of Stokes 
parameters. The numerical integration is done using the 

42 

Runge-Kutte-Gill method. The series solution, (Timoshenko ), 
is used to find the values of the shear stress. The compu- 
tations were done using IBM 7C44 computer.* Also, the 
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fractional part of the characteristic retardation of the 
Phase difference at the exit along a particular set of axes 
are found directly by experiment. In most of the cases, 
the integral fringe order is one. This can be seen by coun- 
ting the fringe number at the middle point of the model frcm 
the edges where the stress is known. (Along the edges, the 
secondary shear stress is zero). By numerical Integration 
of the optical equations also, the same value for the inte- 
gral fringe order is obtained. 


We shall use the following abbreviations in this 

chapter. 

Th; Theoretical 

Ex* Experimental 

Towards right from the vertical (direction of 

tensile loading) when looking towards the source 

L: T ° Ward left from vertical when looking towards 

the source 

*ets Retardation 


Pri.ch.dir: Primary characteristic direction 
Sec .ch. dir: Secondary characteristic direction 
P: Tensile load 

^ : Torque 

Pri.str.dir: Principal-stress direction 



257 


TABLE 6.1 

Characteristic Parameters Bor Various Loadings 
Model dimensions: a = l"; b = 1/2 n ; c = 6 £ tt 
Wavelength of| light used: 5893° A 
Material fringe constant: 1 39 Psi. in/fringe 



Load 

Pri.ch.dir . 

Sec .ch.dir . 

Ch 

. ret. 

Sec .pr . 

T en sion Torsion Th. 

P T 

lb . lb . in . 

Ex. 

Th. 

Ex. 

Th. 

Ex. 

st r. dir 
at en- 
trance 
Th. 

96 

8 

37°L'~ 

31 °L 

37°R 

35 °R- 

-338° 

332° 

27°L 

96 

12 

>5+5° L 

41 °L 

45°r 

41 OR 

4010 

3940 

32°L 

96 

16 

58 6 L 

52°L- 

' 58 'ur% 2 0 R 

444° 

445° 

7 35°L 

96 

24 

1 2°L 8.5°L 

12°R 

8 °r 

470° 

45 70 , 

38ol 

90 

15 

?f°L 

50°L 

55°H 

53 °R 

436° 

433° 

35° L 

100 

18 

71 °L 

63 °L 

71 °R 

63 °R 

•F 

CD 

o 

472° 

36°L 


In Table 6.1, the agreement between theoretical and experi- 
mental values for the characteristic retardation is excellent, 
the maximum discrepancy being • In Chapter 3, we noted 
that in anti -symmetric case of stress distribution (like 
the present case), the primary and secondary characteristic 
directions will be inclined equally, but in opposite direc- 
tions to the secondary principal -stress directions at the 
entrance and exit (of the model) respectively. The values 
of the characteristic directions in the table, support 
this statement. The theoretical values for these directions 
computed numerically are symmetrical about the direction of 
tensile loading with 1 to 3 degrees difference. This may be 
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due to the effect of round off errors involved in the numer- 
ical integration. However in the tables, the average values 
of these are presented. The values of secondary principal- 
-stresses presented are calculated using the value of shear 
stress computed from the series solution. However, for some 
cases, these values are cheered using strain gages on the 
surface of the model. The theoretically computed -values for 
the secondary principal-stress directions exactly coincided 
with the experimental values. The maximum difference bet- 
ween the theoretical and experimental values of the charac- 
teristic directions is 8°. This can be attributed partly 
due to the errors involved in the numerical integration 
and partly to the possible errors in the loading alignment. 
However maximum care was taken to minimise the latter. 

The governing optical equations are thus shown to 
predict the optical phenomena in the presence of the rota- 
tion of the secondary principal-stresses fairly accurately. 
The experimental techniques for the determination of charac- 
teristic parameters are successfully employed to determine 
these parameters. From the values obtained, we note that 
difference in orientation of the characteristic and secondary 
principal-stress directions increases as the torsional load 
is increased. 

(2) Experimental Determination of Characteristic 
Parameters For Various Wavelengths of Light 
Experimental Procedure: 

The quarter -wave plates of the polarise ope are for 
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the green filter (5461 °A) attached to the camera of the pola- 
rise ope. For other wavelengths of light these wave plates 
cease to be quarter-wave plates. For green light (5461 °a) 
the characteristic parameters are found by the methods des- 
cribed earlier. For other wavelengths of light, the two 
quarter -wave plates are removed and the compensator is fixed 
just after the model. The model is taken out of the field, 
the polarizer and analyser are crossed and the compensator 
axes are rotated for minimum intensity (with some arbitrary, 
but nonzero retardation in the compensator). The axes of 
the compensator are then parallel to the axes of the polari- 
zer. After locating the axes of the compensator thus, its 
axes are kept at a certain angle to the axes of the pola- 
rizer and the retardation in it is varied until the intensity 
is minimum. The corresponding reading of the compensator 
corresponds to a zero for that particular wavelength of light. 
Few other successive readings are found in the same way which 
correspond to zero retardation or a multiple of 2tt retarda- 
tion. Thus the compensator readings are calibrated for 
that particular wavelength of light. Mow the retardation 
in the compensator is made 90° ♦ The primary and secondary 
characteristic directions are found using the method 2. The 
polarizer axes are kept at 45° to the primary characteristic 
directions. The ellipticity tan to of the exit light ellipse 
is found by Tardy's method. The characteristic retardation 
<5 = 2co . The values of characteristic parameters for diff- 
erent wavelengths of light in Table 6.2 are found experimen- 
tally in this manner. 
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TABLE 6.2 

Characteristic Parameters For Different 
Wavelengths of Light. 

Load: P = 109 lb . T = 19 lb. in. 

Dimensions: a = 1 " , b = O.jpPS", c = 6-§-'' 


Wavelenght 
of light 
used 

Material Fringe 
Constant F 
?si. in/fringe 

Pri. 

Th. 

ch. dir. Sec. 

Ex. Th. 

ch.dir 

Ex. 

. Ch. ret 
(fractio- 
nal value) 
Th. Ex. 

5950 

158 

34R 

36R 

34L 

3SL 

89° 75° 

5520 

147 

28R 

29R 

28 L 

26 L 

67° 62° 

5004 

128 

1 L 

0 

1R 

2R 

50. 5° 41° 


In Table 6.2 there is excellent agreement between 
the theoretical and experimental values of the characteris- 
tic directions. The discrepancy in the theoretical and ex- 
perimental values of the characteristic retardation is 4$. 

(if the integral fringe order of 1 is added to the fractional 
fringe orders presented in the table. The integral fringe 
order is found experimentally and by numerical integration 
to be 1 . When the wavelength of light is changed from 
5004°A to 5950°A f the orientation of the characteristic direc- 
tions changed by 34°, which is quite a marked variation. 

Thus the data (characteristic parameters) for different 
wavelengths of light can be used to obtain more information 
about the stress distribution along the optical path. The 
values in Table 6.2 a&o show that the governing optical 
equations ban predict the optical phenomena sufficiently 
accurately. It may be noted that the wavelength and the 



261 


material fringe constant donot follow any linear relation- 
ship. The material fringe constant is found experimentally 
within an experimental error +2 (i.e. about + 1.5$). 

(3) Tests in Support of the Existence of an Optically 
Equivalent System 

It has been proved theoretically, that for a parti- 
cular wavelength of light, any optical path (or a part of 
it) in a photoelastic model is equivalent to a system with 
a single retarder and a pure rotator. The characteristic 
parameters have also been defined for this optically equi- 
valent system i 

So, the following statements hold true: 

(i) When a certain amount of retardation is added along 
the primary characteristic directions at the entrance 
to the model, it appears as an additional retardation 
along the secondary characteristic directions at the 
exit of the model. 

(ii) During this process of addition of retardation, the 
amplitudes of light vector components along the secon- 
dary characteristic directions will not vary. 

(iii) Using the experimentally determined characteristic 
parameters, one can predict the exit light ellipse 
parameters for any type of incident light. 

(iv) If the characteristic retardation of the optical path 
under consideration is nullified with a compensator, 
the optical path behaves as a pure rotator. If the 
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i , 

rotatory power is nullified then the optical path 
behaves as a single retarder. 

Conversely, experimental results which confirm these 
statements form supporting evidence for the existence of an 
optically equivalent system. 

For a number of cases (of various loadings) the 
characteristic retardation is nullified with a compensator 
and then the optical path is then found to be acting as a 
Pure rotator. This confirms that the optical path is equi- 
valent to a single retarder and a pure rotator. 

Tables 6.3 and 6.4 show the variation in the inten- 
sity and retardation along the secondary characteristic dire- 
ctions, when the retardation along the primary characteristic 
directions is varied. For this the characteristic parameters 
are found for the particular loading. The axes of the first 
quarter -wave plate are kept at 45° to the primary characteris- 
tic directions. If the angle between the axes of the polari- 
zer and first quarter-wave plate is e , the retardation along 
the primary characteristic directions at .the entrance to the 
model is 2 0 . The analyser is kept with its axes along the 
secondary characteristic directions. The polari 2 er is rota- 
ted so that the retardation along the primary characteristic 
directions is varied and the variation of the photometer 
reading is noted. To find the retardation along the secon- 
dary characteristic directions for a particular value of 2 ©, 
the axes of the second quarter -wave plate are kept at 45° 
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to the secondary characteristic directions so that we obtain 
a line ellipse before the analyser. The analyser is then 
rotated for minimum intensity. If the angle between the 
axes of the analyser and second quarter -wave plate is w , 
the retardation along the secondary characteristic direc- 
tions is 2o5 . Thus the a values are f ound experimentally 
for various values of e ♦ 

In Table 6.5 the parameters of light ellipse at 
the exit for various types of incident light are presented. 
Two sets of values are presented in this table. One set 
corresponds to the values found directly by experiment. The 
second set of values is computed from the characteristic 
parameters which are found experimentally. The light ellipse 
parameters that are tabulated are the phase difference and 
amplitude ratio along the secondary principal-stress direc- 
tions at the exit of the model. These directions are com- 
puted using the series solution for the shear stress; To 
find these parameters experimentally, the following procedure 
is adopted. 

Remove the second quarter-wave plate and keep the 
analyser axes at some inclination to the secondary principal- 
-stress axes at the exit. Introduce a compensator between 
the model and the analyser with the axes of the compensator 
along the secondary principal-stress axes. Vary the retar- 
dation <$ c of the compensator until the intensity of light 
from the analyser becomes a minimum. The phase difference 
along the secondary principal-stress axes at the exit is 
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then - ^ c * Now, we have plane polarized light after the 

compensator, the direction of which can be found by rotating 
the analyser for minimum (zero) intensity. If the angle 
between the axes of the secondary principal-stresses at the 
exit and the analyser is « , the amplitude ratio along the 
secondary principal-stress axes is tan a . The proof for this 
is as follows: 


Let the Stokes vector before and after the compensa- 
tor with respect to its axes be given by V-j and V 2 which 


are: 







S 1 


S 1 

V 1 £ 

s 2 

and = 

S 9 Cos 8 - S_ Sin q 


_ s 3. 


S 3 Coss c + S 2 Sin 5 c 


when V 2 is referred to the axes of the analyser it becomes 


S-j Cos 2 © + Sin 20 (S 2 Cos <5 Q - S 3 Sin s c ) 


— —r 
Si 

(S 2 Cos <5 C - S 3 Sin 6 C ) Cos 20 - S-| Sin 2e 

= 

s 2 

S 3 Cos § c + S 2 Sin 5 C 

: 

h_ 


where 6 is the angle between the axes of the compen- 
sator and the analyser. 


96, 


= Sin 29 (- S 0 Sin 5 - S~ Cos ) = 0 

2 C 3 c 

- S„ 


if Sin 29 = 0 or tan 8 = — 2 = tan£ i.e, 

s 2 


_ e 


where £ is the phase difference along the secondary 
principal-stresses of the light ellipse from the model. 

Also Sin 20 f 0. 
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TABLE 6.3 

Variation of Intensity and Retardation Along the 
Secondary Characteristic Directions When the R etar - 
elation Along the Prima ry Char acter i stic Directions 
is Varied 

Dimensions: a = 1 ", b = 0.528", c = 6f" 

Load: P = 109 lb. T = 19 lb. in. 

Wavelength of light used: 5461 °A 

Material fringe constant: 14-7 Psi. in. /fringe 

Characteristic Parameters: 



Ex. 

Th. 

Pri. Ch. Dir. 

63° L 

62°L 

Sec. ch . dir . 

60°R 

62°R 

ch. ret. (fractional value ) % 56° 

67° 


The intensity variation along one of the secondary charac- 
teristic directions when the retardation along the primary 
characteristic directions is varied is 2 divisions on 100 


scale of photometer. 


Ret .along 

Ret. observed difference 

Ret. along 

Ret.obser- 

diff- 

pri. ch. dir . 

along sec. 


pri. ch.dir 

. ved sec. 

ere- 

(deg. ) 

ch.dir . 

(deg. ) 

( deg . ) 

ch.dir. 

nee 


(deg.) 

( deg . ) 

( deg . ) 






24 

0 

56 

18 

100 

162 

18 

20 

74 

20 

120 

184 

18 

40 

$4 

21 

140 

202 

18 

60 

115 

23 

160 

220 

20 

80 

138 

180 

240 
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TABLE 6.4 

Variation of Intensity and Retardation Along 
the Secondary Characteristic Directions When 
the Retardation Along the Pr imary Characteristic 


Directions is Varied 

^ it 


Dimensions: a = 1 , b = 0.528 , c = 

Material fringe constant: 158 Psi. in. /fringe 
Wavelength of light used: 5893°A 

Load: P = 96 lb . T = 16 lb. in. 
Characteristic Parameters: 

Ex. Th. 

Pri.ch.dir. 38.5°R 38.5°R 

Sec. ch. dir. 36°L 38.5°R 

Ch. ret .(fractional value) -44° -44° 


Ret. along Ret. along diffe- Ret. added Ret. observed diffe- 
pri.ch. sec. ch. dir. rence along pri. along sec. rence 



With analyser axis along one of the secondary principal - 
-stress directions at the exit, the intensity reading on 
100 scale of photometer varied from 21 to 23, •when retarda- 
tion along the primary characteristic directions varied* 
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TABLE 6.? 

P arameters, of Lleht ^ £g ^ £ 

t he Optically g a ujvalent Modal 
Load: P = 100 lb. T = 18 lb. in. 

Wavelength of light used: 5461 ° A 
Dimensions: a = 1 " , b = 0.528*' , e = ^ 

Material fringe constant: 147 Psi. in. /fringe 
Ckaracteri.qi-.io Painters; 


Pri, eh. dir. 

Sec. ch. dir. 

Ch. ret. 

Sec. pr.str.dir. ) 
at entrance ) 

Sec. pr.str.dir. ) 
at exit . ) 


- 116 ° 

65°L 

65°R 

35°L 

35°R 


- 100 ° 
61 .5°L 
62°R 


Ret .at Ret Talcns sec ~ — 7_-: 

entrance pr.str.dfr.at & 

along sec. exit p * 

fde| t )* dlr ‘ , calculated Ex. 

(.deg.; (deg.) from ch.par 

(deg. ) 


Amplitude ratio along sec • 
pr.str. dir. at exit 


calculated from 
ch.par. 


1 .54 
2.05 
2.4l 

3.08 

4.01 

4.01 

2.75 

2.14 

1.73 

1 .28 
1.00 
0.78 
0.63 


1.50 

1.93 

2.54 

3.34 
4.00 
3.78 

2.97 

2.25 

1.73 

1.35 
1 .07 
0.85 
0.67 
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Prom Tables 6.3 and 6.4, we note that the intensity 
variation along the secondary characteristic directions when 
the retardation along the primary characteristic directions 
is varied is 2 divisions on the 100 scale of photometer. 
Since the photometer pointer by itself shows a fluctuation 
of 2 divisions on 100 scale, this can be taken to mean that 
the amplitudes along the secondary characteristic directions 
remain constant. When the retardation along the primary 
characteristic directions is varied in steps of 20°, the 
retardation along the secondary characteristic directions 
varied in steps of 18° to 24°. The error involved in the 
measurement of phase difference is 5° in the present experi- 
mentation. i.e. the experimental error involved can be as 
high as 5° with the present experimental set up. So, what- 
ever retardation is added along the primary characteristic 
direction the same appears as an additional retardation 
along the secondary characteristic directions. In Table 6.5 
the two sets of values of the phase difference € and the 
amplitude ratio e_ are compared. One set of values is 

-X. 

obtained directiy by experiment afld the other set of values 
is calciilat ed from the (experimentally determined) charac- 
teristic parameters. The agreement between the two sets 
of values is good. Thus all these experimental results 
confirm the statements made earlier regarding the optically 
equivalent system. These form a supporting evidence for 
the existence of an optically equivalent . system for an 



optical path in a three-dimensional photoelastic system for 
a particular wavelength of light. 

CO Experimental Check for the Validity of Assumptions 
Involved in the Previous Scattered Light Methods: 

We note that the assumptions involved in the pre- 

10 11 12 

vious scattered light methods ’ * are satisfied when the 

secondary principal-stress axes at the entrance and exit 
coincide with the primary and secondary characteristic? 
directions . 

Tables 6. 6 - 6.10 show the amplitude ratio e and 
the phase difference € along the secondary principal-stress 
axes at the exit when the retardation is varied along the 
secondary principal-stress axes at the entrance. Also shown 
are the -values of e„ and £ which are computed from the 

-X. 

(experimentally determined) characteristic parameters. These 
sets of values are presented for 5 loadings. 

To vary the phase difference along the secondary 
principal-stress axes at the entrance the axes of the first 
quarter -wave plate are kept at k-5° to these secondary 
principal-stress axes. The polarizer is then rotated to 
vary the retardation along these axes. If the angle between 
the axes of the polarizer and the first quarter-wave plate 
is , the retardation introduced along the secondary prin- 
cipal-stresses at the entrance is 2 <*>. The azimuth of the 
light ellipse at the exit is found with the analyser. The 
axes of the second quarter-wave plate are kept parallel to 
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the axes of this light ellipse so that it reduces to a line 
ellipse. The orientation of this line ellipse is found with 
the analyser. From this orientation the ellipticity and the 
azimuth of the light ellipse from the model are found, The 
values of e x and £ are calculated from the parameters of 
this light ellipse. The characteristic parameters are found 
experimentally by the methods described earlier. The values 
of £ and e^ for various types of incident light are also 
computed using these parameters. Sodium vapour lamp (5893°A) 
was used as the source of light. Only the quarter -wave 
plates of the polariscope were used for the measurement of 
ellipticity etc. These introduce only 85° of retardation 
for wavelength 5893° A and not 90°. 

Wavelength of light used: 5893°A 

Material fringe constant: 139 Psi. in. /fringe 

— i) 

In the tables to foolow, the following abbreviations are 

: Retardation added along secondary 
principal-stress axes at entrance. 

: Intensity of light along on§ of the 

secondary principal -stress axes at exit 
as given by the reading on 1000 scale 
of photometer. 

: Azimuth of the out coming light ellipse 
with respect to the vertical (line of 
tensile loading) 

Q.w.p. : Second quarter-wave plate position 


used: 

Ret. ent. 

Intensity 


Initial azimuth 
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Final azimuth: 


Ph.diff . 
(Ex.) 


F,i. diff. 
(Calc. ) 


Ellipticity 

(Ex) 


Ellipticity 
(Calc . ) 


Ampl. ratio 
(Ex. ) 


Ampl. ratio 
( Calc . ) 


Final azimuth of the outcoming line ellipse 
after the introduction of the second quart er- 
-wave plate. 

Phase difference along the secondary princi- 
pal - stress axes at exit calculated from 
the exit light ellipse parameters which are 
found ex perimen tally. 

Phase difference along the secondary princi- 
pal-stress axes at exit calculated from 
characteristic parameters which are found 
experiment ally. 

Ellipticity of the outcoming light ellipse 
found directly by experiment. 

Ellipticity of the outcoming light ellipse 
calculated from the characteristic parameters 
which are found experimentally. 

Amplitude ratio along the secondary princi- 
pal-stress directions calculated from the 
exit light ellipse parameters which are 
found experimentally. 

Amplitude ratio along the secondary princi- 
pal-stress directions calculated from the 
characteristic parameters which are found 
experimentally. 
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TABLE 6.6 


? = 90 lb. T = 15' lb* in. 


Ret. ent . 

Initial azimuth 

* •w* p 

Final azimuth 

0 

a 4 

0 

Jr 

50° L 

7°R 

1+5° 

72 °r 

27°R 

46 °R 

60 ° 

91 °R 

46°R 

72°R 

90° 

79° L 

56°r 

90°L 

1 20 ° 

76°L 

59°R 

72°L 

135° 

71°L 

64 °R 

58°L 

150° 

65°L 

70°R 

44°L 

1800 

55°L 

80 °r 

23 °R 


TABLE 6.7 

P = 96 lb • T = 1 2 

lb . in . 


Ret . ent • 

Initial azimuth. 

Q.w.p, 

Final azimuth 

0 

76 °L 

59°R 

520L 

30° 

67°L 

68 °R 

35°L 

1+50 

55 °L 

8 o°r 

17°L 

O 

O 

vO 

40°L 

85°r 

1 °L 

90 P 

93°R 

48°R 

71 °L 

120 ° 

3°L j 

48°L 

11 °h 

t35° 

3°R 

42°L 

10 °fl 

150° 

2°R 

44° L 

4°L 

180 ° 

5°R 

4o°L 

8 °L 



TABLE 6.8 



P = 96 lb. T = 

8 lb. in. 


Ret. ent. 

Initial azimuth 

Q.w.p. 

Final azimuth 

0° 

85 °R 

40°R 

71 °R 

30° 

85 °R 

40 OR 

95° L 

45° 

86°R 

41 °R 

8?°L 

60° 

83 °R 

38°R 

81 °L 

90° 

84 °R 

39° l 

63°L 

120° 

1 2°R 

33°’L 

56°r 

135° 

9°L 

54°l 

S°R 

150 ° 

8°L 

53°t 

20°R 

180° 

7°L 

52°L 

4°r 


TABLE 6,9 


P = 96 lb. T = 16 lb. in.. 


Ret.ent. Initial azimuth Q.w.p* Final azimuth 


0 

38°L 

83°L 

1 0°R 

30° 

17°L 

62° L 

5°R 

O 

O 

\D 

5°L 

50°L 

f°R 

90° 

4°R . 

4l°L 

4°R 

120° 

16°R 

29°L 

3°R 

135° 

1 7°R 

28 °L 

0 

150 0 

28°R 

1 8°L 

3°R 

180° 

56°R 

11°R 

28°R 



•TABLE 6.10 


2 75 


p = 96 Lb. T = 2k lb. in. 


Aet . ent • 

Initial azimuth 

Q.w.p. Pi 

nal azimuth 

0 

4-Or 


4-1 °L 

12°L 

30° 

22°R 


24-0 l 

5°R 

4-5° 

28 °R 


+7°L 

14- or 

60° 

330R 


7°L 

27°R 

90° 

4-8°r 


3°R 

4-5°r 

1200 

6 6 °r 


21 or 

73°R 

1500 

?2°R 


27°R 

84- °r 

180° 

84-°l 


51 °R 

65° L 


1 - — 



... .... _ 



TABLE 6.11 



W £-\ 4 - s-w*. 4 * 

P = 90 

lb. T 

= 1 5 lb. in. 


•ii 0u • ent • 

Ph 

Ex. 

. <Stff . 

Calc. 

Elliptic it y 

TSV * rt « -t ^ 

0 

77° 

84-0 

0 .65 

0.53 

4-50 

•660 

52° * 

0.4-9 

0.4-9 

60° 

199 

•3° 

0.34- 

0.38 

90° 

22° 

21° 

0.I9 

0.15 

120° 

8° 

14-° 

+O.07 

+0.07 

135° 

28° 

38° 

+0.24- 

+0.20 

1500 

4-60 

620 

+0 #38 

+0.31 

1800 

73° 

84-0 

4*0 * 63 

40*53 
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TABLE 6 .12 

P = 96 lb. T = 12 lb. in. 


Ret. ent. ' 

Hi' 

Ex. 

.diff. 
cal* ... 

S 

Ex. 

llipticity 

0 

1+8° 

39° 

0.4-5 

vdlC * 

0.29 

300 

650 

71° 

0.63 

0.59 

4-5° 

800 

85° 

0.78 

0.76 

600 

87° 

32° 

0.81 

0.84- 

90° 

34-o 

57° 

0.29 

0.53 

120° 

-32° 

+29° 

0.25 

0.24- 

135° 

15° 

+1i+o • 

0.12 

0.11 

1500 

110 

3° 

0.09 

0.02 

130° 

27° 

39° 

0.23 

0.29 


T M: uS 6.1 4 

P = 96 lb. t = 8 

lb. in. 


Ret. ent- 

Ph. 1 

lip f jr 

Jii /v m 

diff. 

Calc . 

Elliptic it 7 

0 

37° 

30° 

0.25 

0.25 

30° 

20 

20 

0.00 

0.02 

1+50 

22° 

19° 

0.14- 

0.15 

60° 

4-6° 

34-° 

0.30 

0.29 

900 

7+° 

63° 

0.66 

0.60 

1 20° 

880 .. 

89° 

0.97 

0.93 

135° 

74-0 

75° 

0.64- 

0.73 

1500 

67° 

61° 

0.53 

0.55 

180° 

31° 

300 

O.19 

0.25 



T^BLE 6.14 
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■p> t ""L »> — 

P = 

96 lb. 1 

Bet. ent. 


•dif. f. 
Calc, 



0 

86 ° 

73° 

30° 

5 9° 

7+° 

60 ° 

27° 

30° 

O 

O 

O 

1 ° 

6 ° 

120 ° 

28° 

30° 

135° 

36° 

O 

O 

150 ° 

5 6 ° 

50 ° 

180 ° 

1 

» 

00 

MO 

O 

73° 


TABLE 6 .I5 


16 lb. in. 


_ Ellipticity 
■M*. Calc. 

0.54 

0.53 

0.4o 

0.39 

0.21 

0.17 

0.00 

0.04 

0.24 

0.26 

0.32 

0.37 

0.47 

0.46 

0.53 

0.53 






p = 

96 lb . T 

= 24 lb . in . 


xt Go # Gnt » 

Pn 

Ex. 

• diff. ~ 
Calc • 

Elliptic it y 

_ -1 _ 

0 

32° 

3^° 

O .29 

uaj,n 

0.28 

30° 

37° 

52° 

0.30 

0.27 

45 0 

36° 

70° 

O .25 

0.24 

60 ° 

4l 0 

73° 

O .19 

0.19 

900 

4l 0 

130 

0.06 

0.06 

120 ° 

31° 

10 ° 

0.13 

0.08 

•Jh 

vn 

0 

0 

35° 

23° 

0.21 

0.20 



TABLE 6.16 


2?8 


Bet . ent . 


0 

450 

60 ° 

90° 

120 ° 

135° 

150° 

1800 


P - 90 lb . T = lb • in . 


Intensity 

Ampl. " 
Ex. 

ratio. 

Calc. 

23 

1 .44 

1 .87 

13 

1 .76 

O.99 

11 

1 .28 

0.81 

6 

0.97 

0.53 

4 

O.87 

0.37 

3 

0.75 

0.34 

5 

0.68 

0.37 

8 

0.71 

O.54 

TABLE 6.17 

96 lb. T = 

12 lb, in. 



Pet » ent * 
0 

30° 

^50 

60 ° 

90° 

120° 

135° 

150° 

I800 


Intensity’ 
8 

10 
12 
14 
■ 20 
23 

23 
25 

2-3 


Ampl. ratio. 
-M* Calc. 


0.98 

0.86 

0.84 

0.82 

1.35 

0.60 


h .1 

1 .49 
1 .30 
’ 1.13 

rggg 

0*64 


0*71 0.58 
0.67 O.54 
0.78 0,55 
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TABLE 6.18 



P = 96 lb. T 

= 8 lb. in. 


Ret . ent . 

Int en sit y Am pi . 

Ex. 

ratio 

Calc . 

0 

11 

2.11 

1 .48 

30° 

11 

2.53 

1 .51 

1+5° 

11 

2.29 

1 .45 

60° 

13 

2.18 

1.36 

90° 

17 

1 .36 

1 .14 

120° 

24 

0.99 

0.93 

135° 

2? 

0.70 

0.8+ 

350° 

29 

0.62 

0.77 

180° 

32 

0.4l 

0.67 


TABLE 6.19 



P = 96 lb. T 

= 16 lb. 

in • 

Ret .ent . 

Intensity 

Atnpl. 
_ Es • 

ratio 

Calc. 

0 

8 

0.55 

0.57 

30° 

14 

0.51 

0.41 

60° 

20 

0,61 

0.39 

90° 

25 

0.81 

o.53\ 

120° 

28 

/ 1 .20 

' | /■ 

0.79 \. 

135° 

27 

j ( 1.22 

0.96 ! 

150° 

27 

1.47 

1.17 : 

18o° 

23 

1. 38 

1.^+/ ■ 
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TABLE 6.20 




P = 96 lb . T = 24 

lb . in. 


Ret. ent. 

Intensity 

Ampl. ratio 

Ex. Calc. 

0 

26 

O .92 

1.53 

30° 

32 

1 .54 

2 .% 

450 

34 

1 .97 

3.9^ 

60 ° 

35 

3.17 

5.23 

90° 

33 

10.76 

3.40 

120 ° 

28 

3.65 

1 .80 

150 ° 

19 

2.38 

1 .08 

180° 

.12 

1.03 

0.65 

In 

Tables 6.6 - 6.10, the 

experimental values of 

the azimuth 

of the light ellipse at 

the exit 

and its orien- 


tation after it is reduced to a line ellipse are presented* 
From these values, the ellipticity, the amplitude ratio and 
phase difference along the secondary principal-stress direc- 
tions at the exit are calculated. In Tables 6.11 -6.15, 
these values of amplitude ratio and phase difference are 
compared with the values calculated from the characteristic 
parameters (which are found experimentally). In Tables 6.16 

■ ; ,vS' 

6.20, the intensity along the secondary principal-stress 
directions (as given by the photometer reading) arid the amp- 
litude ratio along the secondary principal-stress directions 
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at the exit-(£ ^calculated from the light ellipse para- 
meters at the exit and (ii) calculated from characteristic 
parameters.,,., are presented. Fran these tables we note that 
both the sets of values of the ellipticities match very veil, 
whereas for phase difference and amplitude ratio, the match-*, 
ing is not good. The reasons for the discrepancy in the later 
case are; (i) the wave plates used introduce only 85° of re- 
tardation for wavelength of light 5893^A (ii) errors in 
measurement of the orientation of light ellipse severely 
afiect the values obtained for phase difference and amplitude 
ratio since the constant elliptic it y curves are very steep 
for azimuth values near 0 and 90°, Fig. Lpf • (iii) errors 
in measurement of orientation of light ellipse do not have 
a significant effect on the values of ellipticity obtained^ 

Also, we note from these tables that the assump- 
tions of the previous scattered light method! 0 * 11 are not 
valid for the loadings considered. When the retardation 
along the sec. pr. str. directions at the entrance is varied, « 
at- -. the exit, the intensity along the secondary principals 
-stress directions does vary and the phase difference along 
these directions does not increase by the same amount as the . 
retardation added at the entrance. 

In this chapter, we have presented experimental 
evidence for a number of cases in support of the. following 
stat ement ss 



The total light path in a three-dimensional photo- 
elastic model for a particular wavelength of light 
is equivalent to a system with a retarder and a pure 
rotator. The three parameters of this equivalent 
system called the characteristic parameters can he 
found by the techniques described earlier. 

The governing optical equations (derived in Chapter 2) 
describe, the optical phenomena in three-dimensional 
cases where rotation of the secondary principal-stress 
axes along the optical path is present. 

The variation of the characteristic parameters with 
wavelength of light can be quite marked in some cases. 
In the case considered, the characteristic directions 
varied by 3*+-° when the wavelength of light is varied 
by 9 1+ 6°A.. 

The assumptions involved in the previous scattered 
light methods are not valid when the rotation of 
the secondary axes is large. These assumptions 
lead to large errors in such cases. 



CHAPTER 7 


COLIC LOS ION 

This thesis which consists of six chapters in 
addition to the present one deals ’- T ith the problems of 
stress -induced optical activity and determination of the 
optically equivalent model in three-dimensional photo- 
elasticity. Chapter 1 , the introduction, describes the 
general problem when the secondary principal-stress axes 
and their magnitude continuously vary along the light path 
in a stressed photoelastic model. It gives a brief review 
of the stress-optic laws in two - and three-dimensions and 
attempts to modify these in the presence of optical activity 
It also discusses briefly some of the scattered light tech- 
niques with specific relevance to the effects of optical 
activity. Finally, the objectives of the present investiga- 
tion are stated along the methods of approach. 

In Chapter 2, the concept of light ellipse is 
introduced. It is shown that studying the properties and 
transformations of light ellipse and then using this infor- 
mation in photoelastic analysis has a definite advantage 
over the usual procedure of resolving the light vector into 
its components trigonometrically along various directions. 



ox' f “ 13 SlmPler aGd P^ents the Physical aspect 

tnS Pr ® 1 « “ed readily. The u , ht * ^ 

«">* in terms of Sto v, s / ^ *“«•• ls ^ 

^Oa, s sector and throughout th-?« • 
sation the concept of lleh x M1 . “ lnVesti - 

bsen used. 1PSG 3 " d St3kes V90t °r Have 

». „,C«z to ” *™.. « 

tion of secondary principal-stress ^ PPeSenCe of tha rota- 

t hat f ho , . SS axes is presented to show 

hat the kinematically derive . h 

y derived equations nrorHof 

Phenomena sufficiently P ^ * t optical 

lc iently accurately. Thpn«t* 

are thGn “ kinematically m teims of 

t tura out to be a - - --I: 

linear different-? -si « . . irst order 

These ^nations with variable coefficients 

These are very convenient for analvsx* a 

ITsi ncr XV, aalysis and application. 

lng these equations the analvMr»»i 

naiyt ical solution for 

where ratio H of the rat, of „ , CaS6 

arbitrarv f x- v " ^ 0 (thls ° an ie 

7 unction) to the rate of retardation }' • 

tant. The case ax ls a cons_ 

, hiacussed by Bructer and Mindlin where 

both , and r are constants is a particular case of this. 

lence of ^t** 16 46441163 Pr °° f the eqUiva - 

to a s ^ stom containing a 

e ar der and a pure rQtator ^ presentea< ^ ^ Polnoare , 

P ere representation of Polarised f orms for thls proof ±a 
also demonstrated. The proof is then extended to any i ength 

0P 1Cdl ln a three - dim ^sional Photoelastic model 
that is, any len g th of opt icai path in a three-dimensional 
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photoelastic model is again equivalent to a system contain- 
ing a single retarder and a pure rotator. A discussion of 
some special cases (case of a pure rotator, anti -symmetric 
case, axi -symmetric case, reversal of optical path, reflec- 
tion polarise ope) in which there exists a relation between 
the three characteristic parameters, is presented. 

Two simple and direct methods of determining experi 
mentally the characteristics parameters are presented. It 
is to be noted that there donot exist any such methods to 
find the characteristic parameters directly in a simple my. 
However Kent and Lawson ^ suggested an ingeneous device to 
measure the ellipticity directy. This has been used by 

14,15,1$ 

Robert and Guillemet to find the characteristic para- 

meters. But this method requires the rotation of the analy- 
ser at a particular speed and some additional instrumenta- 
tion. 

In Chapter 4-, a review of some of the previous 
scattered light techniques, the assumptions involved and the 
analysis, is made. A critical analysis of these methods is 
made making use of the analytical solution (obtained in 
Chapter 2) for the case where R = J0'/p» is a constant <0' 
can be any arbitrary function) and numerical solutions (of 
the governing optical equations) for some cases. The limi- 
tations of the previous methods are discussed. 

The scattered light techniques essentially involve 
the following steps: (1) measurement of scattered light 
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ntensity along different direction taking into account, 
he variation of the lengths of light path involved, (2) 
finding the parameters of the light ellipse and the oharact- 
ameters at the point under consideration (3) e va- 
-ting the functions fi and ? using these values and the 
g ernmg optical equations. Some new scattered light tech- 

are presented which donot involve any assumptions or 

he rotation of the analyser at a Dart • i 

• y aT: d Particular speed. 

Chapter 6 deal'' iru-v, o 

" lth 3 Problems, (i) graphical 
methods in Phot oelasticitv +-v, 

, , s.icity, (ii) theory of compensation and 

(in) determination of integral fringe order in photoelas- 

ticity. Graphical methods are Very convenient to find the 
integrated optical effect in the presence of the rotation 
of the secondary principal-stress axes. A review of the 
earlier graphical methods is made. Of these the stereo- 
graphic and parallel projections are derived from Poincare' 
sphere representation of polarised forms. The methods using 
quaternions needs geometrical construction for each problem. 

In our discussion a new approach to the graphical methods 
IS described which doesnot involve the use of Poincare' 

Ihe light ellipse has essentially four parameters 
the phase difference, elUotieity, amplitude ratio and 
azimuth. These are connected by specific relations which 
can be used to plot two of these four Parameters keeping a 
third parameter constant. Based on these, two graphical 
methods are presented. One of these is very useful in find- 
ing the effect of errors in various polariscopes. 
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The second method Is very simple to construct since it 
involves only straight lines. Etocedures are described to 
determine the integrated optical effect and the characte- 
ristic parameters using these graphical methods. 

The second problem deals with the theory of com- 
pensation. nse Of (a) a systeffi of four ^ 

mg the same retardation with a polarizer or an analyser 
and (b) a system of eight retarders each having the same 
retardation, as a large field compensator is described, 
ccncept of an optically equivalent system is made use 
in deriving this. The angles between the axes of diff- 
erent retarders in either system (a) or in (b) are to be 
adjusted such that the rotatory power of the equivalent 

system for this system of rpfaHono v „ 

^ 01 retarders becomes zero. The 

principle involved is essent ini i-o- a-n 0-0-4- • „ 

essentially an extension of the work 

uy Pancharatnam 28 and Tuzi 2 ? . 


The third problem deals with the determination of 
integral fringe orders in two-dimensional photoelastic 
models under load. The drawbacks of the present methods 
are discussed. Construction of a nomogram and a ready 
reckoner table with which one can determine the integral 
fringe oraers virtually upto any number. This procedure 
doesnot have the limitations of the previous methods. This 
is essentially an extension of Pant 1 s work in which the 

integral fringe orders are determined from the fractional 
fringe orders for two wavelengths of light when the two 
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total fringe orders donot differ by more than unity. 

Experimental results are presented in Chapter 6 
as supporting evidence to the theory developed in the pre- 
vious chapters. The experiments were conducted with a rec- 
tangular model under combined tension and torsion in a 
transmission polariscope . These form the supporting evi- 
dence for the following statements: 

(1 ) The governing optical equations in terms of Stokes 

parameters derived in Chapter 2 predict the optical 
phenomena in the presence of rotation of secondary 
principal-stress axes (along the optical path) 
sufficiently accurately. 

(2) The techniques described in Chapter 3 can be 
successfully used to determine the characteristic 
parameters . 

(3) The change in the characteristic directions is 

% ' 

found to be 3 1 )- 0 when the wavelength is changed 
from 50d>+ o A to 5950°A (the corresponding change 
in the material fringe constant is from 128 -158 
Psi. in./fringe). This is quite a significant 
variation and in cases such as this, the data for 
different wavelengths of light can be used to 
obtain more information about the state of stress 
along the optical path. 
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Any part of optical path in a three-dimensional 
Photoelastic model is equivalent to a system 
containing a retarder and a pure rotator. 

The assumptions on which the existing scattered 

light techniques are based, are not valid in a 
general case. 




APPENDIX A 


(i) Determination of the Characteristic Parameters 
When the Direction of Observation is Fixed 

From Eqs.- 0+.4-2) 

S 1 = S 1 Cos 2i i + (S 2 Cos p^- S 3 Sinpv, ) Sin 2* 

If = Cos 26, S 2 = Sin 26 CosC , S 3 = Sin 26 Sine 

S^ - Cos 2© Cos 24> + Sin 2ij» (Sin 2© CosC C°s P-^ 

- Sin 2© Sine Sin P 1 ) 

d Sj| _ sin 24* (Sin 2 e Cos Sin€ + Sin 2© Cos£ Sin P-| ) 

d e 

The above expression is equal to zero when 

tan€ ~ - t an p ^ 
or e = - P-j 

This means that if the retardation e is varied 
along the primary characteristic direction, the intensity 
along any arbitrary direction at the exit of the optically 
equivalent system becomes extremum when e = - p 1 where p 1 
is the characteristic rctardati.cn of the equivalent optical 
system. Let the angle between the axes of the polarizer n 
the quarter-wave plate be « . Mote the intensity I, • Now 
rotate the polarizer through 90° and note the intensity I 2 < 
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Change the angle between the axes of the polarizer' and 
quart er -wave plate (keeping the quarter -wave plate fixed) 
and ' repeat this procedure until l 1 = i 2 * i n this position, 
the major axis of the light ellipse at the point under consi- 
deration makes an angle of 4-5° with the direction of observa- 
tion. (i) One possibility is that this light is circularly 
polarized. If it is so, the amplitudes along the secon dar y 
characteristic directions are equal and the amplitudes along 
the primary characteristic directions are equal. But the 
light ellipse after the quarter -wave plate will have its 
major axis along the quarter-wave plate axes so that along 
directions making 4-5° with these axes, the amplitudes will 
be equal. That is, the quarter-wave plate axes make 4-5° with 
the primary characteristic directions. The phase angle along 
these primary characteristic directions after the quarter- 
-wave plate, is 2* . So, we have 2°( + p -| = ir/2. All 
these statements will be true only if the light at the point 
under consideration is circulirly polarized, (ii) The second 
possibility is ‘that the direction of observation is along 
■one of the secondary characteristic directions, (iii) Neither 
of the cases (i) and (ii) is true. 

To check these, rotate the polarizer and note the 
intensity I. (i) If I doesnot change when the polarizer is 
rotated this falls into the second case so that the direc- 
tion jf observation corresponds to one of the secondary cha- 
racteristic directions. The axes of the quarter -wave plate 
are at 4-5° to the primary characteristic directions. 
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The characteristic retardation cannot be found unless the 
direction of observation is changed, (ii) If I changes when 
the polarizer is rotated, find that position of the polari- 
zer for which the condition 1^ = is satisfied. In this 
position let the angle between the axes of the polarizer 
and the quarter-wave plate be *=Cj . Rotate the polarizer for 
minimum intensity and lot this angle be « 2 . When the angle 
between the axes of the polarizer and quart er -wave plate is 
^ for obtaining 1^ = I^ and if the light is circularly po- 
larized at the point under consideration 2*-, + p.j = V2. 

When the polarizer is rotated for minimum intensity 2*2 + P-] 

TT 

= 0 so that (=*-1 -<*2 = jj-* Change the orientation of the 

quart or -wate plate until this condition is satisfied. When 
this is satisfied, the axes of the quarter-wave plate coincide 
with the primary characteristic directions and ^ = ir/. 2 - 2 . 

= - 2^2* To find the rotatory power, nullify the characteris- 
tic retardation with a compensator and rotate the polarizer 
for minimum intensity. The angle between the axes of the 
polarizer and the direction of observation will give the ro- 
tatory power. 
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and = - SgP’Cc^^ - S-j p* San20 ( 2.44 ) 


where 0 is the orientation of the secondary 
principal-stress directions with respect to the fixed axes. 
From Eqs. ( 2*4-2 ) and (2 A3) 

S.J = 0 

when = 0 or p « = 0 or Sin20 = 0 (i.e. 0 = 0 or 7C/2) 

Also 


8.. = Cos2o: Cos2 Dt 


Similarly, 


when 


S' = 0 

S 0 — 0 orp* = 0 or Cos20 = 0 (i.ei 0 = ^"A or 3 /6 A 


and 

S 2 = Cos2w Sin2 0( 

S becomes S n when o(is replaced by ( ^A) 

1 tL 

(when w remains constant) 

. Si becomes S ’ when U is replaced by ( 

• # T 

(when (ii remains constant) 

For all points where the intensity is extremum along 
the Z-axio (i.e. along the optical path) for a particular 
direction of observation 


dS 1 

dT 


= 0 


I | 

If at a point on the optical path in the model, Sg , 

either p*= 0 or S 3 * 0 or both p* and S 3 are zero.^ To eli- 
minate the case where S = 0, the incident light is to be 
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modified. If = 0 for any type of incident light we 

can conclude that at the point under consideration p* = 0. 

At this point the directions of the secondary principal-stress 
axes can not he found since no retardation is "being added at 
this point. If = 0, hut f 0, the direction of obser- 
vation is along one of the secondary principal-stress axes. 
(Because j- 0, p’ j. 0 otherwise S-j' = = 0. But 0 = 0). 

This suggests a method of finding the directions of the secon- 
dary principal-stresses directly from the scattered light 
fringe pattern. 

The scattered light fringe pattern shows the in- 
tensities of scattered light along different directions 

13 

at different points on the optical path, Swinscn et al . 

At each point find the directions of observation °<i y <*2, ^3 
etc. along which S-j =0. Find the value of S,j along direc- 
tions ^ 2 , ** 3 e "tc. inclined at an angle of 45 ° to the 
directions <*1 <*2 ,*3 the corresponding points on the 

optical path. Let the value of S-j be nonzero corresponding . 
to the direction *< n say. Then the direction is along 
one of the secondary principal -stress directions. A laser 
beam of a very small diameter will be an ideal source for 
obtaining the required scattered light fringe pattern. The 
of foots Of varying lengths of (scatt.rod or transmitted) 
light path can bo eliminated in the following way. For this 
take two scattered light fringe photographs one with right 
and the other with left circularly polarised light at the 
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entrance. At any point and along a particular direction, 
the intensities in the two patterns are I 1 and 1 ^, then the 
required intensity I (in which the effect of varying light 
paths is eliminated) in the first pattern is 


h +1 2 


Because , 


I. = K a 2 
1 X 


I^i = K a £ 


I^+I 2 = K 

2 h x i 
= — = 

x K +I 2 


After knowing the value of 0 (one of the secondary 
principal -stress orientation) at each point, the value of p* 
can be calculated by using the relation 

P . -jl 

S 3 Sin 20 


The value of at any point can be found in the following 

way. 

At a point on the optical path and along a particular set of 

reference axes X and X, 

s i = 

i 1 - 1 2 

so s 1 



, ^ o 

where I. and I ar ~ -t-v, 

1 ± 2 a ^e the intensities of 

when the incident • t 8 scattei *ed light 

xgnt: is right and i n f+ 

respectively. circularly polarized 

- ind the value of d n i^ 

X direction. Thi s , 1 a ' dlreotl ' : ® ^5° to the 

lowing « and S th “ ^ *° S 2 ^ th ° Section. 

I 2* Lne v alue of s 

relation 3 J -°wnd from the 

s 3 = 0 - si 2 - s |)i 
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